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Abstract 


A detailed Monte Carlo calculation of the phase diagram of bosonic massdeformedIKKT Yang- 
Mills matrix models in three and six dimensions with quartic mass deformations is given. 
Background emergent fuzzy geometries in two and four dimensions are observed with a 
uctuation given by a noncommutative U(l) gauge theory very weakly coupled to normal scalar 
fields. 

The geometry, which is determined dynamically, is given by the fuzzy spheres 5^ and 5^ X 5^ 
respectively. The three and six matrix models are efectively 

in the same universality class. For example, in two dimensions the geometry is completely 
stable, whereas in four dimensions the geometry is stable only in the limit M —> oo, where M is 
the mass of the normal uctuations. The behavior of the eigenvalue distribution in the two 
theories is also different. 

We also sketch how we can obtain a stable fuzzy four-sphere 5^ X in the large N limit for 
all values of M as well as models of topology change in which the transition between spheres of 
different dimensions is observed. 

The stable fuzzy spheres in two and four dimensions act precisely as regulators which is the 
original goal of fuzzy geometry and fuzzy physics. Fuzzy physics and fuzzy theory on these 
spaces are briey discussed. 

Key words : 

Emergent geometry, Matrix Models, Fuzzy, Gauge Theory, Mulitrace, Monte Carlo, hybrid 
Monte Carlo. 



Resume 


Un calcul de Monte Carlo detaille pour le diagramme de phase des modeles matriciels de type 
Yang-Mills IKKT bosonique deformee en masse dans trois et six dimensions avec des 
deformations de masse quartique est aperues. Les geometries floues emergentes de fond dans 
deux et quatre dimensions sont observees avec une fluctuation donne par une theorie de gauge 
non commutative U(l) tres faiblement couplee aux champs scalaires normaux. La geometrie, 
determinee dynamiquement, est donnee par les spheres floues 5^ et 5^ X . 

Les trois et six modeles matriciels sont effectivement dans la meme classe 

d'universalite. Par exemple, en deux dimensions, la geometrie est completement stable, alors 

que dans quatre dimensions, la geometrie est stable seulement pour la limite M -» oo , ou M est 

la masse des fluctuations normales. Le comportement de la distribution des valeurs propres dans 

les deux theories est egalement different. Nous decrivons egalement comment nous pouvons 

obtenir une 5^ X 5^ floues a quatre spheres stable dans la grande limite N pour toutes les 

valeurs de M ainsi que des modeles de changement de topologie pour lesquels la transition 

entre les spheres de differentes dimensions est observee. Les spheres floues stables 

en deux et quatre dimensions agissent precisement comme des regulateurs qui est l'objectif 

initial de la geometrie floue (fuzzy geometry) et de la physique floue (fuzzy physics). 

Mots cles : 

La geometrie emergente, Modeles matriciels, la physique floue, Theorie de jauge, Multitrace , 
Monte Carlo, Hybride Monte Carlo. 





^cLail 4 Vi^ll-sjCCa (jjjjjjii ^ Wl<a\ 4L-*aLa ^^lijj-a ClAjLuiaj lla3 d-a*J! !c& ^^3 

C .ll m.Vlg 1' 4_^j.}]! ^ja 4_}lifa Clll& jdl) Clad 4 jlui j 4j!)Ij ^^3 J>lixs £-jl-l IKKT Cilljj'U.'-i.all 

4j^)Jaij ol}-»».a Cl'1 jISj ^a 1 g'i^-y^^La CllCi cl*ji 4juji j (jjCaj (j-a 4 .u«Ijl\I 4 >>».al»\l 4_}jLCa&i! 

.4_u»lJaj 4 .iaLm (Jjij-v ) 4 q i» >>i 43j-aJ 4jjjLa XJ( 1) 4 iLcCi jJC-^JjLla-a 
Sj\j X. Spj J Spj 4 Clll^Sllj oUaa-a 1 jfLj-alj,p o3,W.a\l 4 m.Vig \l 

^alLuia a-LCaill ^^3 ? JllLa .4_jj^a]l (_y&j l^J Lll*3 Claj'lH 4 _LuiI.1ui j 4_ll}Ij ClllS^aj-axAl! ^cLaj 

.4 i.al}-»ill dlLlilll 4 ^/[ (J j j^l d?j O A -laS3 ^piuus a-LCaaJI 4 jljjI LailJ LaLal) 

C!j£ Jj. 1\ 4jqj*s dlc£ .LCaji <■ q\'i-v.a (jCL^iai]! ^ 4jjlc]l ^oa]l ClAjujjj djiuj 

4ja^.Xs 4iUCaVlj M {£ <-£ idd Q* . M (J «CCblg-iJl X S^j 4-Cxalc. 

4 q\'i-y.a\l Clllj^]! (jjj dlVj~vMI (jjSj \g \-0 qa ( ^jk\ _JJ*j ^cLaj 

4 m.'ng 1] CllVc*-aS 43Cj (J-a*J Claji 4jujI j (jjCaj ^ o^>a!Lai«a]l 4jCa-al*J! Clll^)£]l 

.iLCsusliJI p.Ll jjsll J 4 xC-»<q\ »\1 p-LCail! 


4_i^llLdl ‘*M aKII 


j}j){£i ^lij^aJI ( ^lijj-a]l 4j^)lal ? 4 >>»,al»\l e-Ll^ia]! ; Ci)l3j) l U.'~».a\l ^iLaj ? 4 nlilill 4 m.Vig \l 



Emergent Geometry and Gauge 
Theory in 4 Dimensions and The 
Noncommutative Torus 


by 

KHALED RAMDA0 

Physics Department , Badji Moukhtar University 
DISSERTATION 


Submitted in partial fulfillment of the requirements for the 
degree of Doctor of Physics 

ADVISOR 

Professor Badis Ydri 


1 Email: khaled .ramda@univ-annaba. org 



Abstract 

A detailed Monte Carlo calculation of the phase diagram of bosonic mass- 
deformed IKKT Yang-Mills matrix models in three and six dimensions with 
quartic mass deformations is given. Background emergent fuzzy geometries 
in two and four dimensions are observed with a fluctuation given by a non- 
commutative U( 1) gauge theory very weakly coupled to normal scalar fields. 
The geometry, which is determined dynamically, is given by the fuzzy spheres 
S 2 n and x respectively. The three and six matrix models are effectively 
in the same universality class. For example, in two dimensions the geometry 
is completely stable, whereas in four dimensions the geometry is stable only 
in the limit M —y oo, where M is the mass of the normal fluctuations. The 
behavior of the eigenvalue distribution in the two theories is also different. 
We also sketch how we can obtain a stable fuzzy four-sphere S 2 N x S 2 N in 
the large N limit for all values of M as well as models of topology change 
in which the transition between spheres of different dimensions is observed. 
The stable fuzzy spheres in two and four dimensions act precisely as regula¬ 
tors which is the original goal of fuzzy geometry and fuzzy physics. Fuzzy 
physics and fuzzy held theory on these spaces are briefly discussed. 
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Chapter 1 

Introduction and Summary 


Understanding how our universe began is one of the most fundamental 
themes in theoretical physics, Gauge theories provide the best known descrip¬ 
tion of the fundamental forces in nature. At very short distances however, 
physics is not known, and it seems unlikely that spacetime is a perfect con¬ 
tinuum down to arbitrarily small scales. 

Gauge theory is one of the ways to try understand this dilemma, and pre¬ 
cisely Non-commutative gauge theory, Since its emergence in the 80’s 0 0 
10 0 . noncommutative geometry has helped to reveal deep mathemat¬ 
ical relationships between ordinary geometry and other structures, among 
them differential algebras and normed algebras. In particular, noncommu¬ 
tative geometry has shed new lights on gauge theories. Indeed, a theory 
of connections can be defined in great generality using the noncommuta¬ 
tive language of associative algebra, Different approaches have been pro¬ 
posed to study noncommutative spaces. The theory of spectral triples, de¬ 
veloped by Connes, emphasizes the metric structure i E 0- On the other 
hand, many noncommutative spaces are studied through differential struc¬ 
tures [8] [ 9 ] 10 11 [121. However, all the noncommutative gauge field 


theories studied so far use the same building blocks, even when there are 
defined through different approaches. 

A gauge interaction is an implementation of the principle that the theory 
should be invariant under some local symmetry. In particle physics, these 
local symmetries take the form of functions g : A4 —> G on the space-time 
M. with values in a structure group G. Electromagnetism is associated to the 
group G=U(1), the electroweak theory by Glashow, Weinberg and Salarn uses 
the group G = 1/(1) x SU{ 2) and chromodynamics relies on G = SU(3 ) fl31. 
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Noncommutative geometry being an extension of differential geometry, 
it naturally generalizes this theory of fiber bundles and connections. What 
is astonishing is that this generalization is very elegant, very powerful and 
very effective, not only from a mathematical point of view, but also in its 
applications to physics. 

To describe noncommutative spaces, the noncommutative geometry is 
now investigated by many authors and using this framework one can even 
consider the differential geometry of singular spaces like, for example, a 2- 
point space which has been shown to provide a geometrical interpretation of 
the Higgs mechanism |6| 14]. In general, such noncommutative spaces can 
be obtained by quantizing a given space with its Poisson structure. Further¬ 
more, if the original space is compact one obtains a finite dimensional matrix 
algebra as a quantized algebra of functions over this space. 

In the matrix model, matter and even spacetime are dynamically emerged 
out of matrices 15 16 . Spacetime coordinates are represented by matri¬ 

ces and therefore noncommutative geometry appears naturally. The idea 
of the noncommutative geometry is to modify the microscopic structure of 
the spacetime. This modification is implemented by replacing fields on the 
spacetime by matrices. It was shown [l7| 18 19 20 that noncommutative 


Yang-Mills theories in a flat background are obtained by expanding the ma¬ 
trix model around a flat noncommutative background. The noncommutative 
background is a D-brane-like background which is a solution of the equation 
of motion and preserves a part of supersymmetry. Various properties of non¬ 
commutative Yang-Mills have been studied from the matrix model point of 
In string theory, it is discussed that the world volume theory on 


view 


21 


D-branes is described by noncommutative Yang-Mills theory 22 


A different kind of noncommutative backgrounds, a noncommutative sphere, 
or a fuzzy sphere is also studied in many contexts. In |23j| it is discussed in the 
framework of matrix regularization of a membrane. In the light-cone gauge, 
they gave a map between functions on spherical membrane and hermitian 
matrices. In BFSS matrix model 


considered in 25 26 27 


24 membranes of spherical topology are 


in string theory context is discussed in 28 29 


A noncommutative gauge theory on a fuzzy sphere 

The approach to construct 


a gauge theory on the fuzzy sphere were pursued in 
The fuzzy sphere 


30 


30 31 32 33 34 35 


can be constructed by introducing a cut off param¬ 
eter N for angular momentum of the spherical harmonics. The number of 


independent functions is J^ 0 (2/ + 1) = (Y + 1) . Therefore, we can replace 
the functions by (N + 1) x (N + 1) hermitian matrices on the fuzzy sphere. 
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Thus, the algebra on the fuzzy sphere becomes noncommutative. 

One of the principal goals of the study of field theories on fuzzy spaces is 
to develop an alternative non-perturbative technique 36 . To date, this new 


approach in the case of four dimensional field theories has been limited to 

and S 4 


studies of Euclidean field theory on S 2 x S 2 37 , CP 3 38 


39. A 4 


is really a squashed CP 3 and includes many unwanted massive Kaluza-Klein 
type modes. Even S 2 x S 2 is not ideal since it has curvature effects that drop 
off as power corrections rather than exponentially as in the case of toroidal 
geometries 


40 


The fuzzy approach does, however, have the advantage of preserving con¬ 
tinuous symmetries such as the SU (2) symmetry of a round S 2 and does not 
suffer from fermion doubling |4l] 

In fact, this work divided into two main themes, first is Noncommutative 
scalar phi-four theory, and the other is Noncommutative gauge theory. In 
both subjects we focused on phases diagrams. 


• The first theme is phase diagrams of the multitrace quartic matrix 
models of noncommutative $ 4 , this theory is enjoys three stable phases: 
(1) disordered (symmetric, one-cut, disk) phase, (2) uniform ordered 
(Ising, broken, asymmetric one-cut) phase and (3) non-uniform ordered 
(matrix, stripe, twocut, annulus) phase. This picture is expected to 
hold for noncommutative/fuzzy phi-four theory in any dimension, and 
the three phases are all stable and are expected to meet at a triple 
point. The triple point is identified as a termination point of the one- 
cut-totwo-cut transition line and is located at (6, c) = (1.55, 0.4) which 
compares favorably with previous Monte Carlo estimate. 

• The seconde theme, we study a six matrix model with global SO( 3) x 
SO (3) symmetry containing at most quartic powers of the matrices. 
This theory exhibits a phase transition from a geometrical phase at 
low temperature to a Yang-Mills matrix phase with no background 
geometrical structure at high temperature. This is an exotic phase 
transition in the same universality class as the three matrix model 
but with important differences. The geometrical phase is determined 
dynamically, as the system cools, and is given by a fuzzy spherical 
background S 2 N x Sjy, with an Abelian gauge field which is very weakly 
coupled to two normal scalar fields. 
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Chapter 2 

Noncommutative Field theory 

2.1 The caconical case 

Nocncommutative gauge theory in the canonical case,where the commu¬ 
tator of two coordinates is a constant. We will start with the most commonly 
used ^-product for the caconical case. 

2.1.1 The Moyal-Weyl ^-product 

In the canonical case, the noncommutative coordinates fulfill commuta¬ 
tive relations 


[£*, x j ] = iO lj 


( 2 . 1 ) 


with the constant noncommutativity parameter 6 G R. The noncommutative 
algebra generated by the noncommutative coordinates can represented on 
the space of functions on R n by introducing a noncommutative product, the 
Moyal-Weyl ^-product 


42 43 


/ * 9 = ®g = fg + % -B li diSd,g + 0(2) 


( 2 . 2 ) 
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2.1. THE CACONICAL CASE 


with m.(f <E> g) = fg and d i = The product is associative,as 

(f*g)*h = m.e^ kldk ® dl {mJ diidi ^f®g)®h (2.3) 

= m m e h ekl< ^ d k® 1 ® d i+ l ® d k®di)di®djm f ® g & fr, 

= m.m.e5 fly(ai ® 1 ®^ +ai ® s >®<g> g 0 h 
= m.e ^ 6 ' 3 di ® dj (f ® [rn.e^ 9Udk ® dl g ® h ) 

= f -k (h-k h) 

and obviously reproduces (1.1). Furthermore, as 9 is antisymmetric, usual 
complex conjugation is still an involution 

J^~g = m .e-^ d ^J®g = gkf (2.4) 

and integration has the trace property 

J d n xf k g = j d n xgkf (2.5) 

if the function f and g vanish suffciently fast at infinity because f k g has to 
be integrable in the firs place 

Differentiation on this space is an inner operation, i.e. we have 

i() ,n '(K = [.x •] ( 2 . 6 ) 

2.1.2 Commutative gauge theory 

A non-abelian gauge theory is based on a Lie groupe with Lie algebra 

[T a ,T b } = iffT c (2.7) 

Matter fields transform under a Lie algebra valued infinitesimal parameter 

A = A a T a (2.8) 

in the fundamental represntation as 

Sx'i/j = iX'ip (2.9) 


It follws that 


— 5 \ 5 ^) = 


( 2 . 10 ) 
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2.1. THE CACONICAL CASE 

The commutator of two consecutive infinitesimal gauge transformation closes 
into an infinitesimal gauge transformation. 

As differentiation isn’t a covariant operation, a Lie algebra valued gauge 
potential a* = a ia T a is introduced with the transformation property 

5 x ai = <%A + i[A, a*]. (2.11) 

With this the covariant derivative of a field is 

'Di'ip — ditj} — iaitjj. ( 2 . 12 ) 

The field strength of the gauge potential is defined to be the commutator of 
two covariant derivatives 

fij = i\Pii T>j] = didj - djdi - i[di, dj] (2.13) 

For nonabelian gauge theory, f VJ is not invariant under gauge transformation, 
but rather transformation covariantly, i.e. 

d x f = i[\,f}. (2.14) 

The same is true for the f V) f 13 (Lagrangian density).In order to get a gauge 
invariant action, we have to use the trace over the representation of the gauge 
fields. As the trace is cyclic, the commutator with the gauge parameter 
vanishes and the action 

S= j dx"tvf l3 r (2.15) 

becomes invariant. 

2.1.3 Noncommutative gauge theory 

To do noncommutative gauge theory in the ^-product approach, we can 
simply mimic the commutative consrtuction, replacing the ordinary product 
with the ^-product. 

Field should now transform as 

<5 A 'L = (2.16) 

The commutator of two such gauge transformations should again be a gauge 
transformation, i.e we want 
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2.1. THE CACONICAL CASE 


(y A -y A )^-%: A] $, (2.i7) 

which is only possible for gauge groups U (N), as for A = A a T a and 5 = 5 a T a 
the commutator 

[2, A] = ![=„ * A b ]{T\T b } + t{2„ * A b }[T\T b ] (2.18) 

As coordinates do not transform under gauge transformations, multiplication 
from the left with coordinates no longer is a covariant operation, i.e. 

5\(xi * T) = Xi * A * T 7 ^ A * Xi * T (2.19) 

This is very much like the situation in commutative gauge theory, where 
acting with a derivative from the left isn’t a covariant operation. Following 
the procedure there, we introduce covariant coordinates X 1 by adding a gauge 
field Ai as 


X i = x* + 0 ij Aj (2.20) 

To make the X 1 covariant, i.e. 5\X l = z[A * A"*], the gauge field has to 
transform as 

Sa^Aj) = -ifr * A] + i [A * e ij Aj] ( 2 . 21 ) 

and therefore 

S A Ai = diA + i[A* Ai] (2.22) 

in exact analogy to the commutative case. The commutator with the co¬ 
ordinate produces the derivative on the gauge parameter, as [x\ f} = iQ 1 ' 3 f. 
More generally we can introduce a covariantizer T> that applied to a function 
f renders it covariant [44] 

S A {V(f))=i[\*V(f)] (2.23) 

We can now go on to formulate noncommutative gauge theory much in the 
same way as we formulated commutative gauge theory. 

The covariant derivative T>, can be introduced as 

= dyk - iAi * T, (2.24) 
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2.2. THE GENERAL FORMALISM 


the filed strength F\j as 

Fij = i[Di * Vj] = diAj - djAi - i[Ai * Aj\. (2.25) 

The relation to the covariant coordinates subsists at this level with 

-i([ X i * X j ] - i9 ij ) = 9 ik 6 jl F kl (2.26) 

In noncommutative gauge theory,the field strength F is not gauge invari¬ 
ant, even for gauge group 1/(1). It rather transforms covariantly under gauge 
transformations, i.e 

= -i [A * F llv *F liV ]. (2.27) 

Just inserting a trace over the representation of the gauge group no longer 
guarantees gauge invariance. To get gauge invariant expressions, we have to 
use the trace property of the integral. If we set the action for non- commu¬ 
tative gauge theory as 

S= f d n xtrF IUJ * F IW , (2.28) 


2.2 The general formalism 


2.2.1 Seiberg-Witten gauge theory 


Noncommutative transformation properties are determined by the trans¬ 
formation properties of the commutative fields they depend on. Therefore 


the fields again transform as 45 


S a ^^[a] = iA a [a} * \ER,[a], 

leading to the same consistency condition for the gauge parameter 

i(d a A/3 (^Aq,) T [A q , A^] ^A_i[a )( 3] 

The transformation law for the covariantizer is now 


(2.29) 

(2.30) 


5 a (D[a](f)) = i[A a [a\ * D[a\(f)]. (2.31) 

The Seiberg-Witten-map can be easily extended to the derivations 5\ of the 
*-product.The noncommutative covariant derivation Dx[o] can be written 
with the help of a noncommutative gauge potential Ax [a] now depending 
both on the commutative gauge potential a* and the Poisson vector field X 


8 a A x [a\ = JxA Q [a] + i[A a [a\ * A x [a}}. (2.32) 
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2.2. THE GENERAL FORMALISM 


2.2.2 Commutative actions with the frame formalism 

we recall some aspects of classical differential geometry. Suppose we are 
working on a n-dimensional manifold M with metric r/ /iiy .Then there are lo¬ 
cally n derivatives d p which form a basis of the tangent space of the manifold. 
We can always make a local basis transformation to frame 

e„ = (2.33) 

with e£e“ = 5% , where the metric is constant 

Vab = « (2.34) 

Since forms are dual to vector fields, they may be evaluated on the frame. 
For the gauge field we get 


o(e tt ), (2.35) 

leading to the covariant derivate 

A/, = {Dijf){e a ) = e a ^> - ia a A (2.36) 

The field strength becomes 

fab i\D a , Db\ iD(fc aj Cfe]) ®([^aj £&]) ®b] (2.37) 

Locally this means that 

a a = <a M , D a f) = e^D,, and f ab = (2.38) 

Using these defnitions, the action for gauge theory on a curved manifold can 
be written in the two diferent bases as 

S = ~\f d n x^ V ab V cd facfbd = ~\J d n x^g^f, p f ua (2.39) 

where 

= \J det (g^) = \Jdet(e“elriab) = det e“ (2.40) 
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2.2. THE GENERAL FORMALISM 

2.2.3 Gauge theory on curved noncommutative space- 
time 

In order to formulate gauge theory on a curved noncommutative space- 
time, we need a frame e a and a Poisson structure {., .} p = ir^d^ A d v that 
are compatible with each other. Compatibility means that the frame e a 
commutes with the Poisson structure {., .} p , i.e. 

^a{f 1 9 }p {l/) 3}|) t {/) Lslpi (2-41) 

We can define a covariant derivative of a field by using a derivation 8x 

Dx^i’ = fix y by — iAx * V IV (2.42) 

With this, a field strength could be defined as 

—i I'x.y = [D x * D y ) — D[ X ,y]* (2.43) 

The noncommutative covariant derivative (2.14) and field strength (2.15) 
evaluated on the frame e a then read 

D a § = D ea § = <5 ea $ - iA ea * $ (2.44) 

-iF ab = - iF Cafib = [D ea * D e J - Di ea!eb]ir . (2.45) 

The field strength will transform covariantly under gauge transformations, 
i.e. we have 

<$a = *[A?F]. (2.46) 

To make the action gauge invariant, the integral has to have the trace prop¬ 
erty, a noncommuative gauge action 

S = -\J d n xQ V ab V cd FacF bd (2.47) 

with = y/~g + 0(1 ) 

that goes in the commutative limit 

S ^~\j 9 IJJJ 9 pa .UUa (2.48) 
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2.2. THE GENERAL FORMALISM 


Scalars 

For the noncommutative version of a scalar Lagrangian 

r) ab Dj>D b <\) + m 2 #, (2.49) 

we also need an involution of the ^-product, i.e. 

( f*g)=9*f • (2-50) 

To make the NC Lagrangian invariant under NC gauge transformations, the 
NC gauge parameter A and the NC gauge held Ax have to be invariant under 
this involution to get 

= (A * (ft) = 0*A = < f>*A (2-51) 

and 

(Ax * 4>) = 0* A x = (f>* Ax (2.52) 

Putting everything together, we therefore end up with an action 

S = j d n xQ.(— ^rj ab rj cd F ac * F bd + r] ab D a $ * D b $ - m 2 $ a 4>) (2.53) 

that is invariant under noncommutative gauge transformations 

S A S = 0 (2.54) 

and reduces in the commutative limit 

S J d n x y /g(-^g lxv g prT f lip f v<T + - m 2 ^) (2.55) 

to scalar electrodynamics on a curved manifold 

Spinors 

The commutative spinor action can be written as 

s spinor = \ j d 2 Xy/g^/i'fe^(d ll ~ iA H + m ( 2 -56) 

Usual gamma-matrices { 7 a , 7 6 } = 2// ab and 7 ^ = 7 a e((, we get = 

2g pu . The noncommutative version (2.28) is easily constructed, and we get 

S spinor ^ J d 2 xO,^ij a (S ea — iA ea + m) * 4/ (2-57) 

with e a = e'gdfj, which is invariant under NC gauge transformations 
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2.3. NONCOMMUTATIVE SCALAR FIELD THEORY 


2.3 Noncommutative scalar field theory 

We are able to define a scalar field theory on this geometry. At this point 
we make a change of notation and introduce the short-hand notation Weyl 
operators W[f] —> f 

2.3.1 Noncommutative scalar action 

We start with the ation of an Euclidean commutative A0 4 theory 

S[(j)} = I d d x(^d fl (p(x)d^x) + ^-0 1 2 (a) + (2.58) 

where 0 is a real valued scalar field and d is the dimension of spacetime. 

To transform an ordinary scalar field theory to a noncommutative field 
theory we can use the Weyl quantization via Hermitian operators (j) 

The quantum field theory written in terms of Weyl operators <f>, correspond¬ 
ing to a real scalar field 4>{x) on [R d 

Z — J d0exp ^ — S[(j>]^ (2.59) 

5[01 = + W) (2.60) 

The measure d<p is here the ordinary path integral measure for scalar fields 
V<p 

This theory may be formulated in coordinate space by applying the map 

f{x) = Tr(W[f] A{x)) (2.61) 


W[f] = / d d xf(x)A(x) with A(x) = 


d d k 

(2tt) c 


^ikpXng—ikpXn 


(2.62) 


where A (a;) is a Hermitian operator that can be understood as a mixed basis 
for operators of fields, and using bL[/]bL[g] = W[f * g ], 


A[0] = / d d x 


1 Tfl A 

2 i d M x )) 2 + * 4>{x ) * 4>{x) * <t>{x) 


.(2.63) 
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2.3. NONCOMMUTATIVE SCALAR FIELD THEORY 


the kinetic term and the mass term do not contain the star-product, because 
of this property f d d fi(x) * f 2 (x) = f d d f 1 (x)f 2 (x) . 

This implies that the free propagator of the quantum held theory also 
remains unchanged, and only the interaction part, i.e., the vertex of Feynman 
rules, gains extra contributions from the noncommutativity. (£g) 

In particular, the commutative vertex function iX gains an extra phase 
factor of the form 


V(ki) = J2 e ~^ kj (2.64) 

i<j 

and ki are the momenta flowing into the vertex. V ( ki ) is not invariant under 
arbitrary permutations of the momenta but only under cyclic permutations, 
so one has to keep track of the order in which propagators are connected 
to the vertices of Feynman diagrams. We can be drawn on a plane without 
intersecting prooagators,’planar’. Namely, we may replace every line in a 
planar Feynman diagram by a double line.which has only non-intersecting 
solid lines and loops. [?] (seeFig.(l)). due to the momentum conservation 
at the vertices and the planarity, we may label the lines of the double line 
notation by momenta which correspond to the original momenta via the 
relation ki = — l i2 . Accordingly,when /q i = 1,4, are the incoming 

momenta for a vertex in cyclic order, the phase factor (3.7) becomes 

g—f E£=i kjAiij+i (2.65) 

where each of the expressions l tj A l tj +\ corresponds to one of the incoming 
propagators. The over-all phase factor for a diagram is then the product 
of the phase factors corresponding to each of the vertices of the diagram. 
Therefor, by the expression (3.8), we End that for a planar diagram the 
factors corresponding to the internal propagators cancel out, since they con¬ 
tribute opposite ±/j. A U j+ 1 to the exponent of the over-all phase factor. 
Consequently, we are left with the over-all phase factor 

V( Pi ) = e -^<i pihpj (2.66) 

Pi are the momenta associated to the external lines of the original single- 
line diagram in the cyclic order [?]. It immediately follows that the UV- 
divergencies of the commutative quantum scalar held theory, which arise 
from the integrals over the internal momenta of Feynman diagrams, are also 
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2.3. NONCOMMUTATIVE SCALAR FIELD THEORY 


present in the planar diagrams of the noncommutative theory, and therefore, 
due to this example, the noncommutativity of soacetime does not seem to 
help to naturally regularize the divergencies of quantum held theory 


2.3.2 UV/IR mixing 

We consider the one-particle-irreduciblc two point function T of the non¬ 
commutative A0 4 scalar held theory. 

OO 

rfp) =< >= Y. A ” r ‘(p)- P- 67 ) 

n =0 


At lowest order the two-point function is given by T°(p) = p 2 + m 2 . 

Where the hrst term is the planar contribution and the second term the 
nonplanar one. In fact, since we have only one external momentum in the 
case of two point function, the planar phase factor (3.9) equals unity, and thus 
the planar diagrams give exactly the same correction as in the commutaive 


case. 


The one loop contribution splits topologically into two parts, one planar 


and one non-planar diagram 46 


T 1 

p 


T 1 

np 


A 


d A k 


3(2vr) 4 
A 


k 2 + m? 

d A k 


} ikf\p 


6(27t) 4 J k' 2 + m? 


( 2 . 68 ) 

(2.69) 


In Refs [?] 46 it is shown that the contribution of planar diagram to non- 


commutattive perturbation theory is propotional to the commutative case 
(hg planar and non-planar) 


Regularizing the momentum integrals at the energy scale A we hnd 46 

A 


T 1 = 

p 


T 1 = 

np 


487T 2 

A 

967T 2 


A 2 — m 2 In ( + 

rri 2 J 

A 2 
A e// 


A e// -m 2 In 




+ ... 


Here we introduced the effective cutoff A e ff given by 

1 


A eff ~ 


1/A 2 + p o p 


pop := -p^Q^y > 0 


(2.70) 

(2.71) 

(2.72) 
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k 



+ 



Figure 2.1: Planar and nonplanar loop diagrams, respectively. 
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2.3. NONCOMMUTATIVE SCALAR FIELD THEORY 


When we take the UV-limit A —» oo of the internal momentum, we find 


-pi A—>00 

^ np 


x 

96n 2 



m 2 In 


m 2 (p op) 2 



(2.73) 


This expression diverges at the low energy limit p —> 0 (IR limit), The UV 
limit does not commute with the IR limit. 

At small momenta or small non-commutativity parameter the two-point 
function reads 


T(p) ~ p 2 + m 2 + 3AT^(0) + 0(A 2 ) (2.74) 


Taking now the UV limit leads to the standard mass renormalization of the 
A0 4 theory. Taking these limits vice versa, the effective cutoff is given by 



(2.75) 


and A e ff diverges and therefore also T^Qo) either in the limit 9 —> 0 or in 

the infrared limit when the incoming momentum p —> 0 

The planar one loop contribution of T(p) by defining the renormalized mass 

through 


M 2 ef f = m 2 + 2AT^(0) (2.76) 

Removing the cutoff while keeping M 2 ^ fixed, then leads to a finite T(p) 
for finite incoming momenta p. For zero momentum T(p) diverges and the 
divergence at one loop is given by 

r (p) = p 2+ M 2 ff + + subleading terms, with £ = (2.77) 

This is not too surprising, since the nonplanar phase factor e tkAp in (3.12), 
which dampens the singularity of the momentum integral, approaches unity 
as p —* 0 

This exotic mixing of the high energy (UV) and low energy (IR) scales in non- 
commutative theories, which does not have a counter-part in commutative 
theories, is called LTV/IR mixing 
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2.3. NONCOMMUTATIVE SCALAR FIELD THEORY 


2.3.3 Phase structure of non—commutative A0 4 

The UV/IR mixing is one of the most intersting properties of non—commutative 
field theory and has no connterpat in the commutative case. 

Gusber and Sondhi studied the phase diagram of 4d A0 4 theory 
based on an action of the Brazovskiia from 


48 


47 


We try summarize their resuls. At small non—commutativity parameter 
6 they obtained an Ising type (second order) phase to an uniformly ordered 
phase with < 0 > 7 ^ 0 . 

At large 6 , Leads ao an ordered phase. In this phase < 4> > varies in space, 
which involve some non-uniform patterns like stripes. 

According to Ref 


47 


the pahase diagram in the — A plane is then given 
by Figure (|2.3l) , where A is a momentum cut-off and the A is the coupling. 


I 11 another approach renormalization group technique were used to study 
the phase diagram of the A0 4 model 49 


Chen and Wu obtained in d — 4 - 


a new IR stable fixed point, and therefore a striped phase exists. In contrast 


to the resuls in Ref 47 , this implies that in d — 4 there is no stiped phase. 
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Figure 2.2: he phase diagram for the Ising scalar held theory, in d > 1. 
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2.4. NONCOMMUTATIVE GAUGE THEORY 


2.4 Noncommutative gauge theory 


2.4.1 Star-gauge invariant action 

To define a YangMills theory on a noncommutative plane we have to 


generalize the map (2.62). Let AJx) be a Hermitian gauge field on [R , 


which corresponds to the unitary gauge group U(n). We can introduce the 
Weyl operators corresponding to A^x) by taking the trace of the tensor 
product of A(Y) and the gauge field [50 


(2.78) 


Afj,= d d xA(x) <g) A^x) 


the derivative of Weyl operators is equal to the Weyl operator of the usual 


derivative of the functions 50 


[S„,w[!\\ = / d d xd„f(x)A(x) = w[a„f\ 


(2.79) 


Based on this equation a noncommutative version of the Yang-Mills action 
can ce defined 

S[A] = --^Tr trAr ^[<9 m , Ay] - [d v ,Aj\ (2.80) 

where the term in brackets is the operator analog of the field strength tensor. 
Here Tr is the operator is given by an integration over space-time 

Tr W[f] = I d d x f(x) (2.81) 

and tr^ denotes the trace in color space. This action is invariant under 
transformations of the form 

A^UA^ -iU[d^,U], (2.82) 

where U as an arbitrary unitary element of the algebra of matrix valued 
operators, i. e. 

UU* = U'U = 1 ® t n (2.83) 

The symbol 1 is here the identity on the ordinary Weyl algebra and l n is a 
n x n unit matrix. 
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2.4. NONCOMMUTATIVE GAUGE THEORY 


To set up the action in coordinate space we can construct an inverse map 


of (2.78). The Yang-Mills action in coordinate space the reads 


1 


S[A\ = ~— J d d xtr N (F^(x) kF^x)), 


(2.84) 


where we introduced the noncommutative field strength tensor F„ u given by 


Fp,(x) = dpA^x) - dvA^x) - i[A^{x),A u {x)\ 


(2.85) 


The index V indicates that the products in this commutator are star-products. 


According the (2.85) the simple gauge group £7(1) we have Yang-Mills type 


strucure. herefore there exist three and four point gauge interactions and 
noncommutative £7(1) theory is asymptotically free. 

The action ( |3.54[ ) under star-gauge transformation given by 

A^ix) —» £7(x) ★ A^ix) -k U(x)^ — iU(x) k d^U(x)\ 


UiW > w y^j, " x 

£7(x) is a star-unitary matrix field, 

U(x) k U{ X y = £7(x) t * £7(x) = 1 


( 2 . 86 ) 


(2.87) 


On the classical level we considered non-commutative U{n) theories which 
reduce to the ordinary £7(n) theories in the limit 9 —» 0. 

In Ref. 51 it was shown that the gauge groups SU(N), SO(N) cannot be 


realized on a flat non-commutative manifold, while it is possible for U(N) 

The interaction of the SU(N) gauge bosons with the £7(1) gauge boson 
plays an important role in the consistency check. In particular, the SU(N) 
theory by itself is not consistent, and when 0—^0 limit of the U(N) theory 
does not converge to the ordinary SU(N) x £7(1) commutative theory, even 
at the planar limit [521 

The U(n) group is closed under the starproduct; the product of two 
starunitary matrix fields is again starunitary. In contrast to £7 (n) the special 
unitary group SU(n) is not closed, since in general [50] 


det(£7) x det(H) ± det(£7 * H ). 


( 2 . 88 ) 


2.4.2 Gauge-invariant observables 

Let us say a few words about observables in noncommutative Yang-Mills 
gauge theories, since it is such a subtle subject, which clearly highlights the 
nonlocaliy of these theories, and is also deeply related to UV/IR mixing. 
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2.4. NONCOMMUTATIVE GAUGE THEORY 


We consider an arbitrary ariented smooth contour C C lZ d with smooth 
parametrization £(t) : [0,1 ]—> lZ d , and endpoints £(0) = 0, £(1) = v in lZ d . 
Introduce the noncommutative parallel transport operator, construct gauge 
invariant observables using open Wilson lines, which are nonlocal operators 
defined as 


53 54 55 


W(x;e) = J Pexp*| * 


dC Ai{x + £) 


>c 


(2.89) 


P is path ordering. The index V at the exponential function indicateds 
that in the expansion of this function the star-product has to be used. The 
operator W(x;£) is an n x n star-unitary matrix field and transforms under 


the star-gauge transformation (2.86) like (it is important property of an open 
Wilson line) 


W(x; 0 -)■ U(x) * W(x- 0 * U(x + vY (2.90) 

under finite noncommutative gauge transformations (54 

W{k-Y) ■= I d A tr[W(x;0] *e ik ' x (2-91) 


The transformation of W(fc;£) under a gauge transformation according to 


(2.90) is 


>V(fc;0-> d i tT[U{x)*W(x-,g)*rf(x + v)*e ik - x ] (2.92) 


In the traditional commutative use, gauge invariance would force us to 
close the contour C. But in the noncommutativr case, gauge transformation 
can affect translations of space-time (the translations can be arranged by 
(star-) multiplication with plane waves): 

U(x + v) = e ik ^ * U(x) * e~ ikpXp (2.93) 


where k = 9 .v is the total momentum of path C. With the definition of the 
noncommutative parallel transporter and equation (5.1) we can associate a 


star-gauge invariant observable with any arbitrary contour C v by 54 


0(C V ) 


trjv(W(: v-C v )*e ikpXli ). 


(2.94) 
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We notice here that UV/IR mixing manifests itself in the property that 
k —> oo as v —> oo 

Gauge-invariant operators, generalizing the standard local gauge theory 
operators in the commutative limit, are local in momentum space and are 


given by a Fourier-type transformation 54 


6(k) = / d d OU{x ; C k ) * e i 


(2.95) 


where O(x) is any local gauge invariatn operator of ordinary Yang-Mills the¬ 
ory. In the commutative limit 9 = 0,u = 0; there are no gauge-invariant 
quantities associated with open lines in ordinary Yang-Mills theory. In that 
case, the totarn momentum of a closed loop is unresricted, and we can rem- 
place e lkflXfl (the momentum eigenstate) by an arbitrary function f(x) paric- 
ular, taking f(x) = 8 d (x — a) recovers the standard gauge-invariant Wilson 
loops of Yang-Mills theory. But for 9 ^ 0, closed loops have 0 memntum k, 
and only e lkljXlJ = 1 is permitef in O(C) above there is no local star-gauge 
invariant dynamics, because evrything has to be smeared outby the Weyl 
operator tarce Tr ^ f d d x. Hence the gauge dynamics below the noncom¬ 
mutativity scale is quite different from the commutative case. [;53| 


2.4.3 Application 

Noncommutative QED and UV/IR Mixing 


Despite the growing understanding in formulating noncommutative gauge 
held theories, The problem of UV/IR mixing still shows up, as expected, 
when one calculates the higher order diagrams of noncommutative gauge 
held theories. 

The propagators of noncommutative gauge theories are again equal to 
those of their commutative counter-parts, but the vertex functions contain 
nontrivial phase factors, as in the case of noncommutative scalar held theo¬ 
ries, which ultimately lead to the UV/IR mixing 56 . For noncommutative 


QED one obtains the vertex functions for Feynman diagrams 57 


For the correction n((/ given by a massless fermion loop to the photon 


propagator one finds 58 


d 4 ~ 


-sin 2 (-p A k) 


*n^ = —4 g 2 


27r 4 (p — k) 2 p 2 


(2.96) 
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Using 


sin 2 (-p A k) — -[1 — cos (p A k)] 


(2.97) 


we can isolate the planar and non-planar contributions. The planar part gives 
then the usual logarithmically UV-divergent but renormalizable contribution, 
whereas the nonplanar part with the dampening phase factor cos(pAfc) gives 
the leading order term 


56 58 




k^k h 

y 


(2.98) 


at the IR-limit of the external momentum, which clearly diverges quadrat- 
ically as k —y 0. Therefore we again encounter the UV/IR mixing, where 
an IR-divergence of the external momentum arises from the UV-limit of the 
integral over the internal loop momentum. Similar IR-divergencies arise also 
from other higher order corrections to propagators and vertices 156 


The physical interpretation of terms like the (5.2) is very interesting. For 


small noncommutative momentum, the one-loop inverse propagator is given 
by 


58 




(*, 


k 2 3 






, k tJ ‘k t 

y 


(2.99) 


where K (fcg = k 2 +K' 2 ) represents the projection of the spatial momentum on 
the plane.From this one-loop inverse propagator we can read the dispersion 
relation for the two physical. Suppose k is along the 2direction so that K is in 
the 1-direction. Then the photon polarized in the direction perpendicular to 
K satisfies the same dispersion relation as a photon would in the commutative 
theory 


UV/IR Mixing in Noncommutative QED via Seiberg-Witten Map 


In the paper 59 on the connection between noncommutative geometry 


and String Theory, Seiberg and Witten introduced a mapping, which relates 
gauge field theories in noncommutative spacetime to ordinary commutative 
onesknown as the Seiberg — Witten map. This mapping has virtues, since 
some aspects of gauge theories, such as observables and gauge fixing, are 
more easily understood and dealt with in the commutative theories. 
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for example in [601 , that the UV/IR mixing is absent in the Seiberg- 
Witten formalism. However, we will find that this is presumably due to the 
expansion in the noncommutativity parameter matrix 9 in the 0-expanded 
Seiberg — Witten map. In the 0-exact Seiberg-Witten map for noncom- 
mutative QED the UV/IR mixing reappears, as we will demonstrate. This 
same argument was expressed by Schupp and You in 61 , where a noncom- 


mutative model with a gauge held coupled with a spinor held in the adjoint 
representation was considered. 

The adjoint representation of the gauge group, however, corresponds to 
a chargeless particle with an electric dipole moment proportional to 9 , and 
therefore in their model the interaction vanishes at the commutative limit 
9 —> 0. Accordingly, their model does not correspond to a noncommutative 
theory of electrically charged fermions, which should reduce (classically) to 
the connnu- tative QED in the commutative limit. 
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2.4. NONCOMMUTATIVE GAUGE THEORY 


Noncommutative QED via Seiberg-Witten Map 

We consider exclusively the gauge group [/*( 1). We want to express the 
action of noncommutative QED, 


Sncqed = I d 4 x T(z$ - m)T - T * A * T - * F l ‘ ,u 


( 2 . 100 ) 


in terms of the commutative fields up to the first order in A, so that we 
can calculate the photon propagator correction coming from the one-loop 
photon self-energy diagram. Denoting the noncommutative fields by hats 
and dropping the lower index from 6\. The gauge held 

A^ — + 0(A 2 ), 

and for the spinor held 


vp = vp _ -Q a P 
2 


A a (<9/34/) + ~(dpA a ) T 


+ 0(A 2 ), 


( 2 . 101 ) 

( 2 . 102 ) 


where we use the notation 


(/*i g)(x) ■= 


e 2 d\ A d 2 — 1 
|<9i A d 2 


f(xi)g(x 2 ) 


(2.103) 


X\=X2=X 


Since (/ * gy = g' * P for any functions (or matrices) / and g, we hnd that 


T = vj = - 

2 


(dpty A a + -T (dpA a ) 


+ C(A 2 ), (2.104) 


Substituting (2.101), (2.102) and ( |2.104 ) into the action (2.100), we hnd the 
fermion- photon interaction term to be, up to hrst order in A, 

(2.105) 


/hi) _ 


l 


(<9/3 v b) A a + -T (dpA a ) 


(i<j — m)\ 1 / 


— -0 a/3 4/(0 — m) 


A a *i (<9/34/) + ~(dpA a ) T 


For the corresponding vertex function we get the expression 

j pffclAfc2 1 

V^h, k 2 ) = -iYe^ k > - -(£ - k 2 n^ + y 2 ) - 

2 k\ A k 2 

where k x and k 2 are the incoming momenta of the outgoing and incoming 

fermions, respectively. 


(2.106) 
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Photon Self-energy and UV/IR Mixing 

Now, using the vertex function ( |2. 106 ), we find the first order fermion 
loop correction to the photon propagator given by the one-loop photon self¬ 
energy diagram in Fig. (4.2) to be 


n %(k) = 


d A P 

Wf 


X 


j sim( z p/\k) 
2 A k 


(P ~ -kYk p T pv e--^ k 


~ ip- -k) u k p T plx e* pAk 


+ - 


2 ' - p 
i sim 2 i\p A k ) 
4 (ipAfc) 2 


(P - \kr(p - 1 -kYk ll k„T l ” 


(2.107) 


where 


T^iKp) : = 


{p — kYp v + p^ip — k) u + [m 2 — ip — k).p}i 
[ip — k) 2 — m 2 } \p 2 — m 2 ] 


(2.108) 


which is the only term we get in the commutative case. Therefore, the first 


term in (2.107) is naturally understood to correspond to the planar part of 


the diagram, and in fact follows straightforwardly from the first terms of 


the vertex functions (2.106) as the phase factors cancel each other, in the 


same way as they do for the planar diagrams of a noncommutative scalar 
field theory. The other terms, on the other hand, clearly correspond to the 
nonplanar part with nontrivial phase factors that give rise to UV/IR mixing. 


Indeed, the second term in (2.106) can be shown to yield the leading order 


contribution (The third term leads also to similar IR-divergent terms) 


^ np ik) 


8 U L k v 4 k) l k v + W 

71-2 ^4 71-2 £,4 


(2.109) 


at the IR-limit of the external momentum. The first term in (2.109) is similar 


to (2.98) found in the naive formulation above, whereas the second term is 


gauge variant and should cancel, when all the second order contributions in 
the coupling constant are taken into account. Therefore we conclude that 
also gauge field theories defined via Seiberg — Witten map appear to fail to 
be renormalizable because of UV/IR mixing, which further shows that this 
is a generic property of noncommutative theories. 
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P -k 



Figure 2.3: One-loop photon self-energy diagram. 
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2.4.4 Reduced Models and Emergent Phenomena 


In this part we will work out the nonperturbative, construction definition 
of noncommutative Yang-Mills theory, this can be completely described in 
the language of matrix models (arising here as reduced models). This will also 
reveal some beautiful features of the vacuum structure of noncommutative 
gauge theories. 

Will remove derivative operators <9,- or d t from the noncommutative gauge 
theory action. There is no analog of this manipulation in ordinary Yang-Mills 


theory 53 


Let us introduce covariant coordinates: 


Ci = (e-% xi + Ai (2.110) 

Then Q —* UCiW under gauge transformations. We can be represented the 
adjoint actions 

(r%[F, -] (2.111) 

i.e. di are inner derivations of the algebra 1Zg. Then the entire noncommu¬ 
tative gauge theory can be rewritten in terms of the C u we may rewrite the 
covariant derivative as: 



(2.112) 

where d[ = <9* + ij& . Then using [d^'F] = 0, We compute: 


[Vi,/] = -i[4 /] 

Fij = f[Vi,Vj] 

= -i[c' i ,c' i ] + (r% 

(2.113) 

(2.114) 

and consequently, 


s™ = Y T r£([4c/+;(0«) 2 
g 

(2.115) 

Cj are element of the 7 Zg, so spacetime derivatives have completely dis- 
appeared in this rewriting of noncommutative Yang-Mills theory. 53 

Flat connections F tJ = 0 give: 

5> 

i 

<S» . 

1 

(2.116) 
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Ci are like the momentum operators idi. In particular, X 1 = Q^Cj formally 
represent noncommuting position operators in the ground state, wherein 
X 1 — x l ( Ai = 0). Then the noncommutative gauge degrees of freedom 
Ci are fluctuations around this canonical (Moyal) noncommutative space- 
time. More generally noncommutative spacetime [X’ l 1 X^\ = (X) are 

obtained as non-vacuum solutions of Yang-Mills equation of motion: 

[C i ,[C i ,C j }] = 0 (2.117) 


Thus (noncommutative) spacetime emerges as a dynamical effect in the ma¬ 
trix model. This is the essence of its relation to the so-called IKKT matrix 
model for the non-perturbative dynamics of type IIB superstrings, and also 
in the more recent models of emergent gravity which clarify the origin of 
gravity in noncommutative gauge theory. In this setting, gravity is related 
to quantum fluctuations Ci of spacetime at the Planck scale, while noncom¬ 
mutative held theory arises from held dependent fluctuations of spacetime 
geometry (determined via 0 VJ {x)). In particular, UV/IR mixing arises due to 
a non-renormalizable gravitational sector in the IR with G 


A 53 


This large N matrix model is is called a twisted reduced model. The 


twist is (6 


—i\ 


nr 


The noncommutative spacetime 7 Zg is effectively hidden 


in the inhnitelymany degrees of freedom of the large N matrices C t and it 
reappears from expanding the matrix model around its classical vacuum (this 
is a dynamical emergent of spacetime). It is formally gotten by reduction of 
ordinary Yangs-Mills theory (with background hux) [53:] 

The (straight) open Wilson line has a particularly simple form in this 
matrix model formulation: 


0(C k ):= j d d xtr(U(x j , C k )) * e * fc - a: )Tre ifc ' e ' < ^ (2.118) 

which is manifestly gauge-invariant under Ci —>■ gCig _1 
This follows easily from: 

0{C k ) = Tr {U{C k )b{C^e ik - x ) (2.119) 
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2.5. THE CLASSICAL LAGRANGIAN AND HAMILTONIAN 

DYNAMICS OF MATRIX MODELS 


2.5 The classical Lagrangian and Hamiltonian 
dynamics of matrix models 

The fundamental idea is to set up an analog of classical dynamics in which 
the phase space variables are non-commutative, and the basic tool that allows 
one to accomplish this is cyclic invariance under a trace. Since no assump¬ 
tions about commutativity of the phase space variables (such as canonical 
commutators/anticommutators) are made at this stage, the dynamics that 
we set up is not the same as standard quantum mechanics. 

2.5.1 Bosonic and fermionic matrices 

We shall assume finite-dimensional matrices, although ultimately an ex¬ 
tension to the infinite-dimensional case may be needed. The matrix elements 
of these matrices will be constructed from ordinary complex numbers, and 
from complex anti-commuting Grassmann numbers. Just as a complex num¬ 
ber can be decomposed into real and imaginaty parts 

C = C R + iCj (2.120) 

a complex Grassmann number can be decomposed into real and imaginary 
parts 

X — X R + iXi (2.121) 

With Cri.j and X R j is real. Real Grassmann numbers are built up as products 
of a basis of real Grassmann element X 1? X 2 ,.... which obey the anticommu- 
tative algebra {X r ,X s } = 0 

Let B\ and B 2 be two N x N matrices with matrix elements that are 
even grade elements of a Grassmann algebra over the complex numbers, and 
let Tr be the ordinary matrix trace, which obeys the cyclic property 

Tr LmD 2 = D B l)mit(^2)nni =e<* 2 )nm(£i)mn = Tr B 2 B 1 (2.122) 

m,n m,n 

Similarly, let X\ and X 2 be two N x N matrices with matrix elements that 
are odd grade elements of a Grassmann algebra over the complex numbers, 

, which anticommute rather than commute, so that the cyclic property for 
these takes the form 

TtX\X 2 = ^^(T 1 ) mn (A’2) nm = ~'^^ / (X 2 )nm(Xl)mn = ^TrX 2 X 1 (2.123) 

m,n m,n 
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Final bilinear cyclic identity 


Tr BX = Tr XB. 


(2.124) 


We shall refer to the Grassmann even and Grassmann odd matrices B, X as 
being of bosonic and fermionic type, respectively. Clearly, operators that are 
of mixed bosonic and fermionic type can always be linearly decomposed into 
components that are purely bosonic or purely fermionic in character. 

The extra minus sign that appears in the odd grade case of Eq (2.123) has 
implications for the adjoint properties of matrices. Letting O 9 be a matrix 
of grade g, we define the adjoint by 


(0% n = {O' 


(2.125) 


Letting now Of 1 and O a 2 be two matrices of grade gi and g 2 respectively, 
this definition implies that 


(ofop)L = (of = Ew 

k 

= (-i) M2 £(0?)U0f );* = Y* 0 ? )L(of )L 

k k 

= (-ir»(Of)*Of) , ) m „ ( 2 . 126 ) 


The cyclic/anticyclic properties of Eqs ( 2.122|2.124 ) are the basic identities 
from which further cyclic properties can be derived. For example, from the 
basic bilinear identities one immediately derives the trilinear cyclic identities 


TvBi[B 2 , B 3 \ 
Tr Bi{B 2 , B 3 } 
Tr B{X 1 ,X 2 } 
TtX 1 {B,X 2 } 
Tr ,X[B u B 2 \ 
Tr ,X{B 1 ,B 2 } 
TvX 1 {X 2 , X 3 } 
TyX\ [X 2 , X 3 ] 


TvB 2 [B 3 ,Bi] = TtB 3 [B 1 ,B 2 \, 

TiB 2 {B 3 , B ]} = Tr B 3 {B 3 , B 2 }, 

Tr X 1 [X 2 ,B] = TrX 2 [Xi,B], 

Tr {X 3 , B}X 2 = Tr[X u X 2 ]B, (2.127) 

TrB 2 [,X,B 1 ]=TrB 1 [B 2 ,,X], 

TiB 2 {,X,B 1 } = TrB 1 {B 2 ,X}, 
TrX 2 {X 3 ,X 1 } = TrX 3 {X 1 ,X 2 }, 

TrX 2 [X 3 ,X 1 ] = TrX 3 [X 1 ,X 2 ], 


2.5.2 Derivative of a trace with respect to an operator 

The basic observation of trace dynamics given the trace of a polynomial 
P constructed from non-commuting matrix or operator variables, one can 
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define a derivative of the complex number TrP with respect to an operator 
variable O by varying and then cyclically permuting so that in each term the 
factor SO stands on the right. This gives the fundamental definition 

5TrP 

STyP = Tv——SO (2.128) 

60 

which for arbitrary infinitesimal SO defines the operator 5TrP /SO. In general 
we will take O to be either of bosonic or fermionic (but not of mixed), and 
we will construct (iTrP to always be an even grade element of the Grassmann 
algebra. (When P is fermionic, we can always make it bosonic by multiplying 
it by a c-number auxiliary Grassmann element a). With these restrictions, 
for SO of the same type as O, the operator derivative 5TrP /SO will be of the 
same type as O, that is, either both will be bosonic or both will be fermionic 

Suppose that P is a bosonic monomial containing only a single factor of 
the operator O , so that P has the form 

P = AOB , (2.129) 

with A and B operators that in general do not commute with each other 
or with O. THen when O is varied, the corresponding variation of P is 
SP = A(SO)B, and so cyclically permuting B to the left we have 

= e s Tr PASO, (2.130) 

= e B BA, (2.131) 

{ e B = 1 when the operator B is bosonic 
e B = — 1 when the operator B is frcmionic. 

The operator product AO is of the same bosonic or fermionic type as B, 
so we have e B = e A o an d could equally well write 

STtP 

— = e AO BA , (2.132) 

Suppose that P is a bosonic monomial containing two factors of the operator 
O that is being varied, and so has the general structure 

P = AOBOC, (2.133) 


STtP 

STyP 

SO 
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with A, B, and C operators that in general do not commute with each other 
or with O 

The variation of P is 


6P = A{60)B0C + A0B{60)C (2.134) 

Thus we have in this case 


5TtP 

6TtP 

60 


Tr (e AO BOCA(60) + e c CA0B(60 )), 
e A oBOCA + e c CAOB, (2.135) 


{ ec = 1(—1) According as whether C is bosonic (fermionic) 

6ao — 1(—1) According as whether the product AO is bosonic (fermionic). 

Let us expand 6TtP/60 is the form 


STvP 

60 


J2 C " R n 


(2.136) 


with the K n distinct Grassmann monomials that are all c-numbers (i.e., mul¬ 
tiples of the N xN unit matrix), and with the C n complex matrix coefficients 
that are unit elements in the Grassmann algebra 

Let us choose 60 to be an infinitesimal a times Cj, , with a a real number 
when O is bosonic, and with a a Grassmann element not appearing in K p 
when O is fermionic. There must be at least one such element, or else K p 
would make an identically vanishing contribution to Eq (2.128), and could 
not appear in the sum in Eq (|2. 136) We then have 


]T TrClC n K n a = 0 (2.137) 

o 


the coefficients of all distinct Grassmann monomials must vanish separately, 
we have 


TrCpCp = 0 


(2.138) 


We conclude that 


<5TrP 

~60 


= 0 


(2.139) 
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When O is bosonic, a useful extension of the above result states that the 
vanishing of STrP for all self-adjoint variations SO. The C n split into self- 
adjoint and anti-self adjoint parts, 


c n = cr 


+ c: 


(2.140) 


with C™ = and C% sa = C. For self-adjoint 50 implies that Tr C^ a 50 
is real, and Tr C/f a SO is imaginary. The vanishing STrP implies that both 
of these traces must vanish separately. Taking 50 = Cp a then implies the 
vanishing of (7®°, while taking 50 = iCp Sa then implies the vanishing of Cp Sa . 

In our applications, we shall often consider trace functionals STrP that are 
real, which will be true when the adjointness properties of the operators from 
which P is constructed imply that P — P^ is either zero or is an operator with 
identically vanishing trace. Real trace functionals STrP have the important 
property that when O is a self-adjoint bosonic operator, then STrP/50 is 
also self-adjoint. To prove this, we make a self-adjoint variation O , and use 
the reality of STrP to write 


0 


IrnTrdTrP oc Tr 


Ti 'SO 


STrP 

SO 


so - (soy 


STrP 

f STrP\ t_ 

SO 

V SO ) \ 


/ JTrP\ t_ 


(2.141) 


2.5.3 Lagrangian and Hamiltonian dynamics of matrix 
models 

Let L [{qy}, {qy}] be a Grassmann even polynomial function of the bosonic 
or fermionic operators {q r } and their time derivatives {q r }. The discrete 
index r labels the matrix degrees of freedom for a general matrix dynamics. 

Just as a classical dynamical system can have any number of degrees of 
freedom, the numbers tib and rip- of bosonic and fermionic operators {q r \ 
are arbitrary, and are unrelated to the dimension N of the matrices that 
represent these operators. From L, we form the trace Lagrangian 

L [{q r }, {<A}] = TrL [{q r }, {g r }] (2.142) 

The corresponding action 


S 



(2.143) 
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Using the definition of Eq (2. 128|) for get variations of the action 


0 = S = 


<5l . 

y * Tr X (^ A -+^<h 


(2.144) 


or after integrating by parts in the second term and discarding surface terms 

( SL d SL \ 


0 = S = 


/dSTr WA * ItW 

J ^ \5q r dt Sq r J 


(2.145) 


For this to hold for general same-type operator variations Sq r , the coeffi¬ 
cient of each Sq r in Eq (2.145) must vanish for all t, giving the operator 

(2.146) 


Euler—Lagrange equations 


hL d SL 

— -- — = 0. 

oq r dtoq r 


by the definition of Eq. (2.128), we have 

SL 


SL 

5q r 




S{^qr)i 


= 0 , 


(2.147) 


for each r the single EulerLagrange equation of Eq. (2.146) is equivalent to 
the N 2 Euler-Lagrange equations obtained by regarding L as a function of 
the N 2 matrix element variables (q r )ji ■ 

Let us now define the momentum operator p r conjugate to q r by 

P r =^ (2.148) 

d(jr 

p r is of the same bosonic or fermionic type as q r . We can now introduce a 
trace Hamiltonian H by analogy with the usual definition 


H = Tr ^2 p r q r — L 


(2.149) 


The variation is 


hH = Tr ^2 {{$Pr)<ir + p r Sq r ) - Tr ^ + 


5q r 


TT^2((Sp r )q r — p r Sq r ) 

r 

Tr ^ { e rQr)Spr ~ PrSq r ) • 


(2.150) 
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Therefore the trace Hamiltonian H is a trace functional of the operators {q r } 
and {p r } 


H = H [{q r },{p r }], 
with the operator derivatives 

AH AH 

— Pn ~ ^rQr 

dq r dp r 


(2.151) 


(2.152) 


e r — 1 According to whether q r ,p r are bosonic 

e r = — 1 According to whether q r ,p r are fremionic 


2.5.4 Trace dynamics models with global supersymme¬ 
try 

The Wess-Zumino model 

We begin with the trace dynamics transcription of the Wess-Zumino 
model. The trace Lagrangian for the WessZumino 


L = J d 3 x Trf - l -(d,Af - - xYd„X + + l -G 2 

- rniAF + BG - xx) ~ A[(A 2 - B 2 )F + G{A, B] -2 X(A- 


•y 1,2,3 are real symmetricand 7 0 , iy 5 are real skew-symmetric. A,B,F,G self- 
adjoint N x N lorentz matrices and x a Grassmann 4-component column 
vector spinor. The notation x is defined by y = y T 7 °. The numerical 
parameters A and m are respectively the coupling constant and mass. 

Taking operator variations of Eq (2.153) by using the recipe of Eq (2.128) 
, the Euler-Lagrange equations of Eq. (2.146) take the form 


d 2 A = mF + X({A,F} + {B,G}-2xx), 
d 2 B = mG + \(-{B,F} + {A,G}-2ix l5X ), 
l tl d/iX = m X + A({A, x} ~ {B, 75 %}, (2.154) 

F = mA + X(A 2 ,B 2 ), 

G = mB + A{A, B}. 
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Transforming to Hamiltonian form, the canonical momenta of Eq (2.148) 


P x = ~X1° = X T , 
Pa = d Q A, 

Pb = d 0 B , 


(2.155) 


and the trace Hamiltonian is given by 
H - 1 


= J rfVlW-[Pj + P| + (VA)‘ + (VB) 2 ] + P x 7 “fV X 
+ \(F 2 + a 2 ) - mxx - AB x7 °{A - *75B, M. 


(2.156) 


The trace three-momentum P. wgich together with H forms the trace four- 
momentum P CT , is given by 


P 


^xTt^PaVA + PbVR + P x Vy). 


(2.157) 


Let us now perform a supersymmetry variation of the fields given by 


5A = ey, 

5B = ie 75X, 

Sx = ^[F+ i-f 5 G+ x fl d fl {A + ix 5 B)}e, (2.158) 

5F = e'fd^x, 

5G = itrf*n*d v x, 

(2.159) 


with e a c-number Grassmann spinor (i.e., a four-component spinor, the spin 
components of which are lxl Grassmann matrices). Substituting Eq (2.158) 
into the trace Lagrangian of Eq (2.153). The variation of L is given by 


SL = J d 3 Tr(,P<9 M e) (2.160) 

r = ~xr[(Ydu + m)(A + i l5 B) + X{A 2 -B 2 + i l5 {A,B}pA61) 
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which identifies the trace supercharge Q a as 

Q a = I d 3 TrJ°a (2.162) 

= I d 3 Ti±(P x + x T )l(.rd» + m)(A + i l5 B) (2.163) 
+ \(A 2 — B 2 + i'y s{A, R})]cr. 


2.5.5 The supersymmetric YangMills model 


As a second example of a trace dynamics model with global supersym¬ 
metry, we discuss supersymmetric YangMills theory. 

We start from the trace Lagrangian 


L = / d 3 x Tr 


- xi J 'D„x + \d 2 


(2.164) 


with the field strength and covariant derivative D /t constructed from the 
gauge potential 7L according to 


F/xu — d^A u — d u An + [Ap, A u ], 

D^O = d„0 + [A„, O] (2.165) 

with 


D,F vX + D u F Xfi + D X F^ — 0. 
The Euler-Lagrange equations of motion are 


D = 0, 

i k F> v x = 0, (2.166) 

= 2g 2 xYx 


as usual for a gauge system, the u — 0 component of Eq. 
dynamical evolution equation, but rather the constraint 


(2.166) is not a 


F>iF 10 = 2g 2 xi°X. 


(2.167) 


Going over to the Hamiltonian formalism, the canonical momenta are given 
by 

PA, = -j 2 F ol , P x = X T , (2.168) 
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and the axial gauge trace Hamiltonian is 

H = H^ + H y , (2.169) 

with 

Hi = J d 3 x1t(^fpp%- gLf£ (2.170) 

— 2^2^ 1J ^ 2 ~ + [^1.7l2]) 2 — --^{d^Ai + (t?3^l2)“]") . 

1 

Fqs = -g 2 J dz'e(z- z')[-(p x x + x T Px) + D iP Al + D 2 pa 2 ]\ z ' (2.171) 


H x = / d 3 xTv(p x x°XiDix) 


(2.172) 


where we have taken care to write H in a form symmetric in the identical 
quantities p x and y T , and where e(z) = 1(—1) for z>0 (z<0). The trace 
three-momentum is 


/ d 3 xTr f F ri 
•* ' 1=1 


ip Al +P x Dix 


(2.173) 


and the conserved operator C of Eq. (2.6) is given by 

2 


C= f d 3 x(Y^[A llPAl }-x,P x 

J ^ z=i 


(2.174) 


with a contraction of the spinor indices in the final term of Eq. (2.174) un¬ 
derstood. By virtue of the constraint of Eq. (2.167), the conserved operator 
C can also be written as 


C = 


3 

f d3x dipA ‘ = - [ 

J l=1 J Sf 


d?Sip Al 


(2.175) 


sphere at 00 


which vanishes when the surface integral in Eq. (2.175) is zero. The corre¬ 
sponding conserved current C M , of which C is the charge, is given by 


(> = -[4,,FH + 2X7 M X- 


(2.176) 
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The conserved trace quantity N and the corresponding conserved current 
have the same form as in the WessZumino model 

N = — i J d 3 x Try T y (2.177) 

NT = iTrxYx 

Making now the supersymmetry variations 

$Ap = ige^^x (2.178) 

&X = (§“[%*» 7 *]^ 75jD ) 6 ’ 

SD = it'is'fDpX 

in the trace Lagrangian, with e again a c-number Grassmann spinor, we find 
using cyclic invariance under the trace and the 7 matrix identities D that 
whene is constant, the variation vanishes. When e is not a constant, the 
variation of L is given by 

*L = -ijteT r{m , (2,79) 

f = {- g xYF„ P [Y,Y ho } 
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2.6 Matrix model approach 

Field theories and especially gauge theories admit many classical solu¬ 
tions: solitons, instantons and branes, which play important roles in non- 
perturbative physics. Using this approach, it was possible to study many 
nonperturbative features of noncommutative held theory such as solitons 
and instantons (for more details see 


62 in references). 


We will call this approach matrix model approach, as the gauge theory 
can be described as a matrix model having the noncommutative space as its 
ground state, the fluctuations creating the gauge theory. Therefore we are 
looking for spaces that can be represented as fnite-dimensional matrix alge¬ 
bras, where everything is well defned. The space on which we will base our 
constructions will be the fuzzy sphere 63 , an N-dimensional matrix algebra 


corresponding to a truncation of the spherical harmonics on the sphere at 
angular momentum N — 1. In four dimensions, using the product of two 
such fuzzy spheres S'fj x Sj^, generated by W-dimensional matrices. When 
the two spheres they have a same limit, this fuzzy space goes over to the 
product of two commutative spheres, but in a diferent limit, it also goes to 
noncommutative R 4 with canonical commutation relations. 


2.6.1 The Heisenberg algebra 

In two dimensions, the coordinate algebra with canonical deformation 

[x,y]=i6 (2.180) 

The noncommutativity isn’t between the coordinates and momenta, it is 
between the coordinares themselves. We can use the usual Fock space rep¬ 
resentation for this aglbebra. 

The Fock space is given by 



T-L — { \n >, n G N 0 } 

(2.181) 

And defining 


x± x ± iy 

(2.182) 

with 


[x+, xJ\ = 29. 

(2.183) 
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Now define the creator and annihilation (or destruction) operators usually 
introduced in the process of second quantization 

x _ | n >= V29Vn + 1 | n + 1 >, x + | n >= V29Vn | n - 1 > (2.184) 

One can then consider a Fock space with a basis |n > (n ^ 0) provided by 
the eigenfunctions of the number operator N, 

N = x_x+ (2.185) 

N\n >= n\n > (2.186) 

and the vacuum state |0 > defined so that 

x_|0>=0 (2.187) 

For 6 ~ 0 one can write 

N = x_x + « ( x 2 + y 2 )/29 = r 2 /2d (2.188) 

so that conhguration space at infinity can be connected with n —> oo in Fock 

space. 

This can be generalized to higher dimensions. Any 2n-dimensional alge¬ 
bra with commutation relations. As we will mostly be concerned with the 
4-dimensional case in the following. 

The most general noncommutative 7 Zg is generated by coordinates subject 
to the commutation relations 

[x l x,x v \=i6 lw , fjt,u e {1, ...,4} (2.189) 

9^ can always be cast into the form 

0 0i2 0 0 \ 

-0i2 0 0 0 

0 0 0 0 34 

0 0 -034 0 

To simplify the following formulas, we restrict our discussion from now on to 
the selfdual case 



0 


flU 


1 

2 


per 


(2.190) 
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and denote 

9 := 9 12 = $34 

The complex coordinates 

x± L := Xi ±ix 2 , x±R := x 3 ± ix 4 , 


the commutation relations (2.189) take the form 


[x +a ,X- b ] = 29 Sab, [ X+a,X +b ] = [X- a , X- b ] = 0 
and the standard basis 

'H = {\n 1 ,n 2 >, ni,n 2 e N 0 } 


(2.191) 


(2.192) 


(2.193) 


(2.194) 


with 


X-L | 7 l 1; 7 l 2 > 

x +L I ni,n 2 > 
X-r | ni, n 2 > 
x +R | n 1 ,n 2 > 


\f29\Jn\ + 1 | ni + 1, n 2 > 
'J29^fn{ | ni - l,n 2 > 
v / 2$\/n 2 + 1 | ni,n 2 + 1 > 
\[29^Jn~2 | ni, n 2 - 1 > 


(2.195) 


2.6.2 Noncommutative gauge theory 


Let us start by considering the noncommutative version of a pure U( 1) 
gauge theory. We can introduce gauge theory by using a matrix action 




(2.196) 


where the AT are infinite-dimensional matrices, and the trace is over the Fock 


space (2.194) 


The action is invariant under unitary transformations 

-»• rtXpU 

The fluctuation A^ around the ground state x M as 

-Afi x^i T A^i 


(2.197) 


(2.198) 
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The fluctuations A M are understood as infinite-dimensional matrices acting 
on the Fock space (2.194) as well. They have to transform as 


U] + rfApU 


(2.199) 


to make the X M gauge covariant. The gauge covariant field strength then 
reads 


iF^ u = ([X^, X v ] - = [x^ A v \ - [x v , A + [A^ A v \ 


and the action (2.196) becomes 

s = ^tr(F^F^) 

We can also use the complex covariant coordinates 


X +L = X 1 +iX 2 , X +R = X 3 + iX 4 

X_ L = X l -iX 2 , X_ R = X 3 -tX 4 

and the corresponding field strength 

Feta,/3b [Acm, A^gj,] 2 06 a f/jS a l} 


( 2 . 200 ) 

( 2 . 201 ) 

( 2 . 202 ) 

(2.203) 

(2.204) 


with a,b G {L, R} and a, f3 E {+, —The action (2.196) can xow be written 
in the form 


S = “^2 tr £ F+ a ,- a F+ a - a - F +a,+b F -a,- 


(2.205) 


a,b 


and the equations of motion are given by 

iF aa ^] = 0 (2.206) 

a,a 


2.6.3 U( 1) instantons on TZj 


We will look for solutions of the equation of motion (2.206) which can be 
understood as instantons of the gauge theory 

In Ref 64 they schowing it’s a in noncommutative lZ 2 e that the noncom- 
mutative gauge theory contains the classical and quantum dynamics of all 
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U(N) gauge theories and that classical solutions are labeled by the rank of the 
gauge group and the magnetic charge. Also the BFS solutions describe vari- 
ouns D-l string attached. They can be interpreted as localized flux solutions, 
sometimes called fluxons. 

The situation on IZj is more complicated, and there are diferent types of 
non-trivial U (1) instanton solutions on Tig 

Assuming that 6^ is self-dual, there are two types of instantons: 


1 -There there exist straightforward generalizations of the two-dimensional 
localized fluxon solutions with self-dual held strength. As in the two-dimensional 
case, we will refer to these 4-dimensional solutions as fluxons. 
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2 -There are other types of U (1) instantons on 7 lg, in particular anti- 
sclfdual instantons which are much less localized than the fluxon solutions. 
For the construction of the fluxons, let us consider a finite dimensional 


subvectorspace V n of the Fock-space H of dimension n spanned by finite set 
of vectors |ni,n 2 >G TL 

V n =<{\i k ,jk>] k — 1 ,...., n} > . (2.207) 

We introduce a partial isometry S mapping TL to TL\V n , which has 

S*S = 1, (2.208) 

SS* = 1 -P Vn (2.209) 

where P n : TL —>■ V n is the orthogonal projection. By unitaty gauge transfor¬ 
mation we can assume that V n is spanned by the vectors |0 >,. |n — 1 > 

with the projection operator onto the subspace V n 

n 

Pv n • ^ ^ \^kijk iki jk\- (2.210) 

k= 1 

According 64 the solutions of the equation of motion given by 

n 

X\ L := Sx +L S ] + I i k ,jk >< ik,3k\ (2.211) 

k =1 
n 

Xl R ■■= Sx +R s t + J2^\ik,jk >< tk,Jk\ (2.212) 

k =i 


51 






2.6. MATRIX MODEL APPROACH 


With 


\Kl, Kr) = \Kl, X-r\ = \*-l, Kr\ = X-r] = 0 (2.213) 

and X” a = (X” a )t, ^’ R <g C determine the position of the fluxons. The field 
strength fo this solution is 

iv = Ay V (2.214) 

The action corresponding to the instanton solution is proportional to the 
dimension of the subspace V n 

87T 2 87T 2 

= ~^n. (2.215) 

’’pseudoparticle” solutions we mean the long range fields A^ which minimize 
locally the Yang-Mills action S and for which S(A) < oo. The space is 
euclidean and four-dimensional. All fields we are interested in satisfy the 
condition: 


Fpv — d^A v — d^A^ + [A^, A u \ 


(2.216) 


This field —> 0 when x —> oo. 

Let us start consider S 3 is very large sphere in 4-dimensional space. From 

flT 216 t 



s 3 



9g(x) 

dx„ 


s 3 


(2.217) 


Where g(x) are matrices of the gauge group. Hence every field A^x) produce 
a certain mapping of the sphere S 3 onto the gauge group G. The phase space 
of the Yang-Mills fields are divided into an infinite number of components, 
each of which is characterized by some value of q, where q is a certain integer. 

We search for the absolute minimum of the given component of the phase 
space. In order to do this we need the formula expressing the integer q 



s 3 



dg(x) 

dXfj, 


s 3 


(2.218) 
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Chapter 3 

Yang-Mills Matrix Theory 


Modern particle theories, such as the Standard model, are quantum Yang- 
Mills theories. In a quantum held theory, space-time fields with relativistic 
held equations are quantized and, in many calculations, the quanta of the 
holds are interpreted as particles. In a Yang-Mills theory these helds have 
an internal symmetry: they are acted on by a space-time dependant non- 
Abelian group transformations in a way that leaves physical quantities, such 
as the action, invariant. These transformations are known as local gauge 
transformations and Yang-Mills theories are also known as non-Abelian gauge 
theories. 

Yang-Mills theories, and especially quantum Yang-Mills theories, have 
many subtle and surprising properties and are still not fully understood, 
either in terms of their mathematically foundation or in terms of their phys¬ 
ical predictions. However, the importance of Yang-Mills theory is clear, the 
Standard Model has produced calculations of amazing accuracy in particle 
physics and, in mathematics, ideas arising from Yang-Mills theory and from 
quantum held theory, are increasingly important in geometry, algebra and 
analysis. 

Consider a complex doublet scalar held (j) a ; a scalar held is one that has 
no Lorentz index, but, as a doublet, <p a transforms under a representation of 
SU(2), the group represented by special unitary 2x2 matrices: 

4>a{x) -»■ gabfox (3.1) 

where g E SU( 2) and the repeated index is summed over. If this is a global 
transformation, that is, if g is independent of x, then derivative of «p a have 


53 



(3.2) 


the same transformation property as a itself: 

d(f) a x dg ab (j) b _ d(j) b dg ab 

r-t ' r-t {Jab 

ux ^ ox ^ ox ^ ox ^ 

In order to construct an action which includes derivatives and which is invari¬ 
ant under local transformations, a new derivative is defined which transforms 
the same way as (f) a : 

D lt 0 a = ^ + (A,) ab( j> b (3.3) 

where is a new two-indexed space-time field, called a gauge field or gauge 
potential, defined to have the transformation property 

do 

(AOm —>■ g a b{A^) cd g d ^ - -^r9cb (3-4) 

Now, under a local transformation 


D t g a bD 


(3.5) 


and so, D^a transforms in the same way as (j) a . This derivative is called 
a covariant derivative. A physical theory which includes the gauge field A 
should treat A M as a dynamical held and so the action should have a kinetic 
term for A^ . In other words, the action should include derivative terms for 
A m . These terms are found in the held strength 

f)A 

F,w= d&~ dA ^ XV + [A ^ A ^ (3 ' 6) 

which has the covariant transformation property 

{F^a/) a b t gab{F^i/) c d9db (3-7) 

where [A M , A v \ is the normal matrix commutator. In fact, the simplest Yang- 
Mills theory is pure Yang-Mills theory with action 

S[A] = ~ f d 4 x trace F IJU F IW . (3.8) 


and corresponding held equation 


dF 

F' 1 fus 

dx M 


0 


(3.9) 
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Solutions to this equation are known as instantons. 

More generally, Yang-Mills theories contain gauge fields and matter fields 
like (j) and fields with both group and Lorentz or spinor indices. Also, the 
group action described here can be generalized to other groups and to other 
representations. In the case of the Standard Model of particle physics, the 
gauge group is SU (3) x SU (2) x U(l) and the group representation structure 
is quite intricate. 

Yang-Mills theory was first discovered in the 1950s, at this time, Quan¬ 
tum Electrodynamics was known to describe electromagnetism. Quantum 
Electrodynamics is a local gauge theory, but with an abelian gauge group. 
It was also known that there is a approximate global non-Abelian symmetry 
called isospin symmetry which acts on the proton and neutron fields as a 
doublet and on the pion fields as a triplet. This suggested that a local ver¬ 
sion of the isospin symmetry might give a quantum held theory for the strong 
force with the pions fields as gauge fields [ORaifeartaigh, 1997]. This did not 
work because pion holds are massive whereas gauge helds are massless and 
the main thrust of theoretical effort in the 1950s and 1960s was directed at 
other models of particle physics. 

However, it is now known that the proton, neutron and pion are not 
fundamental particles, but are composed of quarks and that there is, in 
fact, a quantum Yang-Mills theory of the strong force with quark helds and 
gauge particles called gluons. Furthermore, it is now known that it is possi¬ 
ble to introduce a particle, called a Higgs boson, to break the non-Abelian 
gauge symmetry in the physics of a symmetric action and give mass terms 
for gauge holds. This mechanism is part of the Weinberg-Salam model, a 
quantum Yang-Mills theory of the electroweak force which is a component 
of the Standard Model and which includes both massive and massless gauge 
particles. 

These theories were only discovered after several key experimental and 
theoretical breakthroughs in the late 1960s and early 1970s. After it became 
clear from collider experiments that proton have a substructure, the oretical 
study of the distance dependant properties of quantum Yang-Mills theory 
lead to the discovery that Yang-Mills helds are asymptotically free [Gross, 
1999]. This means that the high-energy behaviour of Yang-Mills helds in¬ 
cludes the particle like properties seen in experiments, but the low-energy 
behaviour may be quite different and, in fact, the quantum behaviour might 
not be easily deduced from the classical action. Confinement and the mass 
gap are examples of this. The strong force is a local gauge theory with quark 
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fields. The quark structure of particles is observed in collider experiments, 
but free quarks are never detected, instead, at low-energies, they appear to 
bind together to form composite particles, such as neutrons, protons and 
pions. This is called confinement. It is possible to observe this behaviour in 
simulations of the quantum gauge theory of the strong force, but it has not 
been possible to prove mathematically that confinement is a consequence of 
the theory. The same is true of the mass gap, it is known that particles have 
non-zero mass, and this is observed in simulations, but, there is no known 
way of deriving the mass gap mathematically from the original theory [Clay, 
2002 ], 

The symmetries of Yang-Mills theory can be extended to include a global 
symmetry between the bosonic and fermionic fields called supersymmetry. 
While there is no direct evidence for supersymmetry in physics, the indirect 
case is very persuasive and it is commonly believed that direct evidence will 
be found in the future. Often, supersymmetric theories are more tractable, 
for example, Seiberg and Witten have found exact formula for many quantum 
properties in N — 2 super-Yang-Mills theory [Seiberg , Witten 1994], It is 
also commonly believed by theoretical physicists that the quantum Yang- 
Mills theories in particle physics are in fact a limit of a more fundamental 
string theory. 

3.1 Yang-Mills Theory in Low Dimensional 

3.1.1 Bosonic Matrix Integrals 

We begin our systematic study of Yang-Mills theories on low-dimensional 
with the simplest possible example, namely 0-dimensional Yang-Mills coupled 
to p adjoint scalars with a quartic interaction. There are no Yang-Mills fields 
in 0-dimensions and the effect of gauging is trivial, introducing only an overall 
factor of the gauge group volume. As a result, the partition function that we 
wish to study takes the form of a simple matrix integral 

Z= j D ^expj - j ( 3 . 10 ) 
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We rescaled the make clear that the effective couplings are Ao /m?rrij, the 
partation function becomes 


Y = 



-Y 


A#.,*, 


rrrrn 


3 


(3.11) 


3.1.2 Bosonic Yang-Mills on S' 1 

We now move np in dimension to consider bosonic Yang-Mills on an S 1 . 
The particular model that we consider is 

N r f v p p \ 

S = 5T / tr ( Y + Y M 2 i1 - Yl*- ) ( 3 - 12 ) 

'' A 1=1 i i<j ' 

From A ,M, and the size R of the S 1 , which we will think of as an inverse 
temperature, we can construct two dimensionless parameters 

t = (i?A 1/3 ) -1 m = M\~ 1/3 (3.13) 


Large masses m » 1 

The second regime that one can study corresponds to the limit of large 
mass m » 1 and displays a much more interesting and nontrivial structure. 
Here, the theory is effectively weakly coupled and, since all mode except the 
zero component of the gauge held, Ao , are massive, they can be integrated 
out in perturbation theory. 

We consider a function of a unitary matrix U 


Z = / DUeM~N 2 S eS{u) } U = e 


<f> dtAo 


(3.14) 


The computation of S e ff(U) now proceeds exactly as in the Af = 4 theory 
with the identification 


x = e ~ MR = e~ m/t , (3.15) 

the letter partition function is rather trivial as there is only one bosonic state 
at energy 1 for each scalar held 


z(x) = x 


(3.16) 
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This leads to a simple effective action for U 

OO 

S eS (U) = p V — ti{U n )ti(U~ n ) (3.17) 

n 

71=1 

where the factor of p arises because S e ff receives contributions from p scalar 
fields. 

Writing this in terms of moments p n of the corresponding eigenvalue dis¬ 
tribution, this becomes 


Seff(Pn) ~ ^ ^ m n\p 

n 


m r 


= 1 


px 


where the 1 in (3.19) arises from the Vandermonde measure as usual. 

t c = m(ln p)- 1 


(3.18) 

(3.19) 

(3.20) 


Away from the strict m —> oo limit, (3.18) will receive corrections that can 
potentially alter the order of the transition. In particular, as we have seen 
many times it is the sign of the quartic term in the effective potential for p± 
that is the difference between first and second order. This can be determined 
by computing higher-loop perturbative corrections to (3.18) and integrating 
out the modes pi with i > 1 that remain massive at the transition point. The 
result is an effective potential for p\ 


S(p 1 ) = ml(x,m \pi\ - ~b(x,m 3 )|pi| 4 + ... 

irr 


(3.21) 


The calculation of mj and b to the requisite orders is straightforward but 


tedious, yielding 66 


m\ = (1 — px) —-—-( p 2 — p)(x 2 + 2x)\nx 


b = 


Am 6 
A ip — 1) lnp 


(3.22) 


(ln(p)(9p 2 + 2 p) + Ap 6 + 7 p 2 — Ap — A)m 6 + . 


32 


p6 


Since b < 0 for all p > 1, we see that the phase transition is indeed of 
first order as one moves away from m — oo. From the correction to rn\ in 
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temperature, obtaining 66 


(3.23), we can also compute the leading order shift in the phase transition 

+ ... (3.23) 


t c = m(\np) 1 + 


1 (p — l)(2p + 1) 


4m 2 


plnp 


This result will fit in nicely the picture of the m = 0 theory that we now 
describe. 


3.1.3 Bosonic Yang-Mills on T 2 

We now consider the first of our two-dimensional models, namely bosonic 
Yang-Mills on a rectangular two-torus with radii R \, R 2 
The specific model that we consider has action 

S= J\,f d2 xtT ( F ™ + E + M ^ 2 i] ~ E>'’ $J ] 2 ) ( 3 ‘ 24 ) 

and can be parametrized by the two dimensionless radii 

r\ = RiVX r 2 = R 2 V A (3.25) 

and the dimensionless mass 


M 



(3.26) 


Infinite masses m = oo 


We begin our study of (3.24) by considering the strict limit m = oo. The 
reason we do this is to recall some elementary facts about the exact solution 
of pure Yang-Mills on a two-torus that will be useful when we study large 
but noninfinite masses. 

The exact result for the two-dimensional Yang-Mills partition function 
was first obtained by Migdal 67 using a lattice regularization of the theory. 
We consider putting unitary matrices Ul on the links of the lattice and 
consider the partition function 


Z 


n du l n Zp(u p ) 

L P 


(3.27) 
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where the product is over plaquettes P of the lattice and Zp is a plaquette 
action chosen so that the continuum theory coincides with Yang-Mills. One 
example of a plaquette action is the usual Wilson one 


Z P = exp|— tr(Up + Up 1 ) 


(3.28) 


One can imagine integrating out some of the links to generate an effective 
action for the larger plaquettes that remain. Under this RG flow, the action 


(3.28) is not invariant, instead flowing to the plaquette action 


Zp(U P ) = 


J2d RxR (U P )ex p j - ^C 2 (R) 

R ^ 


(3.29) 


where the representation R has dimension d R and quadratic Casimir C 2 (R) 
and A is the area of the new plaquettes in units of the fundamental ones. It 


is easy to verify that (3.29) satisfies the additivity property 


DUZ P (ViU, Ai)Z p f (rfv 2 , A 2 ) = Zppp! (V 1 V 2 , A 1 + A 2 ) (3.30) 


and consequently is an RG fixed point. To compute the partition function of 


2-dimensional Yang-Mills, we are equally justified in using (3.28) or (3.29) 


The former is easier to work with as the property (3.30) can be used to 
integrate out all of the links except those wrapping nontrivial cycles. Re¬ 
stricting to the torus for simplicity, we arrive at the following expression for 
the partition function 


Z= fDUDV E d R exp( 

■' R 


r 1 r 2 
' 27V 


C^XRiUVU-^- 1 ) 


(3.31) 


where we have used the definitions (3.25) of r l5 r 2 . The links which remain 


are precisely the holonomies that we wish to study. The remaining integrals 


(3.31) are now sufficiently simple that we can perform them exactly to obtain 


Z = l^e 

R 


-WC2(R) 


(3.32) 


This partition function is very well-behaved and exhibits no nontrivial phase 
structure. As an exercise, let us obtain this conclusion in a manner more 
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in line with the sort of analysis that we have used to study other theories 
above. In particular, let us consider integrating out only one holonomy, say 


V , from (3.31). We obtain 


Z = 


I duY,^ WC 2 {r) xr(u)xr(u t ) 

J R 


(3.33) 


Using the generalized Frobenius relations of 


write the integrand in terms of an effective action for U |67 

0 -r ir2 n + 2e- rir2n/2 )ti{U n )ti{U 


z = j 


69 , 1701 it is possible to 

(3.34) 


Introducing an eigenvalue density p as usual, we can obtain an effective action 
for the moments of p 

S'ffiftn) = £ 1(! - e- r,ri ” /2 ) 2 |p»| 2 (3.35) 

n 

We see that for all values of the coupling A A = r\r 2 the effective masses of 
the p n are positive. In the limit t\T 2 —> oo they become arbitrarily light but 
the never become tachyonic. The eigenvalue distributions are thus always 
uniform and there is no phase transition. 
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3.2 Yang-Mills Theories 

3.2.1 Yang-Mills Gauge Theories 

The Yang-Mills matrix models are related to gauge theories by dimen¬ 
sional reduction. We recall the structure of Yang-Mills gauge theories 
A Yang-Mills gauge theory in D dimensions has Lagrangian 

£ = \f 2 - (3.36) 

The fields in this theory are the vector potential X n p , and fermions ipa . Here 
H is a spacetime index running from 0 up to HI, and for the moment we use 
a Minkowski metric. The 7 ^ are Dirac matrices, and ijj is defined 

$ = ^ 7 ° (3.37) 

The fields are in the Lie algebra of a compact semi-simple gauge group G so 
we can write 


X, = x;t a a = rs (3.38) 

where the t a (a = 1,.., g) are the generators of the Lie algebra which we choose 
such that 

Tr t a t b = 2 5 ab (3.39) 

and 

[t a ,t b ]=if abc t c . (3.40) 

The gauge field strength F is defined 

f;„ = d„x: - g„x; + cf^x^xi ( 3 . 41 ) 

and the gauge covariant derivative is 

(V^) a = d,r + cf abc Xlr. (3.42) 

The parameter c is a coupling constant. The theory is invariant under gauge 
transformation 

ip -)• UipU~\ Xp ->• UX^U- 1 - tg-^dU)^ 1 (3.43) 

where U E G. The fermions are optional in this model. We can define a 
purely bosonic gauge theory by simply omitting them. 
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3.3 Yang-Mills Matrix Models 


To obtain a Yang-Mills matrix model, we take the Lagrangian (3.36) and 
assume all the fields are independent of space and time. Effectively, this 
means we drop all the derivative terms from (3.36) 

At this stage, we also move from Minkowski to Euclidean signature. We 
do this by setting 


X% = iX a D , (a — 1,.., g) (3.44) 

and taking the Xfj real . We also set 

7 ° = i-i D (3.45) 

We have been careful to leave this manipulation until last because we wish 
to study the Wick rotation of a Minkowski theory . This leads to a rather 
strange effect in the case of D = 10 when the fermions are Majorana. Since 
the Dirac matrices can no longer all be imaginary, an SO{D) transformation 
would break the Majorana condition. However, after integrating out the 
fermions, full SO(D) invariance is restored since we can analytically continue 
in X D . 

We arrive at the matrix model action 

Sym = -Tr^[Y M , X U \[X„ X v \ + (3.46) 

where we have dropped the coupling constant g since it can be scaled out in 
a trivial manner. In those cases where the fermions were originally Majorana 
(before Wick rotation), we may choose the representation in which the 0“ 
eal. In those cases where the fermions are complex, it will sometimes be 
convenient to rewrite the (0“) for each a as a real vector of double the 
length. We can also absorb the 7 0 which appears in the definition of 0 into 
the 7 m . Thus we shall sometimes write the action in the form 

Sym = -TrD[X„,X„][X p ,X„] + Iv«ryx,,,V7)]) (3.47) 

where the T^ are some new matrices defined in terms of the 7 ^ , and the 0 “ 
are now always real. In the case where the original fermions were complex, 
the range of the indices a and (3 has been doubled. 
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We define a partition function 
d M 


Zd,g = [ f[ dX, fj #«ex p ( J ^ Tr[Y M , AR] 2 + l -T^ a [T^X^ ^-48) 

/X— 1 Q!— 1 \ /XIX ' 

which we shall also sometimes refer to as the Yang-Mills integral. In principle 
one can integrate out the fermions to obtain 


z d,g = I n dXJ> D , G { X>xp ( ^ Tr[A /t , X, 

fl= 1 \ /XIX 


(3.49) 


where the Pfaffian Vd,g is a homogeneous polynomial of degree \Mg. In this 
representation, the gauge symmetry is 


X,, ->• tfXJJ, U eG 


and SO(D) symmetry 


X/j —> 




Y 


Q e SO(D) 


(3.50) 


(3.51) 


In addition, we shall consider simple correlation functions 


< c k (X a ) >= 


[ J] dXYk(X a )V D , G ( A^exp ( l - Tr[Y /t , A, 

/x= 1 ' /xix 


(3.52) 


with Ck a function of the X M which grows like a polynomial of degree k. 
The first question one must ask about these models is whether the in¬ 


tegrals (3.48) and (3.52) which define the partition function and correlation 


functions are well defined. Certainly, we must require at least that the par¬ 
tition function is finite for the theory to make any sense. The difficulty here 
is that the potential Tr [Y /t , X v \ [X fl , X„] has flat directions in which the ma¬ 
trices commute. For example, in the bosonic case, one can move to infinity 
along one of these directions whilst keeping the integrand constant, and thus 
it was widely believed that these integrals may be infinite. However, in the 
case of SU(2) it is possible to perform the integrals for the partition function 
exactly. This was done originally in the supersymmetric cases 71--74j and it 
was found that the partition function does converge at least for D = 4, 6,10. 
Subsequently, eigenvalue densities and some correlation functions have been 
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calculated in 75 . It was believed that the supersymmetric versions should 


be more convergent than the bosonic because the contributions from the 
fermionic integrals would be close to zero near the flat directions. However, 
the 577(2) bosonic partition function was calculated in 76 , and was found 
to converge when D > 5. 

The authors of 1761 were able to use Monte Carlo methods to calculate 


the supersymmetric integrals numerically for 57/(2) and 57/(3), and the cal¬ 
culations have been extended to various other gauge groups, and also to the 
bosonic theories 77,78 . A difficulty with numerical simulations for the su¬ 


persymmetric integrals is in performing the fermionic integrations to obtain 
the Pfaffian, and for this reason, the exact model has only been studied for 
the smaller gauge groups. However, the bosonic models have now been stud¬ 


ied for SU(N ) with N up to 768 78,79]. Analytic approximation schemes 
have also been constructed for the bosonic models in 80 and recently for 
the D — 4 supersymmetric model 


81 


The conclusions of the numerical methods are that the supersymmetric 
partition function converges when D — 4,6,10 and that the bosonic partition 


functions converge at least when D is large enough 78 
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3.4 Introduction of IKKT Models 


Ishibashi, Kawai, Kitazawa, and Tsuchiya (IKKT) have proposed a model 
(0 + 0)-dimensional matrix model should give a Paincare invariant describing 
of type-IIB string theory on a flat background. Beginning with the Schild 
gauge-fixed form of the string action, they found that the matrix regulariza¬ 
tion of the action led to a zero dimensional matrix model equivalent to the 
complete dimensional reduction of ten dimensional super Yang Mills theory. 

What is a matrix theory? Simply put, it is a quantum mechanics with 
matrix degrees of freedom. It is in general comprised of some N x N bosonic 
and fermionic matrices. Matrix theories have the attractive quality that they 
are not quantum held theories, and thus have none of the peculiarities in¬ 
volved with QFTs, such as renormalization. Indeed, since N is finite (though 
taken to be large) there are only a finite number of degrees of freedom. 

Witten found that the low-energy Lagrangian which describes a system of 
N type-IIA DO-branes is equivalent to the dimensional reduction (the theory 
on a point) of ten dimensional super Yang-Mills theory to (0+1) dimensions 



The dimensionally reduced Hamiltonian is 


H = ^Tr jpLP - ^ [X*, X j ] [X\ X j ] + 9 7i [X\ 6] 


(3.53) 


X 1 and P l are N x N with bosonic entries, 9 is a set of N x N matrices 
with fermionic entries. 

Banks, Fischler, Shenker, and Susskind (BFSS) discovery after Witten 
that this Hamiltonian precisely describes M-theory in the light-front coor¬ 
dinate system, in N — oo limit of supersymmetric matrix quantum me¬ 
chanics describing DO branes 83 . The IKKT model is one proposed non- 
perturbative definition of string theory 84 . It is sometimes referred to as 


the IIB matrix model, since it is related to a gauge-fixed form of the Green- 
Schwarz action for the IIB string. 

There are two ways of defining the model. One is as the complete dimen¬ 
sional reduction to 0-dimensions matrix model arising from the dimensional 
reduction in all ten dimebsions of M = 1 super Yang Mills theory. 


s = ~ 

9 2 


Tl '{7 £ tVVP",v ] + W</.rp„pj. 

^ ^ fi,v=0 Z m=0 ' 


(3.54) 
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The action (3.54) looks very similar to the dimensionally reduced Hamilto¬ 
nian of the BFSS conjecture, Eq (3.53).This action still has SO(9.1) symme¬ 
try, as well as the SU(N) gauge symmetry. The A M are N x N hermitian 
matrix components of a vector in an SO(9.1) representation. The fermion 
ip is in a Majorana-Weyl representation and N x N hermitian matrices for 
components 


85 


The partition function for the theory is then defined via euclideanization 
of the action (Wick rotation of A 0 and T 0 ); 


Z 


dAdipe Se 


(3.55) 


The model may also be defined in terms of a grand canonical partition func¬ 
tion or canonical ensemble: 


OO « 

Z\P\ = J2 / dXdVe-^l 

N =1 ■' 


(3.56) 


Here f3 is interpreted as a chemical potential dual to the matrix size N , and 
the action is 


c 2 —_ 

E 2a 2 P 


Tr 


\\A„AA[XrX-] + \^\A^] 


+ pN. (3.57) 


If the large N limit is smooth, we ewpect that the P —>■ p c limit is indentical 


to consider the microcanonical ensemble with fixed N and take N large 85 


From the actions (3.54) and (3.57), it is clear that the action is mini¬ 


mized if the bosonic matrices commute. If that is the case, they may all be 
simultaneously diagonalized and decomposed into a diagonal part X fl and an 
off-diagonal part ; 


A^ — X n + A n 



x. 


\ 


+ A-n 


\ 



The fermion pj is also decomposed into diagonal and off-diagonal com- 
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posents £ and 0 


0 = 


(e 


\ 


\ 


e) 




where and £^ satisfy the constraints = 0 and YuiLiCa = O’ 

respectively, since we may fix the U(l) part by translation invariance. Under 
this transformation, the off-diagonal terms A^ and 0 become massive. They 
may then be integrated out, and an effective action for the diagonal bosons 
(to be interpreted as space-time points) may be obtained by integrating out 
the diagonal fermions £* 85 . 




dAd^e s[A ^ ] = j dXdte~ SBs[x ’® 

= [ dXe- s ° s[x] 


(3.58) 

(3.59) 


Ishibashi et ah interpret the diagonal matrix elements as points of spacetime. 
By expanding the action as a perturbative series and integrating out all the 
other fields, an effective action for the spacetime points is found 

We perform integrations over off-diagonal parts A M and 0 by the pertur¬ 
bation coupling constant g 1 . This scheme is valid if the spacetime points are 
widely separated (\x l — ad |) yfg 85 . Residual gauge symmetries may be 


handled by introducing Faddeev-Popov ghosts and gauge-fixing terms 


Sg.f. + Sf.p. — 


2 g 


1 Tr([X M ,^] 2 )-4Tr([X M ,6][^,c]), (3.60) 


r 


where b and c are the Faddeev-Popov ghost held. 
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The original action (3.54) can be expanded as follows 86 


5 = S 2 + 


(3.61) 

S 2 = 4t r( 
r 

- [x„, i„] [A'", An + [a„, i„] [A", An 



- R"[A„y] - [f.ijr<y -|r"[4,e]), 

(3.62) 

-S'int = 4 tI-( 

r 

- 2{x„An{AnAn - l[i„ a.p'Ui 




(3.63) 


Can be written S- 2 + S g j in terms of the components as 

+ S,.,. = T £((4 - xifAfA^ - ^‘r"(4 - xl)r) 
y i<j 

+ (f - ?)?•>$* Af - e)A'>). (3.64) 


The first and the second terms are the kinetic terms for A and ^ respectively, 
while the last two terms are vertices. The basic building blocks of the 
Feynman rules (see [86]) 
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3.5 The IKKT Model of IIB superstring 

The supersymmetric Yang-Mills matrix theory with D — 10 has been 
proposed as a constructive definition of IIB superstring theory |24], We give 


a very brief introduction here, but for a review see 87 


The idea of the IKKT conjecture is to begin with the Green-Schwarz 
action for the superstring in the Scliild gauge: 


Sgs = / d 2 a 


V$ a ( x u } 2 — " 0 }) + P\/l 


(3.65) 


Here a are 2-dimensional world-sheet coordinates, g = det(g ab ) is the deter¬ 
minant of the world sheet metric, and a, /3 are parameters (which could be 
scaled out). The x^ are target space coordinates, and the Poisson bracket is 
defined 


{x,y} = 4= e ab xd a xd b y 
V9 


(3.66) 


The theory is then regularised essentially following a method of Goldstone 
and Hoppe (for a review, see [881). 

1. A function y on the world sheet is replaced by an N x N traceless 
hermitian matrix Y , with a correspondence 


J d 2 a\/Jy O TrY 


(3.67) 


2. and 


{x,y} -i[X,Y] 


(3.68) 


Performing this regularisation, the action (3.65) becomes 


Sikkt = -a (^Tr^, X v \ [X^ X v \ + ^Tr^r^Y^, ) + /37V (3.69) 

The string partition function is given by 


D[x]D[jp]ex.p(-S GS ) 


(3.70) 
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then becomes the matrix integral 


D 16 


n dx )< n #« ex p (-s IKKT ) 


(3.71) 


fi=l a=l 


which after scaling out becomes the Yang-Mills matrix partition function 
Z10, SU(N ), with an additional factor e& N . 

In their original proposal, IKKT interpreted the integral over the world 
sheet metric f D[g\ as a requirement to sum over N: 



(3.72) 


However, in general, the matrix regularisation procedure outlined above is 
valid in the limit N —> oo, and the partition function is often taken as the 
large N limit 


(3.73) 


ZlKKT ~ 2 W} SU(N) 


which would correspond to a more literal application of the Goldstone Hoppe 
regularisation. The large N limit is not yet well understood, and it is not 
clear exactly how to interpret the model. Nevertheless, an argument relating 
Wilson loops in the matrix model to string field theory in light-cone gauge 
provides additional evidence for the importance of the IKKT model 89 . 
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Chapter 4 
Fuzzy spaces 


Studies of fuzzy spaces cross over a variety of concepts in mathematics 
and physics. The basic idea of fuzzy spaces is to describe compact spaces in 
terms of finite dimensional (N x iV)-matrices such that they give a concrete 
realization of noncommutative (NC) spaces 90,91 . Use of fuzzy spaces in 


92 . Since then, fuzzy spaces 


physics was suggested by Madore around 1992 
have been an active area of research. 

Fuzzy spaces are described by finite dimensional matrices, due to the 
Cayley-Hamilton theorem, there is a natural cut-off on the number of modes 
for matrix functions on fuzzy spaces. So one can use fuzzy spaces to construct 
regularized held theories in much the same way that lattice gauge theories are 
built. Various interesting features of held theories on fuzzy spaces have been 
reported; for example, existence of topological solutions such as monopoles 
and instantons, appearance of the so-called UV-IR mixing, and evasion of 
the fermion doubling problem which appears in the lattice regularization 


4.1 Construction of Fuzzy Spaces 

4.1.1 Construction of fuzzy CP k 

Hilbert space 

A finite dimensional Hilbert space 1-Ln for fuzzy CP k = SU[k + 1)/U[k) 
(k=l,2,...) is given by holomorphic sections of a complex line bundle over 
CP k , the holomorphic sections of the complex line bundle should correspond 
to a unitary irreducible representation G = SU(k + 1). Notion of liolo- 
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morphicity in the representation of G can be realized by totally symmetric 
part of the representation, i.e., (n, 0), where n is the rank of the repre¬ 
sentation (n = 1,2,...). The other totally symmetric representation (0, n) 
corresponds to antiholomorphic part of the SU(k+l) representation and the 
(p, p)-representation gives real representation. 

For SU(2) (corresponding to k=l), the representation is given by a single 
component, say (p), so there is no real representation. The dimension of "Hat 
is then determined by that of the (n, (^-representation for SU(k+l); 

N k = dim(n , 0) = (4.1) 

k\ ni 

Consequently, matrix algebra of fuzzy CP k is realized by x - 
matrices. Operators or matrix functions on fuzzy CP k are expressed by 
linear combinations of N ^ x -matrix representations of the algebra of 
SU(k+l) in the (n,(^-representation. 

We begin with write down a holomorphic U(l) bundle as 

ly'lj) = VWV^l(g), (4.2) 

= <(n,0),m|j|(n,0),JV(*>> (4.3) 

where \(n, 0),m > (m = 1,2 ,, N^) denote the states on the Hilbert space 
"Hat, \(n,0), > is the highest or lowest weight state, g is an element 

of G — SU (. k + 1) and g is a corresponding operator acting on these state. 
'D^l ( g ) is known as Wigner ©-functions for SU (/c+1) in the (n, (^-representation, 
allowing us to interoret the ©-functions as matrix elements. 

Let Ra denote the right-translation operator on g; 

RAg = gtA (4.4) 

where Ia are the generator of G in the fundamental representation (1,0). The 
element g is given by g — exp(itA0 A ) with continuous parameters 9 A . We now 
consider the splitting of i^’s to those of U(k) = SU(k ) x U (1) subalgebra and 
the rest of them, i.e., those relevant to CP k . Let t 3 (j = 1 , 2 ,, k 2 ) and t^+ 2 k 
denote the generators of U(k) G SU(k + 1), t^+ 2 k being a U( 1) element of 
the U(k), and let t± t (i=l,2,,k) denote the rest of LRs. One can consider t±i 
as a combination of rasing-type (t +i ) and lowering-type (L) operators acting 
on the states of R N . Choosing |(n, 0), > to be the lowest weight state, 
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we then find 


= o (j — 1 , 2 , 3 ...,k 2 ), 


kV ( 2'°l(g) = 


nk 


\J 2 k(k + 1) 


D 


O,o) 

mN k 


(9) 


(9) = 0 


(4.5) 

(4.6) 

(4.7) 


(5.1) and (5.2) indicate that $1? ~ V^l(g) is a 17(1) bundle over CP k . 

(4.8) 


One can also check that under the U(l) transformations, 

g -» gh, h = e ttk2 + 2k0k +2k 

4 f m\g) transforms as 


^j n) ( 9 ) ->• = ^J?(0)eap( -i 


nk 


\J 2 k (k + 1 ) ) 


fc 2 +2fc 


(4.9) 


In terms of geometric quantization, equation (5.3) corresponds to the po¬ 


larization condition on a prequantum 17(1) bundle. The Hilbert space is 
therefore constructed as sections of the holomorphic 17(1) bundle 4/^ . The 
square-integrability of %n is guaranteed by the orthogonality condition of 


the Wigner V function 93 




RR' Omm'Okk' 

dimR 


(4.10) 


where ^(<7)2^ffc(<7 _1 ), -R denotes the representation of G — SU{k + 1), 
and dfi(g) is the Haar measure of G normalizsd to unity; J dfi(g) = 1. The 
orthogonality condition of our interest is given by 


d MVn”${9)Vf% 0 lk{9) = 


m' ,N k 


N k ' 


(4.11) 


Symbols and star products 

We define the symbol of a matrix operator A ms (m , s = 1, 2,, N^) on the 
Hilbert space of fuzzy CP k by (94) 


<i> Y. V ^9)A ms V<^\g) 

ms 

= {(n,0),NM\g T Ag*\(n,0),NM) 


(4.12) 
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The star product of fuzzy CP k is defined by 


(AB) = (A) * (B) 


(4.13) 


From (4.12) can be written 


(AB) = (4.14) 


msrr'p 

We use the relation 


5>4”-°>(s)d" p 0> ( 9 ) = ^ 


From (4.11), the trace of a matrix operator A can be expressed as 

TrA = ^A mm = Nk J dfi(g)V^ (k) A mm ,V^ 


.,*(n,0) 

’. JV(fc) 


= A' ( ‘ ) / diJ,(g)(A). 


The trace of the product of two matrices A, B, is also given by 


(4.15) 


(4.16) 


Tr AB = JV<*> / dn(g){A) * (B) 


(4.17) 


Algebraic construction 


Here we present the construction of fuzzy CP k (k=l,2,...). The coordi¬ 
nates Qa of fuzzy CP k can be defined in terms of La which are x N^- 


rnatrix representation of SU(k + 1) generators in the (n, 0)-represent at ion 94 


Qa 



where C k is the quadratic Casimir. With 


Q aQ A l 

dAB cQaQb = Ck,nQ C i 


(4.18) 


(4.19) 

(4.20) 
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where cIabc is the totally symmetric symbol, and 1 is the x iV( fc ) identity 
matrix. For determine the coefficient Ck. n in (4.20), we notice 


Aa = a\{t A ) 


ijdji 


(4.21) 


A a is the SU{k + 1) generators in the (n,0)-representation, t A SU[k + 1) 
generators in the fundamental representation, with normalization 


tr (t A t B ) = 


25 ab ’ 


(4.22) 


and a\,ai are the creation and annihilation operators. Using the completeness 
relation 


1 


(t a) ij (t a) kl ni&il&jk 


k + 1 


5p' Ski, 


(4.23) 


and the commutation relation [aj, a\] = Sij , we can check A a A a = c!f > . We 
hnd (94j|95) 


(Iabc^-bAc = (k — 1) 

= (k~ 1) 


n 


n + 1 


+ 2^ a \^A)ij a 3 


£ T » S 


Representing by La , we can determine the coefficient 




(Ar-l)/ 


n 


C, 


( k ) 


U+i 


+ 2 '• 


(4.24) 


(4.25) 


Wher C^ 4 is the quadratic Casimir for S'U(fc + 1) in the (n,0)-representation 


(it) _ nk{n + k + 1) 


Ct ; = 


(4.26) 


2 (k + 1) 

4.1.2 Construction of fuzzy S 4 

Introduction to fuzzy S' 4 

The idea of fuzzy S 2 has been one of the guiding forces for us io investigate 
fuzzy spaces 92 . For example, CP fc are successfully constructed in the same 
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spirit as the fuzzy S 2 . Those interested in this area, it is of great interest to 
obtain a four-dimensional fuzzy space. CP 2 is not qualified for this purpose, 
since CP 2 does not have a spine structure 96 . S 4 is well motivated, since 
naturally leads to R 1 at a certain limit. 

There have been several attempts in 97,98 to construct fuzzy S 4 from 
a field theoretic. In 99 the construction is through a projection from some 


matrix algebra and, owing to this forcible projection, it is advocated that 
fuzzy S 4 obeys a nonassociative algebra. 

fuzzy S 4 is alternatively considered in a way of constructing a 


In 98 


scalar held theory on it, based on the fact that CP 3 is a CP 1 (or S 2 ) 
bundle over S 4 . Note that the term fuzzy S 4 is also used, mainly in the 
context of M(atrix) theory, e.g., in 100 also showing that the existent fuzzy 
S 2 and S 4 models are natural candidates for the quantum geometry on the 
corresponding spheres in AdS/CFT correspondence. Also the construction of 
fuzzy S 4 is considered through fuzzy S 2 xS 2 . This allows one to describe fuzzy 
S 4 with some concrete matrix configurations 101 . However, the algebra is 


still non-associative and one has to deal with non-polynomial functions on 
fuzzy S 4 . Since those functions do not naturally become polynomials on S 4 
in the commutative limits. In the case of fuzzy CP A , the fuzzy functions 
are represented by full (N x Y)-matrices, so the product of them is given 
by matrix multiplication which leads to associativity of the algbra for fuzzy 
CP /v . The extra constraint is expressed as an algebraic constraint such that it 
enables us to describe the algebra of fuzzy S 4 in terms of the algebra of SU (4) 
in the (n,(^-representation, the algebra of fuzzy S 4 is obtained from S U (4) 
as well with the extra constraint on top of these fuzzy CP 3 constraints. 


102 


Construction of fuzzy S 4 

We begin with construction of fuzzy CP 3 . The coordinates Qa of fuzzy 

CP 3 

Qa = ~^=, (4.27) 

\/cP 

where La are x A^T-matrix representation of SU(4), and the coordinates 
satisfy the following constraints: 

QaQa = 1, (4.28) 

^ABcQaQb = Cs^nQc- (4.29) 
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Where 1 is the N ® x N^ identity matrix 

N (3) = -{n + l)(n + 2 )(n + 3), 
6 


and from (4.25) 


C 3 , n = 


2 (n 1 
- + - 




(3) V 4 2 


(4.30) 


(4.31) 


t ?2 is the quadratic Casimir of 57/(4). From (4.26) 


(3) = 3 n(n + 4) 
2 8 

for k <C n (in this case k — 3), we have 


C\ 


k,n 


C-k 


k(k + l) 


(k- 1 ) 


(4.32) 


(4.33) 


57/(4) is decomposition as, 517(4) —» 517(2) x 517(2) x 17(1), where 517(2) 
and U (1) are defined 


517(2) 

0 


0 

5*7(2) 


0 


Each 5*7(2) is the (2 x 2)-matrix representation. 

In order to obtain functions on fuzzy 5 4 , we need to require 

[T 7 , L a ] = 0 


(4.34) 


where T 7 denote matrix-functions of Q^s and are generators of H repre¬ 
sented by N^ x A7 3 )-matrices, H being relevant to the above decomposi- 
tionof SU(4) (77 = 57/(2) x 7/(1)). Construction of fuzzy 5 4 can be carried 


out by imposing the additional constraint (4.34) onto the functions on fuzzy 
CP 3 . What we claim is that the further condition M.34| makes the func¬ 


tions J-[Qa ) become functions on fuzzy 5 4 . This does not mean that fuzzy 
5 4 is a subset of fuzzy CP 3 . Qa(A = 1,..., 15)’s are defined in R 15 . While 
locally, say around the pole of A = 15, globally they are embedded in R 15 . 
The equation (4.34) is a global constraint in this sence. So the algebra of 


fuzzy 5 4 is given by a subset of SU(4). The emerging algebraic structure 
of fuzzy 5 4 will be clearer when we consider the commutative limit of our 
construction 102 . 


78 














4.1. CONSTRUCTION OF FUZZY SPACES 


A formula for the large N-limit of matrices 


In this subsection, following 103,104 

Lb A, where Lb are the generators of SU(K + 1) 


We now consider the symbol for 
in this case 


the product 

S' 4 , k=3), A being an arbitrary N^ x TV^-matrix. From (4.12), the symbol 
of Lb A is given by 


(L b A) = ((n,0) 1 NW\g T L B Ag*\(n,0),NW) (4.35) 


Commutative limit 

In the large n limit we can approximate Q A to the commutative coordi¬ 
nates on CP 3 ; 


Qa~ 0A = -2 tr(gH A gt 15 ) 

With the constraints for cp 3 according (4.28) and (4.29) is 


4>A<t>A = 1, d A BC <t>A <t>B = \l 2</>C 


(4.36) 


(4.37) 


In (4.36), t A is the generators of S'f/(4), and g is a group element of S'f/(4) 
given as a (4 x 4) matrix 


Matrix-function correspondence 

We studing our construction of fuzzy S' 4 by confirming its matrix-function 
correspondence. We focus on two things: 


(i) Linking between the number of matrix element for fuzzy S' 4 and the 
number of truncated functions on 4-spheres S 4 . 

(ii) The relationship between the product of functions on fuzzy S' 4 and that 
on 4-spheres S 4 . 

In the case of fuzzy S' 2 , S' 2 = SU(2)/U(1). We consider f>|£(f/) is Wigner 
77-functions for SU(2), where the spin-) representations of an SU(2) group 


element g and ( m,n = — j ,...,+_)), according (4.3) 

v mn(9) = (jm\g\jn) 


(4.38) 
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Functions on S 2 can be expanding when n = 0 in terms of particular Wigner 
©-functions U^ 0 (g), which are invariant under a 17(1) right-translation op¬ 
erator acting on g. For definition (4.4). The state |j0) has no 17(1) charge, 
right action of the 17(1) operato, R 3 V^ 0 (g) = 0, the ©-functions basically 
are spherical harmonics 


vl(g) = 


47T 


21 + 1 

and truncated expansion can be written as 

n l 






(4.39) 


(4.40) 


1=0 m=—l 


The number of f l m are computing by Yfi=o(^ + !)• 

This relation implements the condition (i) by defining functions on fuzzy 
S 2 as (n + 1) x (n + 1) matrices. The condition (ii). One can show an exact 
correspondence of products. Let f mn (rn, n = 1, ..., n + 1) be an element of 
matrix function-operator / on fuzzy S 2 . As (4.3) 


(/} = Y (s)V%(g), 

m,n 

where V^fMg) = ©i£(<7 -1 )- The star product of fuzzy S 2 is 


(4.41) 


(fa) = (/) * (a) 


(4.42) 


from (4.57) we can write 


</j} = Y f mn g tl I>5 l (s)I>»to)+“(9)©,“(9) (4.43) 

m,n,k,r,l 

with using the orthogonality of ©-function ©nr (d^krid) — Snk- Let R- 
the lowering operator in right action 


R - V mn(a)V (j + n)(j-n+ l)V%n-l(a) 
by iteration (4.43) [102| 

</9> = E(-ir^7^r^-(f)K(g) 


s=0 


s!(2 j)! 


(4.44) 


(4.45) 
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where we use the relation R* = R + . In the large j limit, the term with s=0 
in (4.45) dominates and this leads to an ordinary commutative product of 
(/) and (g). By construction, the symbols of functions on fuzzy S 2 can be 
regarded as commutative functions on S 2 

The matrix-function correspondence for fuzzy CP 3 can be expressed by 


N (3) x 1V (3) = E dim ( 1 ’ 0 (4.46) 

1=0 

dim(l, l) = ^(2/ + 3)(/ + l) 2 (/ + 2) 2 (4.47) 

where dim(l,l ) is the dimension of SU(4) in the (l, ^-representation. This 
expression indicates that the number of matrix elements coincides with the 
number of coefficients in an expansion series of truncated functions on CP 3 = 
SU(A)/U(3 ) can be ewpanded by V d } l ^(g), Wigner-P-functions of SU( 4) in 
the (l, ^-representation (l = 0, 1, 2,...). The lower index M(M = 1,..., dim(l, l ) 
labels the state in the (l, l) -represent at ion, while the index 0 represents any 
suitably fixed state in this representation. 

Return to the conditions (i) and (ii), we trying 

1. Counting the number of truncated functions on S' 4 

2. A one-to-one matrix-function correspondence for fuzzy 4-sphere S' 4 

3. Proposing a block-diagonal matrix realiazation of fuzzy S 4 


For (1), can be made in terms of the pherical harmonics Y/p 2 z 3m on S 4 


with a truncation at l\ = n 105 


n l\ I 2 

+ 1) = ^(n + l)(n + 2) 2 (n + 3) (4.48) 

Zi =0 Z2—0 ^3=0 

Can be regarded as an N 2 (Z)-degeneracy due to an S 2 internal symme¬ 
try for the extraction of S 4 out of CP 3 ~ S 4 x S 2 


102 


_ . The number 

of truncated functions on CP ' 1 is given by (4.47), the number of those 
on S 4 may be calculated by 


N s ‘(n) = Y, 


1=0 


dim(l, l ) 
N 2 (l) 


E + 4 )(^+ 2 )( 2 ^+- 1 -) 


6 
1=0 

= ^(n + l)(n + 2) 2 (n + 3) (4.49) 
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• For (2). Let ( F)u(I,J = 1, 2, N^) denote a matrix-function on 
fuzzy CP 3 , and can be defined by (I\F\J) as matrix element of the 
function F on fuzzy CP 3 , also we denote = |/), where 

0ii ...i n is the staate of fuzzy CP" 

We need to find an analogous matrix expression (F s ^)jj for a function 
on fuzzy S' 4 , for this, we splitting each i of the 0i 1 i 2 ...i Tl into a and a as 




= {<A 


a\a 2 ..-CLn ? Ycb\a\...a n —\’) ••••> Y / ai...a n _iai ? '+'0,10,2 


...aj. (4.50) 


can obtain the states corresponding to fuzzy S 4 by imposing an addi¬ 


tional condition on (4.50) ,i.e., the invariance under the transformations 
involving any a m (m = 1 ,,n). On the set of states (j>a 1 ...a n , which are 
(n + 1) in number, the transformations must be diagonal 102 , but we 
can have an independent transformation for each state. The number 
of the states is (n + 1), since the sequence of a m = {3,4} is in a totally 
symmetric order. Thus we get (n+1) different functions proportional 
to identity. On the set of states (j) ai a 1 ...a n _ 1 , we can transform the a± 
to b\ = {1,2}, corresponding to a matrix function / oli6l which have 2 2 
independent components. But we can also choose the matrix f ai ,b! to 
be different for each choice of (di...a n _i) giving 2 2 x n function in all, 
at this level. We can represent these as . Where we split i m 

into a m , a m and j m into b m , b m . We find that the set of all functions 
on fuzzy S A is given by 


(F S I = { . . /■(“ l-“n-l)f . . 

V JU \j w ai...d n ,bi...bn'> J a 1 ,b 1 w a 1 ...d rl ,bi ...b n i 

fa^a 2 ,bib2^a 1 ...d n -2,b 1 ---bn-2’ . ’ f a 1 ...a n ,b 1 ...b n } • (4-51) 


The structure in (4.51) shows that F s> is composed of (l + 1) x (l + 1)- 
matrices. The number of matrix elements for fuzzy S' 4 is counted by 


N s \n) = J^(Z + l) 2 (n + l —Z) (4.52) 

1=0 


with the number of these matrices for fixed l being (n + 1 — /) and 
(/ = 0,1,.., n). 
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The symbol of the function F On fuzzy CP 3 can be defined as 
(F) = ^^ i (N\g\I)(F) lJ (J\g\N) 


(4.53) 


|IV) = |(n, 0), N 3 ) is the highest or lowest weight state of fuzzy CP 3 , 
(J\g\N) denote the previous P-function. 

We now consider the product of two functions on fuzzy S 4 , a function 
on fuzzy S 4 can be described by (Z + 1) x (Z + l)-matrices. The symbol 
of a function on fuzzy S 4 is d efined in the same way except that ( F)ij 
is replaced with (F s ) u in (4.53). The star product of fuzzy S 4 is 
written as 


(F^G^) = 

IJK 


)u(G s ) JK (N\g\I)(K\g\N) 


(4.54) 


We can describe CP 3 in terms of four complex coordinates Zj with the 
identification Zj ~ AZj, (A G C — {0}). Following 106 , write Zj 


Z% TTd) (-^ad^"d) a ) 


(4.55) 


where oo,7r is tow spinors, a = 1,2, a = 1,2 and x a a can be dedined 
with coordinate x M on S 4 , x a d = (1^4 — iViX 1 ), (i = 1,2,3), <7* beign 
(2 x 2) Pauli matrices. 

The homogeneous complex coordinates of CP 3 are defined by Z = 
(zi, z 2 , z 3 ) T , T is transposition of the vector or (1 x 3)-matrix. We 
introduce the notion U a = g at4 , g a ,A{oi = 1,2, 3,4) is matrix elements 
are defined the coordinates on CP 3 with UU = 1, Ua s are related to 
Z 


U a = , 

\/l + z.z 


(4.56) 


We can parmatrize U by the homogeneous coordinates Z i: i.e, Ui = 
-4=. Function on S 4 can be considered as function on CP 3 which 

y/Z.Z 

satisfy 


k fcp 3 (Z, Z) — / CP 3(Z, Z) — 0 

UTTa O'Ka 


(4.57) 
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For (3). Let us write down the equation (4.52) in the following form: 


N s Un) = 


+ 1 + 2 2 
+ 1 + 2 2 + 3 2 


+ 1 + 2 2 + 3 2 + 4 2 + ... + (n + l) 2 . (4.58) 

Coordinates of fuzzy S' 4 are then represented by these N^ x N ^ block- 
diagonal matrices, X^ , which satisfy 

XaXa ~ 1 (4.59) 
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The Algebra for Non-commutative 5 4 

At the beginning we will assume 50(5) invariance of the algebra, the 
algebra of the coordinates will be given 

[X a ,X b ] = e abcde X c X d X e (4.60) 

X is (N x N ) hcrmitian matrice and e a b c de is the Levi-Civita symbol, with 
the condition 


(. x a y -C 2 1 = 0 


(4.61) 


The algebra (|4.60|) can be derived as equations of motion from the following 
action 
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5 = / dt Tr<j ±p R (D 0 X a ) 2 + X b ] - ae abcde X c X d X e ) 2 *> (4.62) 


In fact the bosonic part of this action for Matrix theory in the pp-wave 
background has this form |107|. This may be regarded as the bosonic part of 
the M-atrix theory in some background with cencel other bosonic cordinates, 
A" 6 ,..., X 9 , here X a, s (a = 1,..., 5). Ho is the covariant derivative dt + z[A 0 , ] 
and /3r is plaing role the raduis of the R circle. 

A representation of (4.60) in terms of 4 x 4 matrices is given by tensor 
products of Pauli matrices. 


-Aq = -cr 3 <g> err 

-2 1 

X 0 = -cr 3 <g> cr 2 

^ 1 

x 0 = 3^3 ® <?3 
Aq = ~o"i ® 1 2 

X 0 = ^2 ® 12 (4.63) 

1 2 is 2 x 2 identity matrix. We can replacing one of the two Pauli matrices 
by a spin-j represent art in of 50(3), TU (These are four dimensional gamma 
matrices. It is then N = 2(2 j + 1) dimensional representation). When 
satishes 


[Tu),T^]=ze abc T^ 


(4.64) 
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(4.65) becomes 


*o = 

3^3 ® T^) 

>0 
o to 

\° 3 ® T U) 

*o = 

^ 3 ® T (J) 

X 4 o = 

1 

-CTl 0 l 2 j+l 

X 5 0 = 

-cr 2 0 lgj+i 


All Xg’s satisfy (4.60) except 


[X 4 ,X 5 ] = 3 e 45cde X c X d X 

4j(j + 1 ) 

We will construct a non-commutative product on S 
Matrix configuration in terms of the product on S 2 . 


(4.65) 

e (4.66) 

4 corresponding to the 






4.1. CONSTRUCTION OF FUZZY SPACES 


4.1.3 The Fuzzy sphere S% 

Every manifold comes with a naturally defined associative algebra of func¬ 
tions with point-wise multiplication. This algebra is generated by the coor¬ 
dinates of the manifold and is from the definition commutative. As it turns 
out, this algebra contains all the information about the original manifold 
and we can describe geometry of the manifold purely in terms of the algebra. 
Also, every commutative algebra is an algebra of functions on some manifold, 
therefore, what we get is 

commutative algebras <—> differentiable manifolds 

A natural question to ask is whether there is a similar expression for non- 
commutative algebras, or 

non-commutative algebras <—> ???? 

Quite obvious answer is no, there is no space to put on the other side of the 
expression. Coordinates on all the manifolds commute and that is the end of 
the story. So, as is often the case, we define new objects, called noncommu- 
tative manifolds, that are going to fit on the right hand side. Namely we look 
how aspects of the regular commutative manifolds are encoded into their cor¬ 
responding algebras and we call the non-commutative manifold object, that 
would be encoded in the same way in a non-commutative algebra. 

This is going to introduce non-commutativity among the coordinates. 
This notion should not be completely new, as the reader probably recalls the 
commutation relations of the quantum mechanics 

[x\pj] = ih6j (4.67) 

In classical physics, the phase space of the theory was a regular manifold. 
However in quantum theory we introduce non-commutativity between (some) 
of the coordinate and therefore the phase space of the theory becomes non- 
commutative. One of the most fundamental consequences of the commuta¬ 
tion relations is the uncertainty principle. The exact position and momentum 
of the particle can not be measured and therefore we can not specify a single 
particular point of the phase space. Similarly, if there is non-commutativity 
between the coordinates, there is a corresponding uncertainty principle in 
measurement of coordinates. The notion of a space-time point stops to make 
sense, since we can not exactly say, where we are. This introduces a short 
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distance structure to the space, quantities at short distance are not well lo¬ 
calized. Because of this the non-commutative space are sometimes called the 
fuzzy spaces. 

In practice, we often deform a commutative space into its non-commutative 
analogue. In this way we get noncommutative spaces that give a desired com¬ 
mutative limit. An example of such deformation is already mentioned phase 
space of quantum mechanics, but to illustrate the idea better, and since we 
will need the notions later, let show how this works for a two-sphere. 

The fuzzy sphere Sjj is a matrix approximation of the usual sphere S 2 . 
The algebra of functions on S 2 , spaned by the spherical harmonic, is trun¬ 
cated at a given frequency. The algebra then becomes the finite dimensional 
algebra of NxN matrices. More precisely, recall that the algebra of functions 
on the ordinary sphere can be generated by the coordinates of R 3 modulo the 
relation 
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E 

Z —1 


X- 


= R 2 


(4.68) 


This comes with an understood condition on commutativity of the coordi¬ 
nates 


XiXj — XjXi = 0 (4.69) 

Coordinates constrained in this way generate the algebra of all the functions 
on the sphere. Note that this is technically not the easiest way to do so. 
It is easier to introduce only two coordinates 9, 0 on the sphere and define 
the algebra of functions not by the generators, but by the basis, e.g. the 
spherical harmonics. However the two sphere defined in our way is easier 
deformed into the non-commutative analogue. Now we define the fuzzy two 
sphere by the coordinates ah, which obey the following conditions 

= ( 4 - 7 °) 

i= 1 

p the radius of the non-commutative sphere. R did describe the regular 
sphere 


rp . rp . _ rp . rp . 

•Xj ^ J iL/J 




(4.71) 
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The radius of the original sphere was encoded in the sum of the squares of 
the coordinates. We see, that such x's are achieved by a spin-j representation 
of the SU(2). If we chose 

3 J\j2 _ i 

y>, 2 =j(j + 1) = —J—, (4.72) 

i =1 

with 

[Li, Lj] = ietjkLk (4.73) 


x = A N Li, i = 1, 2, 3 


(4.74) 


Ajv is the noncommutative parameter of dimension length 

A - 2p 

An ~7W^T 

with N the dimension of the representation. Matrices x t become coordi¬ 
nates on the non-commutative sphere. The limit N —y oo removes non¬ 
commutativity, since A —* 0, nd we recover a regular sphere with radius 
R 

The algebra of functions Sjj therefore coincides with the simple matrix 
algebra Mat(N, C). The normalized integral of a function / e is given 
by the trace 

l * = ^ 4 ‘ 76 ' ) 

Js% N 

The functions on the fuzzy sphere can be mapped to functions on the com¬ 
mutative sphere S 2 using the decomposition into harmonics under the action 


(4.75) 


Jif=[U,f\ (4.77) 

of the rotation group SU (2) . One obtains analogs of the spherical harmonics 
up to a maximal angular momentum N — 1 . Therefore Sjj is a regularization 
of S 2 with a UV cutff, and the commutative sphere S 2 is recovered in the 
limit N —> oo . 
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4.1.4 The fuzzy 2-sphere Sj Vl x Sj Vr 

The simplest 4-dimensional generalization of the above is the product 
5f x Sjf of 2 such fuzzy spheres, with generally independent parameters 


Nl,r ■ It is generated by a double set of representations of su( 2) commuting 
with each other, i. e. by Af, Af satisfying 


The Casimirs is difines 


[Af, Xf] 

[Af, Af ] = 
[Af, Af] = 


£ AfAf = 

i =1 


Rijk X 


k > 


itijkXf, , 


[Af, Af] = 0 


iVl-1 


iV|-l 


i=l 

This can be realized as a tensor product of 2 fuzzy sphere algebras 

Af = Aj <g> Ia^xa^, 

•Af = 1 jv k xat h Aj. 

The normalized coordinate functions are given by 


L,R 

X- ' = 


2 R 


with 


E( 


f (N L ’ R ) 2 - 1 

E 


A 


L,R 


xh 2 = 


[X 


R\ 2 


= R 


(4.78) 


(4.79) 


(4.80) 

(4.81) 


(4.82) 


(4.83) 


In rinciplc one could also introduce different radii R l - r for the 2 spheres, but 
for simplicity we will keep only one scale parameter R (and sometimes we will 
set R=l). This space 1 can be viewed as regularization of S 2 x S 2 C R 6 , and 
admits the symmetry group SU( 2) L x SU(2) R C 50(6). The generators xf :R 
should be viewed as coordinates in an embedding space R 6 . The normalized 
integral of a function / e 5f L x 5f fl is now given by 



16vr 2 R 4 


M/), 


(4.84) 


where we difine the 167r 2 i? 4 as volume 
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4.1.5 The limit to the canonical case Rq 

We will mainly consider Nl = Nr. It is well-known [?] that if a fuzzy 
sphere is blown up near a given point, it can be used to obtain a (compacti- 
hed) noncommutative plane with canonical commutation relations: Consider 
the tangential coordinates x\^ near the north pole X 3 = R . Setting 

R 2 = N9/2 (4.85) 

they satisfy the commutation relations 

op 9 7? /- 

[xi,x 2 ] = i-jy x 3 = i-^\JR 2 -x\-x\ = i9 + 0(1/N) (4.86) 

Therefore in the double scaling limit with N, R —> 00 keeping 9 fixed, we 
recover the commutation relation of the canonical case, 


[x\,x 2 ] = i9 


(4.87) 


up to corrections of order 1/N. Similarty, starting with Sjj L x and setting 


we obtain in the large Nr , Nr limit 



r9l,r 

2 

(4.88) 

i€ij9 L , 

(4.89) 

ieij9 R , 


[xf, xf] = 0 



This is the most general form of Kg with coordinates (xi,..., xR) = (xf , x %, xf, xf ). 
The integral of a function /(x) then 



xS% 

n R 


f{x) ->• 4vr 4 0 L d R tr((/(x)) =: 



(4.90) 
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4.2 Fuzzification 

V CF 2 acts on a subspace of Acp 2 We can thus conceive of a 

fuzzy Dirac operator D C pa which acts on a subspace of A C p 2 <S>Mati 6 , A C p2 
being obtained from Acpa by restricting “orbital” SU(3 ) IRR’s to ( n,n ), 
n<N. D C p 2 is then obtained from T? C p 2 by projection to this subspace. 
T > cP 2 commutes only with the total 51/(3) Casimir J 2 and not with orbital 
SU(3 ) Casimir £ 2 . This causes edge effects distorting the spectrum of D CP 2 
for those states having (n, n) near (N,N) which T >cpz mixes with (n ,n ), 
n >N. This particular edge phenomenon does not occur for S 2 = CP 1 where 
orbital angular momentum C 2 commutes with the Dirac operator. A way to 
eliminate such problems is suggested by the work of 109 -112 : We introduce 


the cut-off not on the orbital Casimir, but on the total Casimir, retaining 
all states upto the cut-off. That seems the best strategy as it will give a 
fuzzy Dirac operator Dcp 2 with a spectrum exactly that of the continuum 
operator X>cp 2 upto the cut-off point, and which has chirality (chirality r C p 2 
of T>cp '2 commutes with Jf) and no fermion doubling. This approach is the 
same as the method adopted for S 2 in 109-112 . For S 2 , the edge effect 
turned up as the absence of the —E eigenvalue subspace for the maximum 
total angular momentum when the cut-off is introduced in orbital angular 
momentum, and attendant problems with chirality. 

4.2.1 Coherent States and Star Products : The Case 
of S 2 ~ CP 1 


These have been treated in 1111,113,114 . Here we summarize the main 


points so that we can outline the relation of wave functions and those based 
on matrices for fuzzy physics. 

Let us first consider S 2 = CP 1 and its fuzzy versions. The algebra A is 
Mat 2 i+i■ SU( 2 ) acts on A on left and right with generators Lf and — Lf, 
and orbital angular momentum is = Lf — Lf. The spectrum of C 2 is 
K(K + 1), K = 0,1,.., 21. We can find a basis of matrices T/} diagonal in C 2 
and £ 3 (with eigenvalue M) and standard matrix elements for Hi. A acts on 
a (21 + 1 )—dimensional vector space with the familiar basis \l,m). T^ r are 
orthogonal, K(K + 1) and M being eigenvalues of £ 2 and £ 3 : 


( t m, T m') ■= TrT M T M' = constant x 5 KK >5 MM >. (4.91) 
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The above suggests that there is a way to regard A as “functions” on S 2 
with angular momenta cut-off at 21. Such functions are also represented by 
the linear span of spherical harmonics Y KMl K<21. We want to clarify the 
relation of Ykm'h to the matrices T^ r in A. 

Towards this end, let us introduce the coherent states 

\g) = U®(g)\l,l) (4.92) 

induced from the highest weight vector |/; l). g^U^ K \g) is the angular mo¬ 
mentum K IRR of SU( 2). Note the identity 

I ge^ e ) = e iW \g). (4.93) 


It is a theorem 113 that the diagonal matrix elements (g\a\g) completely 
determine the operator a. Further (ge l ^~ e \a\ge l ^~ 9 ) = (g\a\g) so that (g\ci\g) 
depends only on 


g°zg 1 = 


a ■ x, x 2 = 1; xGS 2 . 


(4.94) 


i= 1 


In this way, we have the map 


A—»C°°(S 2 ), 
a — yet, 

where 

a{x) = < g\a\g). 


(4.95) 


In this map, the image of is Ykm after a phase choice: 


Ykm(x) = {g\T*\g). 


(4.96) 


For, under g—>hg , x^R(h)x where h—*R(h) is the SU( 2) vector represen¬ 
tation. Under this transformation, since 

K 

Y KM (R(h)x) = ^ D^ K \h) MM > Y km > (x) (4.97) 

m'=-k 

and 

K 

T’‘^HrW(h)- 1 Tj‘uW(h)= W D IK \h) MM ,T«, (4.98) 

M'=-K 
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where h — vD^ K \h) is the angular momentum K IRR of SU( 2) in a matrix 


representation, we have the proportionality of the two sides. (4.96) and phase 
conventions fix the constant of proportionality. 

The map T^^Ykm is an isomorphism at the level of vector spaces. It can 
be extended to the noncommutative algebra A by defining a new product on 
Ykm’s, the star product. Thus consider (g\T^T^\g). The functions Ykm and 
Yin completely determine T^ f and Ty, and for that reason also this matrix 
element. Hence it is the value of a function Ykm * Y LNl linear in each factor, 
at x: 

(g\T«Tk\g) = [Y km * Y LN ](x). (4.99) 

The product * here , the star product , extends by linearity to all func¬ 
tions with angular momenta <21. The resultant algebra is isomorphic to the 
algebra A. 

m (see also 


’lie explicit formula for * has been found by Presnajder 


114 ). The image of /ha is just —i(x AV)*a. We will use the same symbol C t 
to denote —i(xA'V)i derivation. The * product is covariant under the 577(2) 
action in the sense that 


Ci(a * b) = (Cid) *b + a* (Cjb). 


(4.100) 


It depends on l and approaches the commutative product of C°°(S 2 ) as 
l — kx >. Coherent states thus give an intuitive handle on the matrix rep¬ 
resentation of functions. 

But on S 2 , we also have monopole bundles. Sections of these bundles for 
Chern class n are spanned by the rotation matrices Dmh , j > |n|. They have 
the equivariance property 


DSlse 4 ?") = (sV"". 


(4.101) 


This last equation is essentially a generalization of equation (4.100) , in other 
words one can identify Dmn with < j, m\^n > where |>= D^\j, n > [see 


equation (4.101) 


How do we represent them by matrices? 

In the first instance, let n >0 and consider the coherent states (now with 
an additional label) 


| g]l + n) = U« +n \g)\l + n,l + n) 
\g;i ) = u (l \g)\ij). 


(4.102) 
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They span vector spaces Vi +n and Vi. We can consider the linear operators 
Hom(Vi +n ,Vi ) from Vi +n to V/. They are [2Z + 1] x [2(Z + n) + 1] matrices 
in a basis of VJ+n and Vi, and have U^ l \g) acting on their left (with gener¬ 
ators Lf) and U^ l+n \g) acting on their right (with generators —Lf). We 
can decompose Hom(Vi +n , Vi) under the “orbital” angular momentum group 
(with generators C t = Lf — Lf) into the direct sum ©f-l"(A") 
with the IRR K having the basis Tff, with C 3 Tff = ML As before, we 
choose Tff so that £, follow standard phase conventions. Tff are orthogonal 

Tr ( T M’) ]T M = constant xS K ’ K 5 M > M . (4.103) 


Now consider 


(g-,l\T«\g;l + n). (4.104) 

It transforms in precisely the same manner as (g) under g^rhg and 
g^-ge 1 ^ 6 / 2 and hence after an overall normalisation, 

{g;l\T^\g;l + n) = D^(,g). (4.105) 


Thus Hom(y +n , VJ) are fuzzy versions of sections of vector bundles for Chern 
class n> 0. For n < 0, they are similarly Hom(Vi,Vi + \ n \). This result is due 
to 115| (see also 109-112,116)). An explicit formulae for the fuzzy version 
of rotation matrices can be found in 


111 


It is interesting that Chern class has a clear meaning even in this matrix 
model: It is |V| — |W| for Hom(V, W), where |V| and |W| are dimensions of 
V and W. 

There are two (inequivalent) fuzzy algebras acting on Hom(V, W). Mat\y\ 
A|v| acts on the right and Mat\w\ = A|w| acts on the left, where now a sub¬ 
script has been introduced on A. These left and right actions have their own 
*’s, call them *|y| and *|w|: if Ay, b&A w and a and b are the correspond¬ 
ing functions, then 


bT 


K n- 

M u 


>b *\ W \ Y) 


\w\ * km *|y| a 


(4.106) 


under the map of Hom(V, W ) to sections of bundles. There is also a fuzzy 
analogue for tensor products of bundles. Thus we can compose elements of 
Hom(V, W) and Hom{W,X) to get Hom(V,X) 


HomiV , X) = Hom(V, W)® Am Hom{W, X). (4.107) 
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Its elements are ST, S^Hom(V,W), TEHom(W, X). Its Chern class is 
\V\ — |A|. If S and T are the representatives of S and T in terms of sections 
of bundles, then ST—>S * T. 

Tensor products of two vector spaces Ti and T 2 over an alge¬ 

bra B are defined only if ri(r 2 ) is a right-(left-) T-module JTTT . Hence 
Hom(V, W)®A ]w \Hom(W', X) is defined only if W = W'. So S*T is rather 
different in its properties from the usual tensor product of bundle sections, 
in particular T * S makes no sense if V ^ X. 


4.2.2 Fuzzy Dirac Spinors on Sf. 


We can now comment on the fuzzy . Elsewhere the Watamuras 118,119 


and following them, us 120,121 , investigated the Dirac operator as acting 
on A<g)C 2 = A 2 , A = Mat^i+i- That led to rather an elaborate formalism 
because of the cut-off in orbital angular momentum. So as indicated earlier, 
it seems more elegant to cut-off total angular momentum at some value jo- 


We can now argue such a cut-off leads to the formalism of 109-112,122 


and to supersymmetry. Thus let T^ + ef7om(V) + i/ 2 , V)) with the transfor¬ 
mation property < g;l\T^ + \g;l + | > — >e 1 ^ < g\ Z|T^ + |p; l + \ > under 
g — >ge l ^ a3 . One also has the transformation property 


U {l \g) ] T! n+ U^\g) = V 


(4.108) 


[So j<2l + 1/2 and jo — 21 + 1/2]. Then one can make the identification 


DU(g) = {g\l\T’ m+ WJ + -), 


(4.109) 


since from equation (4.101) it is easy to see that D J m+ (g )— D J m+ (g) under 
g — yge 1 ^ 3 . 

The subscript + in Tij n+ indicates helicity — , i.e T/j + is the fuzzy version 
of < j,Tn\ij}+ > of (??) so that it will be associated with the negative helicity 
part of the wave function. 

For helicity +, but for same j 0 , we have to consider T/ t _e Horn (V), 1 / 2)5 

with 

C <l+ +9) t T+£/ (,) (9) = E+L'(9)+_- ( 4 -H0) 


96 
















4.2. FUZZIFICATION 


Of course now 


D 3 m -(g ) =< g-,l + -\T 3 m _\g-,l >, 


(4.111) 


where both sides will acquire now a phase exp(—i |) under the right 1/(1) 
action, namely under g 

of < j, m\ip J 2 > . 

This is the formalism of 


gexp(i^<Jz) . T^_ is then clearly the fuzzy version 


109 


112 


122 


As we have united V® and 
U( z+1 /2), it is natural to consider OSp( 2,1) or even OSp( 2,2) SUSY as dis¬ 
covered first by Grosse et al in the second paper of (Il2|. 

Because of the mixing of l and l + 1/2, we have to reconsider the ac¬ 
tion of the matrix algebra A approximating A = C°°( S 2 ). Mat 2 i+i acts 
on T^ + (T^_) on the left (right) while Mat 2 i +2 acts on T^ + (T^_) on the 
right (left). So it is best to regard fuzzy functions to act on left (say) of 
T^ + and right of T 3 m _ as Mat 2 i+i■ This suggestion is slightly different from 
that of 109-112,122 where they regard the fuzzy algebra to be Mat 2 i + \ on 


T 3 1+ and Mat 2 i +2 on T^_, both acting on left. However, our proposal does 

not generalize to instanton (monopole) sectors. 

We can restore spin parts to fuzzy wave functions. The spin wave func- 

1 

tions for helicity ± are T£ s ± , where m s denotes the two components of the 
spinor . The positive chirality spinors are defined by 

< 9'J\ T m 3 +\9',l + 7 \ >= D rn s + =< >> (4.112) 

while the negative chirality spinors are defined by 

<9'J+ 9‘J >= D m s — =< \,m s |^i s) > • (4.113) 

So the two components of the total fuzzy wave functions for helicity ± are 


< 7,,m, \ip$ >= 


Ee. ±T U 


T 

1 1 


j,m 


m s d= ? 




^GC,m s = 


(4.114) 


This is the fuzzy version of equation 

The Dirac operator JD 2g is given by the truncated version of (??) 




-, 1/2 


rUh 






OF 






j<2l + 1/2, 


(4.115) 
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being the angular momentum j images of y. 


4.2.3 The Case of CP 2 

Coherent states for CP 2 can be defined using highest weight states. For 
IRR (3,0), we can pick the highest weight state with / = J 3 = 0, Y — —2/3, 
namely the c—quark: |0,0, —2/3) = |0,0, — 2/3; (3, 0)). Then if g — >U^ 3,0 \g) 
defines the IRR, \g; (3,0)) = U^ 3 ’°^(^)|0, 0, — 2/3; (3,0)). For the IRR (IV, 0), 
we can simply replace 10, 0, —2/3; (3, 0)) by its N —fold tensor product 

2 2 2 2 N 

10, 0, —; (3, 0))<E>|0, 0, —(3, 0))<8)...<E>|0, 0, —(3,0)) = |0, 0, ; (IV, 0)), 

(4.116) 

and set 

9 N 

\9\ (IV, 0)) = U^°\g)\0,0,- — -(N,0)). (4.117) 

For (0, IV), we can use the c—quark state \g; (0, 3)) = U ( '°’ 3 \g) 10, 0, +2/3; (0, 3)) 
and its tensor product states. 

The development of ideas now keep following S 2 = CP 1 . Full details can 
be found in |114 . 

General theory confirms that the maps a—>a from matrices in the (IV, 0) 
or (0, IV) IRR to functions on CP 2 , defined by 

«(0 = ((IV, 0); g\a\g] (IV, 0)) 
or 

«(0 = ((0, IV); <7|a|<7; (0, IV)). (4.118) 


are one-to-one so that a *—product on a’s exists. In this map, the £17(3) 
generators £* acting on a become the corresponding CP 2 £77(3) operators 
—ifijk£j-S~- We shall use the same symbol for these operators too. The 

orbital £17(3) action is compatible with * in the sense that C t (a * b) = 
(£jS) * b + a * ( Cjb ). Irreducible tensor operators of SU( 3) are well studied 
El- With their help, fuzzy analogues of D —matrices can be constructed, 
as also sections of 17(1) and 17(2) bundles. 

The fuzzy CP 2 Dirac operator is the cut-off . We omit the details: the 
necessary group theory is already to be found in 116 while the rest is routine. 









4.3. GAUGE THEORY ON FUZZY S 2 x S 2 


4.3 Gauge theory on fuzzy S 2 x S 2 

Now that we have the fuzzy space S% L x Sjf corresponding to (NlNr) 2 - 
dimensional matrices, we want to construct a matrix model having S^ L x 
as its ground state. As in the canonical case, the fluctuations around this 
ground state will produce a gauge theory. But as the matrices are now finite¬ 
dimensional, the model will be well defined and hnite. 

In the fuzzy case, it is natural to construct S% L x S% as a submanifold 
of R 6 . We therefore consider a multi-matrix model with 6 dynamical fields 
(covariant coordinates) and Bf' , with i run from 1 to 3, which are 
(N^Nr) x (N^Nfi) Hermitian matrices. As action we choose the following 
generalization of the action in 


124 


S= -y 

g 2 


2 B i a jb^ia jb + A V^R 


(4.119) 


with a,b = L, R and i,j=l ,2,3; summation over repeated indices is implied. 
Here (Pl,r are defined as 


• = 1 ( b l b l - i 

P L ‘ jj>2 iW B i 4 

1 N 2 — 1 

. _ 1 (r>RuR U 

VR ■ - B i -4—) 


(4.120) 


and R denotes the radius of the two spheres. The field strength is defined by 

HkBi), (4.121) 


FiL jL = 

1 

R 2 

1 

FiR jR = 

1 

R 2 

FiL jR = 

1 

R 2 


41 






jR\ 


This model (4.119) is manifestly invariant under SU{2)l x SU{2)r rota¬ 
tions acting in the abvious way, and U(NrNl) gauge transformations acting 

jL,R 


B; 


U B br u 


-1 


(4.122) 


if the action (4.119) is considered as a matrix model, the radius drops out 
using (4.84). The equations of motion for Bf 1 are 

1 




N 2 


'} + ( B i + RijkB, B k ) 


3 3 4 J 1 ^ * 1 3 

+ ie ijk [Bf, (B r + ie krs B R B R )] + [. B f, [Bf , B R ]\ = 0 


(4.123) 
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and those for B R are obtained by exchanging L -H- R. By construction, 
the minimum or g round sta te of the action is given by F = <p = 0, hence 
B R ' R = X R ' R as in ( 4.8l|4.81 ) up to gauge transformations 
approach on CP 2 


for a similar 

We can therefore expand the covariant coordinates B R 


125 


and B R around the ground state 


P“ = A“ + RAl (4.124) 

where a G {L, R} and is very small, Then A R,R transforms under gauge 
transformations as 

a l, r X RR = U A rr U~ l + U[Xf’ R , U~\ (4.125) 

and the field strength takes a more familiar form (We do not distinguish 
between upper and lower indices L, R .) 

f,lsl = + 

Fmrn = + 

F, Li R = + 

So far, the spheres are described in terms of 3 Cartesian covariant coordinates 
each. In the commutative limit, we can separate the radial and tangential 
degrees of freedom. There are many ways to do R this; perhaps the most 
elegant for the present purpose is to note that the terms f ip L + cp R in the 
action imply that <Pl,r is bounded for configurations with finite action. Using 

‘Pl = ^Af + A r ^ + A r A r , (4.128) 

and similarly for <fR it follows that 

x i A; + A;x i = o(fy, (4.129) 

for finite A a { . This means that A ( - is tangential in the (commutative) large N 
limit. Alternatively, one could consider (f)i = Nip^ , which would acquire a 
mass of order N and decouple from the other fields. The commutative limit 


(4.126) 

(4.127) 
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4.3. GAUGE THEORY ON FUZZY S 2 x S 2 


of (4.119) therefore gives the standard action for electrodynamics on S 2 x S 2 


S 


- f 

2 g 2 Js 2 xs 2 


Tpt -pt 
r ia jb r ia jb"> 


(4.130) 


with a, 5 = L, A , at the north pole xb R = R , one can replace 


r A 


L,R 


R 


—y —6 


d 


lJ dx B ’ R 


(4.131) 


In the commutative limit, so that upon identifying the commutative gauge 
holds Af via 


Acl)L,R _ a l,r 


(4.132) 


the held strength is given by the standard expression 


l a ( d ) R aR a ( C 0 L 


F\l jR — 




(4.133) 


4.3.1 U(k) gauge theory 

The above action gener alizes immed iately to the nonabclian case, keeping 
precisely the same action ( 4.119|4.120 ) , but replacing the matrices B R,R by 
fc(YjjYi) x Ic(NrNl) matrices. The constraint term will then impose as 
ground state 


A f /R ® l fc xfc. 
Expanding the covariant coordinates 

B F R = \f /R ®l kxk + A R/R T a , 


(4.134) 


(4.135) 


in terms of the Heilman matrices T a , the action (4.119) is the fuzzy version 
of nonabelian U(k) Yang-Mills on S 2 x S 2 . 
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4.3. GAUGE THEORY ON FUZZY 5 2 x S 2 


4.3.2 A formulation based on SO(6) 

The above action can be cast into a nicer fo rm by assembling the matrices 

Since it is natural from 


126 


L R • • 

B i ’ into bigger collective matrices, following 
the fuzzy point of view to embed S 2 xS 2 C 1Z 6 with corresponding embedding 
of the symmetry group SO(3)l x SO(3)r C 50(6), we consider 


= (Sf. B 


R\ 


(4.136) 


to be the 6-dimensional irrep of so(6) = su(4). Since (4) © (4) = (6) © (10) , 
it is natural to introduce the intertwiners 


7m = (7 ?, 1 ?) = (7m)“ ,/? ( 4 - 137 ) 

We could then assemble our dynamical fields into a single 4 (NrNl) x 
4 (N r N l ) matrix 


B = B^ + const. 1 

(4.138) 

Of course the most general such 3{N R N L ) xA{N r Nr) matrix contains far too 
many degrees of freedom, and we have to constrain these B further. Since 

577(4) acts on B as 


B ->■ U t BU, 

(4.139) 

the 7 m can be chosen as totally anti-symmetric matrices, 
singles out the (6) C (4) © (4). One can moreover impose 

which precisely 

(it)' = it, 

(4.140) 

(it)' = -it, 

(4.141) 

and 


7 i 7 j = 3ij T itijk'Yk 

(4.142) 

lilf = ~ ie ijklk 

(4.143) 

!■,!?} = o 

(4.144) 


the representation satisfy 




4.3. GAUGE THEORY ON FUZZY S 2 x S 2 


With these we define the 4 -dimensional antisymmetric matrices 

71 = cr 1 (g) cr 2 , ^ = a 2 ® 1, 7 3 = a 3 ®a 2 , 

7 ^ = in 2 ® a 1 , 7 ^ = il (8) cr 2 , 7 ‘ i R = i(y 2 ®a 3 


(4.146) 


Where a 1 , a 2 and a 3 is the the Panli matrices. We can either separate them 
again by introducing two 4 (NrNl) x 4 (NrNl) matrices, 


B L = \+ B.V, B R =‘- + 


R ry R 
l 11 


(4.147) 


we can use the 7 ^ with the above properties to construct the 8 x 8 Gamma- 
matrices 

0 7 ^ 

7O 0 


= 


which generate the S , 0(6)-Clifford algebra 

YY^ + -7 yd 0 


{r^r"} = 


0 




= W lv 


This suggests to consider the single Hermitian 4 (NrNl) x 4 (NrNl) matrix 


C = T^B U + C 0 = 


0 B L 
B l 0 


+ 


0 B R 
-B r 0 


= :C L + C R (4.148) 


where C Q = Ck + G? denote the constant 8 x 8 -matrices 


<7 = 7777 = i / 0 1 


7 = -)rf7rf = 


1 0 

0 1 

-1 0 


in the above basis. Using the Clifford algebra and the above definitions one 
then finds 


c 2 = b„b ii + 1 -- 


(4.149) 


The field strength F /IU coincides with the definition in (4.121) if written in 
the L — R notation, 


Fia jb t \-^ia 1 -Rjfe] T SabtijkBka,’ 


(4.150) 
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Therefore the action 

N 2 N 2 — 1 

S = Tr((c 2 - — ) 2 ) = 8tr (B.B, -— ) 2 + 4tr F, V F, V (4.151) 


with 


(B^B^ 



(■BiLBiL — 



+ {BjjiBjji — 



) 2 ,(4.152) 


and 


N 2 = 


Nr + N 2 


R 


(4.153) 


4.3.3 Stability analysis of the SO(6) 


Consider the action (4.151) We will split off the radial degrees of freedom 
for large N by setting 


BiL — \l + A iL — A iL + AiL + XiL^L 


(4.154) 


and similarly for B iR . The stability of our geometry will depend on the 
behavior of <f> L and Q R . We calculate that 

N 2 - 1 

Bf.B,, --— = N($ l + $ r ) + $ l $ l + $ r $ r (4.155) 

+ A^Afj, — [A m , A^] + 0(1/N), 

where we used that \i a Ai a = 0 and therfore both ,Ai a Xi a = 0(1/N) and 
Aia[Ka, •] = 0(1/N) for a — L, R. Setting 

&L + &R = $i, (4.156) 

&L ~ &R = $2 


we get 


N 2 - 1 

BfjB^ --- = #$! + $!$! + $ 2^2 

+ A^Aft — [A m , . 4 .^] + 0(1/N), 


(4.157) 
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4.3.4 Breaking 50(6) —>■ 50(3) x 50(3) 

To obtain the original action ( |4 .119 ) for S ' 2 x 5 2 , we have to break the 
SO ( 6 )— symmetry down to 50(3) x 50(3). We ca n do th is by using the left 
and right gauge fields C L and C R introduced in (4.148) separately. Their 
squares are 


Of — B iL B iL + - + ^ ^ y 'j iL ijkBjLBkL ), 

Or = B iR B iR + - — f ^ (B iR + ieij k Bj R B k ). 

As both 7 ^, 7 ^ are traceless, we have 

N 2 N 2 

S break := 2Tr((0i-^)(0 2 -^)) (4.158) 

AT 2 — 1 N 2 — 1 

= 16Tr((5 iL R a --^)(5 iR R iR ^-L)). 

With these terms we can recover our action as 

N 2 N 2 N 2 

5 = 5 — Sbreak '■= 2Tr((U 2 —) — 2(0| )(C R -^)) (4.159) 

/V 2 _1 /v 2 _1 1 

= 8tr {(B iL B iL - - L —) 2 + (. B iR B iR - -^—) 2 + -i^F^.160) 
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4.4 Fuzzy Tori 

We consider a toroidal lattice with lattice spacing a and N sites in every 
dimensions. The lattice size is then L = Na. We consider the unitary 
operators 

Z t = exp (i—Xi) , = 1. (4.161) 

The second condition simply restricts the points to Xi G aZ. We have imme¬ 
diately the commutation relations 

[xi,Xj] = i9ij -v=> ZiZj = exp(— 2niQij)ZjZi , O = — 6. (4.162) 

1 j 

We consider the case 9{j = 6Qij in two and four dimensions where 


Q 


o 1 
-1 0 


Q 


0 1 0 0 \ 

-10 0 0 

0 0 0 1 

0 0 - 10 / 


(4.163) 


The momentum in each direction will be assumed to have the usual period¬ 
icity, viz 


27 rnrii 


(4.164) 


The period of rrii is exactly N. The range of m t is [0, N — 1] or equivalently 
[—(N — l)/2 ,+(N — l)/2] and hence we obtain in the large lattice limit 
L —» oo the cutoff 


A = —. 
a 


(4.165) 


The quantization of the noncommutativity parameters 9 and 0 are given by 


« = ^,e = £. 


(4.166) 


In other words, we have rational noncommutativity 0, for N > 2, and hence 
a finite dimensional representation of the algebra of the noncommutative 
torus exists. In general we require N to be odd for 0 to come out rational 
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and thus be guaranteed the existence of the fuzzy torus. The cutoff in this 
case becomes 


A = 


Ntt 

IT' 


(4.167) 


This is consistent with the result of the fuzzy CPU 

The full Heisenberg algebra of the noncommutative torus includes also 
the fuzzy derivative operators 


27iidj 


Dj = exp (adj) , 0 :/ Z,I)J = exp( — - )Zi. 


(4.168) 


J N 

In two dimensions a finite dimensional N x N representation is given in terms 


of the clock and shift operators (with 00 = 

exp(27r?'0)) 


Ti = 

f 0 1 

0 0 1 

\ 

, r 2 = 

/ 1 

UJ 

iO 2 

CO 3 

\ 


V 1 - - - 

0 1 

0 ) 


V 

• / 


(4.169) 


by 


Z\ —T‘2 i Z 2 — Ti , D\ — (Ti 


JV+1 


D 2 = (Tj 


N+l 


(4.170) 


The solution in four dimensions is obtained by taking tensor products of 
these. Thus a real scalar field $ on the fuzzy torus is a hermitean J\f x J\f 
matrix where J\f = N d / 2 , i.e. the space of functions on the fuzzy torus is 
Mat (A/", C). Furthermore, the integral is defined by the usual formula 


= {2-Kd) dl2 Tr. 


'fuzzy torus 


(4.171) 


A basis of Mat (A/*, C) is given by the plane waves on the fuzzy torus defined 
by 


N d / A 


IUr‘IWw« 


i= 1 


%<j 


ij TfliTflj 1 


N d / 4 


exp (ikiXi). (4.172) 
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They satisfy 

<!>%, = > Tr &m' = <W- 

A noncommutative $ 4 theory on the fuzzy torus is given by 


5 = (2Tr6) d/2 Tr 


1 £(a»a + - 4) 2 + " f * 2 + V 


2a 2 


We expand the scalar field <f> in the plane waves (j)^ as 

d* = ^ ^ 

fh 

We compute immediately 

D&D+ = Y^^rnD^Df 


Y, t’mtpm exp( 


2nimi 

N 


Hence 


Tr{D,^D] 


+U 


= 

rh 

= TVd> 2 . 


Thus the action can be rewritten as 

“ 1 


5 = (2 Na 2 ) d/2 Tr 


yy - Ad>A + $) + 


A 


(4.173) 


(4.174) 


(4.175) 


(4.176) 


(4.177) 


(4.178) 


We compute the kinetic term and the propagator given respectively by 

m“ |. (4.179) 


_ / 2 

- y d> rH T+ — Y(1 - cos ah) + 


< dv/ht, >= 


s k s k' 


tv - cos ah) + rn 2 ' 


(4.180) 
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Thus the behavior of the propagator for large momenta is different and as 
a consequence the calculation of on fuzzy tori will be different from the 
result obtained using a sharp cutoff. We get 127 
for d = 2 


< 


d 2 x&(x) >= V 


d 2 r 


' fuzzy torus 


(2vr) 2 JT(1 — cost,) + m 2 a 2 / 2 


, (4.181) 


for d = 4 


/‘ >4 i 2 / \ va 2 r dS- i 

J fuzzy torus ^ J 0 ( 2 tt ) 4 (1 - COS r t ) + m 2 d 2 /2 


.(4.182) 
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Chapter 5 

Emerent Geometry in 
Yang-Mills Matrix Models and 
Fuzzy Sphere 5 2 and 5 2 x S' 2 


5.1 The Model 


A commutative/noncommutative space in Connes’ approach to geometry 
is given in terms of a spectral triple (A, A, H) rather than in terms of a set of 
points 128 . A is the algebra of functions or bounded operators on the space, 


A is the Laplace operator or, in the case of spinors, the Dirac operator, and 
% is the Hilbert space on which the algebra of bounded operators and the 
differential operator A are represented. 

In the IKKT model the geometry is in a precise sense emergent. The 
algebra A is given, in the large N limit, by Hermitian matrices with smooth 
eigenvalue distribution and bounded square trace 129 . The Laplacian/Dirac 


operator is given in terms of the background solution while the Hilbert space 
% is given by the adjoint representation of the gauge group U(N ). 

In this article we will study IKKT Yang-Mills matrix models with quartic 
mass deformations in three and six dimensions with SO(3) and 50(3) x 
50(3) symmetries which will lead naturally to the fuzzy two-sphere S 2 N and 
to the fuzzy four-sphere S% x S 2 N respectively. 

Noncommutative gauge theory on the fuzzy two-sphere 130,131 was in¬ 
troduced in 


132 


moving m 


It was derived as the low energy dynamics of open strings 

This theory 


a background magnetic field with 5 3 metric in 1133 
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5.1. THE MODEL 


consists of the Yang-Mills term which can be obtained from the reduction to 
zero dimensions of ordinary U(N ) Yang-Mills theory in 3 dimensions and a 
Chern-Simons term due to Myers effect 
turbatively in 


134 


135 and 136 and nonperturbatively in 137 


The model was studied per- 
This model 


contains beside the usual two-dimensional gauge held a scalar fluctuation 


normal to the sphere. In 138 a generalized model was proposed and studied 


in which this normal scalar held was suppressed by giving it a potential with 
very large mass. This was studied further in 139,140 where the instabil¬ 


ity of the sphere was interpreted along the lines of an emergent geometry 
phenomena. 

In 


141 an elegant random matrix model with a single matrix was shown 


to be equivalent to a gauge theory on the fuzzy sphere with a very partic¬ 
ular form of the potential which in the large N limit leads to a decoupled 
normal scalar fluctuation. In 142-145 an alternative model of gauge theory 


on the fuzzy sphere was proposed in which held configurations live in the 
Grassmannian manifold U(2N) / (U(N + 1) x U(N — 1)). In 142,143 this 


model was shown to possess the same partition function as the commutative 
model via the application of the powerful localization techniques. 

Noncommutative gauge theory on the fuzzy four-sphere x S 2 N which 
is given by a six matrix model with global SO( 3) x SO (3) symmetry con¬ 
taining at most quartic powers of the matrices was proposed in [l46|. The 
value M — 1/2 of the mass deformation parameter corresponds to the model 
studied m which can also be shown to correspond to a random matrix 
model with two matrices. This theory involves two normal scalar helds plus 
a four-dimensional gauge held. Again the mass deformation parameter M 
is essentially the mass of these normal fluctuations and thus for large M 


these scalar helds become weakly coupled. In 148 an interpretation of these 


normal scalar helds as dark energy as dark energy is put forward. 

Before we give a summary of the main results reported in this article we 
will hrst summarize the main models studied here and their normalization. 
The most general mass-deformed IKKT Yang-Mills matrix model in three 


dimensions up to quartic power in the gauge held X a is given by 139 


S[X] = NTr 


2 ia 


- ~[X a , X b ] 2 + —— e abc X a X b X c + M(X 2 ) 2 + PX 2 
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5.1. THE MODEL 


The Steinacker model 


141] corresponds to M = 1/2 and 


P 


2 a 2 

IT 


(1-4c 2 M). 


The model with M — 0 will correspond to a perfect square action (exact 
Yang-Mills term ~ Fab)- We will study a one-parameter family (labeled 
by M ) with (3 fixed as above. We may also use the parameters m 2 and p 
defined by m 2 = 2 c 2 M and /3 = — ct 2 p especially in one-loop calculations. 
The coefficient c 2 is nothing else but the SU(2) second Casimir, viz c 2 = 
L 2 a = (N 2 - l)/4. 

The above model will lead to an emergent two-sphere S^, i.e. X a ~ L a . 
To obtain an emergent S^ 0 x S^- Q , with N = Nq, we consider two copies of the 
above model with gauge fields X a and Y a respectively, and then couple X a and 
Y a in the usual way, i.e. via commutators squared. The only modification 
required in the normalization explained above is to replace the quadratic 
Casimir as c 2 —» c° = (IVq — l)/4. Thus, we consider the mass-deformed 
IKKT Yang-Mills model in six dimensions given by 


S{x, Y] = S[X] + S[Y] - jTr[X a , n] 2 . 


In the above models a will play the role of the gauge coupling constant. 
More precisely, the inverse gauge coupling constant is found to be given by 
the combination 

5 2 = a 2 N ~ -. 

9 

In order to take the commutative large N limit we will need to study carefully 
the plots of various observables as functions of 5, and determine from their 
behavior under N, the scaled quantities which are stable in the strict N —y 
oo limit. These scaled observables will be termed collapsed and it may even 
happen that the inverse gauge coupling constant itself a 2 will be required to 
be collapsed, i.e. scaled with N, appropriately. 

Indeed, it is found in the case of that scaled observables are actually 
functions of a and not of a, where a the is the so-called collapsed coupling 
constant defined by 


a = a\/N = aN. 


This can also be confirmed by one-loop calculation. 
The main results of this article are as follows: 
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The dynamically emergent geometry, which is given by a fuzzy two- 
sphere S 2 N: in the 3—dimensional mass-deformed IKKT matrix models, 
is found to be stable for all values of the deformation parameter M. 
The c ritical gauge coupling constant a is found to scale as in equation 
(5.24), i.e. as a ~ 1/y/N. The sphere-to-matrix transition line is 


pushed to 0 and only one phase survives. 


• The 6—dimensional mass-deformed IKKT matrix model exhibits a phase 
transition from a geometrical phase at low temperature, given by a 
fuzzy four-sphere x background, to a Yang-Mills matrix phase 
with no background geometrical structure at high temperature. The in¬ 
verse temperature fd is here identified with the gauge coupling constant 
a. 


• The transition is exotic in the sense that we observe, for small values 
of M, a discontinuous jump in the entropy, characteristic of a 1st order 
transition, yet with divergent critical fluctuations and a divergent spe¬ 
cific heat with critical exponent a = 1/2. The critical gauge coupling 
constant is pushed downwards as the scalar field mass is increased. 


For small M, the system in the Yang-Mills phase is well approximated 
by 6 decoupled matrices with a joint eigenvalue distribution which is 
uniform inside a ball in R 6 . This gives what we call the d — 6 law given 
by equation (5.47). For large M, the transition from the four-sphere 
phase x Sj^ to the Yang-Mills matrix phase turns into a crossover 
and the eigenvalue distribution in the Yang-Mills matrix phase changes 
from the d = 6 law to a uniform distribution. 


• In the Yang-Mills matrix phase the specific heat is equal to 3/2 which 
coincides with the specific heat of 6 independent matrix models with 
quartic potential in the high temperature limit. Once the geometrical 
phase is well established the specific heat takes the value 5/2 with the 
gauge field contributing 3/2 and the two scalar fields contributing 1. 

This article is organized as follows. In section 2 we review the construction 
of noncommutative fuzzy gauge theory on the fuzzy sphere from random 
matrix theory, and write down our generalized three matrix model. In section 
3 we show by means of Monte Carlo that the emergent fuzzy sphere in this 
three matrix model is stable for all values of the mass deformation parameter 
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5.2. THE 3-DIMENSIONAL MASS DEFORMED YANG-MILLS 

MATRIX MODEL 


M. In section 4 we write down the analogous six matrix model, and present 
its one-loop quantization, and then discuss in great detail the calculation of 
the phase diagram by means of Monte Carlo. It is shown here that the fuzzy 
four-sphere S 2 N x is only stable for large values of the mass deformation 
parameter M. In section 5 we give a detailed discussion of the phases of the 
model using the eigenvalue distribution. We also revisit in this section the 
effective potential in order to prove the critical behavior of the theory. In 
section 6 various related topics are discussed briefly such as: emergent gauge 
theory in two dimensions, monopoles and instantons, topology change and 
stability of the fuzzy four-sphere x S%, critical behavior, Dirac operators, 
random matrix theory formulation, and generalization to fuzzy four sphere 
S 1 and fuzzy CP™. Section 7 contains our conclusion. 


5.2 The 3— dimensional mass deformed Yang- 
Mills matrix model 


A U(n ) gauge action on the fuzzy sphere can be derived from a simple 
1—matrix model as follows 141 . We introduce Pauli matrices r n and the 


three N x N matrices L a , which are the SU{ 2) generators in the irreducible 
representation of spin s = (N — l)/2, and define the matrix 


C — (- + T a L a )® l r 


It is a trivial exrecise to check that 


A 


C = (j(j + l)-(j) 2 ) 


(5.1) 


(5.2) 


where j is the eigenvalue of the operator J = L + cr/2 which takes the two 
values N/2 and (N — 2)/2. The eigenvalues of C are therefore N/2 with 
multiplicity n(N + 1) and —N/2 with multiplicity n(N — 1) . Hence 


C 2 = (f ) 21 


2 nN • 


(5.3) 


As it turns out this matrix C can be obtained as a classical configuration of 
the following 2nN —dimensional 1—matrix action 


^ = iy v Trtr2 


C 2 - — 


2 \ 2 
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Indeed, the equations of motion derived from this action reads 

N 2 

C{C 2 - —) = 0. (5.5) 

It is easy to see that C solves this equation of motion and that the value of 
the action in this configuration is identically zero ,i.e. S[C = C] = 0. 
Expanding around the vacuum C by writing 

C = — + Co + <7 a C a , (5-6) 

where Co and C a are nNxnN matrices and imposing the condition 

C 0 = 0, (5.7) 


we get 

1 1 

C 2 \ Co =0 — ^ + C 2 + -Cabc&cFab- ( 5 - 8 ) 

The curvature F ab is given in terms of C a by 

Fab = i[C a > Cb\ + £ a bcC c - (5.9) 


Hence, we obtain the action 




T i tv 2 — 1 

-TrFl + -Tr(Cl - -^-T- 


(5.10) 


The normal scalar held on the fuzzy sphere is given in terms of C a by 


$ = 


2^ 

The U (■ n) gauge action becomes 


Cl-L\ N 2 — 1 

, c 2 = L\ = ---. 


(5.11) 


S[C\ 


1 1 

7 2 n 


i 

4 


TrF 2 b + 2 m 2 Tr& 2 


The mass m 2 is given by the value of the Casimir, viz 

m 2 = C2- 


(5.12) 


(5.13) 
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We can bring this action into the form (with X a = 2aC a /3) 


S[X] = NTr 


\[X a , X b } 2 + ^e abc X a X b X c + M(X 2 ) 2 + j3X t 


„ r rri 
M = — 
2c 2 


2 a 2 


(3 = —a/j, = ——(1 — 2m). 
9 


We will also use extensively the parameter a defined by 

4 . 


= a A N z = ( - ) —. 


(5.14) 


(5.15) 


The parameter g is the gauge coupling constant. The classical absolute min¬ 
imum of the model is given by the fuzzy sphere configurations C a = L a . 
Expanding around this solution by writing C a = L a + A a yields a U(n) the¬ 
ory with a 3—component gauge field A where the extra normal component 
is given by <L. In the commutative limit N—>oc the field $ = n a A a is in¬ 
finitely heavy and hence it decouples. The curvature in terms of A a is given 
by F ab i [L a , Afe] A a ] Y ObcA c Y ^[A a , A^] 

The pure gauge action 5'[C'] becomes 


^F a A b iC b A a -\- e abc A c . 


S[C\ = 


1 


;TrF z b + 


2o 


2 rr<L 2 


1 

4 ? 


dn 9 N 2 f . 

4 X- + 2fJ^ W) 


Ag 2 N a0 g 2 N 

This is the action of a pure U(n) gauge theory on the ordinary sphere. The 


action (5.14) defines therefore a pure U(n ) gauge theory on the fuzzy sphere 


at least in the region of the phase space where the vacuum configuration 
C a = L a is stable. In this study we will deal mostly with n = 1 and hence 
Trl = N. 

The action (5.14) should be compared with the action (3.8) of 139 . Here 
m 2 = c 2 whereas m 2 is a free parameter in 139 . Furthermore, /? is fixed here 


as (3 = 2cr/9 for m 2 = 0 whereas /3 is a free parameter for m 2 = 0 in 139 


This action is also slightly different from the action (1) studied in 138 . Here, 
the m 2 = 0 theory is the perfect square action F 2 b /2 whereas the m 2 = 0 
limit of the action (1) studied in |138| is the Alekseev-Recknagel-Schomerus 
(ARS) stringy action given by F 2 b /2 + Chern-Simons or equivalently [l33 


S[X] = NTr 


1 r 2 ia 

- j[x 0 , x b f + — babc X a X b X c 


(5.17) 


The perfect square action F 2 J2 and the ARS action are connected by the 


critical line in the plane a — j3 of the model with M — 0 computed in 1491. 
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5.3 An emergent stable fuzzy sphere S 2 

5.3.1 Results 

The case of a single sphere contained in the three dimensional mass- 
deformed IKKT matrix model which is studied in this section is much sim¬ 
pler and differs from the case of S 2 x S 2 studied in the next section in two 
important aspects. 


1. There exists for every finite N only two distinct phases regardless of 
the value of the mass parameter M. The fuzzy sphere phase S 2 N and 
the Yang-Mills matrix phase. This is conhrmed by the eigenvalue dis¬ 
tribution. 


2. But since the critical value separating the two phases is found in one- 
loop and Monte Carlo to scale as 1/y/N the Yang-Mills phase is van¬ 
ishingly small in the limit N —> oo. 


Thus, effectively the three dimensional mass-deformed IKKT model (5.14) 
contains really a single phase. 

Let us also mention that the critical points are estimated in three different 
independent ways in both the three and the six mass-deformed IKKT matrix 
models. These are: 


1. From the intersection point of the expectation value of the actions for 
different values of N. 

2. From the peak or minimum of the specfic heat. 

3. From the discontinuity or maximum of the radius defined below. 
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5.3.2 Discussion 

We discuss now these issues in more detail. We want then to study, by 
means of the Monte Carlo method , the phase diagram in the plane a — M 


of the model (5.14) which we will rewrite again as 


S[X] = NTr 


\[X a , X b ] 2 + 2 -^e abc X a X b X c + M(X 2 ) 2 + f3Xl\ 


,18) 


The parameter M is taken in the range between M = 0 (perfect square 
action) and M — 1/2 (Steinacker’s action) and even beyond whereas (3 is 
given in terms of M by 


2 a 2 

(3 = —a 2 = —— (1 - 4c 2 M). 

y 


(5.19) 


The parameters f3 is therefore not independent and the two independent 
parameters of the model are a and M. 

The background minimal solution of this model is X a = ct(j)L a . In the 
limit m 2 —» oo, we have /i —» m 2 , and we hnd a critical line separating the 
fuzzy sphere phase solution with 0^0, from the Yang-Mills matrix phase 
solution with 0 = 0, given by the critical line derived originally from a one- 
loop analysis in 135| . Explicitly this is given in terms of the inverse gauge 
coupling constant a 4 rs-/ 1/g 2 by 


a 4 = 


81 

2m 2 


(5.20) 


We use the Metropolis algorithm for the Monte Carlo update. 

For each value of the mass parameter M, we have measured the critical 
point 5* at the minimum of the specific heat 


a, =< s 2 > - < s >' 


(5.21) 


with error bars estimated by the value of the step. See figure (5.2). The re¬ 


sults are shown on table (5.1). We observe from table (5.2) that the collapsed 


coupling constant is not a but it is given by 

a = a\f~N = aN. 


(5.22) 


Another measurement of the critical value a;* is given by intersection point 
of the actions with different values of N. See figure (5.4). The results are 


shown on table (5.3). 
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M/a 

N = 4 

N = 6 

N = 9 

N = 16 

N = 25 

0.5 

1.9 

1.5 

1.3 

1.0 

0.8 

1 

1.7 

1.45 

1.2 

0.9 

0.7 

2 

1.6 

1.4 

1.0 

0.8 

0.6 

5 

1.3 

1.1 

0.9 

0.65 

0.55 

10 

1.2 

1 

0.8 

0.6 

0.5 

50 

0.8 

0.7 

0.5 

0.4 

0.3 

100 

0.7 

0.5 

0.4 

0.3 

0.3 


Table 5.1: The critical values a* for different values of M and N. 


The two different measurements, which give an upper and lower estimates, 
of the critical point d* are included on the phase diagram (5.1). The data are 
fitted to straight lines and compared to theory in table (5.4). The theoretical 
prediction (5.20) in terms of the collapsed coupling constant d is given by 


3 

M 1 / 4 ' 


(5.23) 


This shows explicitly that the fuzzy sphere is absolutely stable in this model 
for all values of M, including the Steinacker value M = 1/2, since the inverse 
of the critical gauge coupling constant behaves as 


3 

ATAM 1 / 4 


0 . 


(5.24) 


The sphere-to-matrix transition line is pushed to 0 and only one phase sur¬ 
vives. This is the result advocated originally by Steinacker in 1141 . 

For completeness we also measure the radius of the sphere R as a function 
of the mass M and N. This is defined by the formula 


1 

R 


4> 2 a 2 C2 


TrX 2 a , c 2 



The results are shown on figure (5.5). 


2 

3 


(5.25) 
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M/a 

N = 4 

N = 6 

N = 9 

N — 16 

N = 25 

0.5 

3.8 

3.7 

3.9 

4.0 

4.0 

1 

3.4 

3.55 

3.6 

3.6 

3.5 

2 

3.2 

3.4 

3 

3.2 

3 

5 

2.6 

2.7 

2.7 

2.6 

2.75 

10 

2.4 

2.4 

2.4 

2.4 

2.5 

50 

1.6 

1.7 

1.5 

1.6 

1.5 

100 

1.4 

1.22 

1.2 

1.2 

1.5 


Table 5.2: The critical values a* for different values of M and N. 


M 

a 

0.5 

3.125 

1 

2.55 

2 

2.1 

5 

1.65 

10 

1.35 

50 

0.9 

100 

0.75 


Table 5.3: The critical values ct* from the intersection point of the actions. 


method 

a 

b 

Minimum C v 

-0.22 ±0.014 

1.29 ±0.036 

Theory 

-0.25 

0.48 

Intersection S 

-0.268 ±0.003 

0.936 ±0.007 


Table 5.4: The slop and the intercept of the straight lines used to fit the data 
with the theoretical prediction. 
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Phase diagram 



Phase diagram 



log(M) 


Figure 5.1: The phase diagram of the three dimensional Yang-Mills matrix 
model. 
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S2,M=2 



S2,N=9 



Figure 5.2: The specific heat in the three dimensional Yang-Mills matrix 
model. 


122 





5.3. AN EMERGENT STABLE FUZZY SPHERE S 2 


S2 



a 


Figure 5.3: C v 
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S2,M=2 



S2,M=2 



Figure 5.4: The action in the three dimensional Yang-Mills matrix model. 
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M=1 


cc 



Figure 5.5: The radius in the three dimensional Yang-Mills matrix model. 
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5.4 The emergent fuzzy S 2 x S 2 is stable only 
in the limit M —> oo 


The six dimensional mass-deformed IKKT model, which will contain S 2 x 
S 2 , is obtained by joining together two copies of the three dimensional mass- 
deformed IKKT model (5.14) for gauge fields X a and Y a respectively, then 

by adding a term 


i.e. 


N x S^ Q , with N = Nq, studied in this section is character- 


coupling the matrices X a and Y a in the usual way, 
proportional to the commutator squared. 

The case of S^ Q x S 2 
ized by three phases as opposed to the effectively one phase characterizing 
S 2 N . We have, for every value of the mass parameter M, the fuzzy S^ o x S^ o 
phase as expected for large values of 5, and a Yang-Mills matrix phase for 
small values of 5. In this case, in contrast with S^, the Yang-Mills phase is 
stable for every value of M as we take N —> oo, and the transition line goes 
as 5* ~ 1/M 0 - 25 , and as a consequence the fuzzy S^ 0 x S^ 0 is stable only in 
the limit M —> oo. 

However, the most distinct difference between the single fuzzy S 2 N case 
and the case of fuzzy S^ o x S 2 Nq is the appearance of a third phase between 
the above standard phases for large M. 

Indeed, there seems to exist another transition in M around the Steinacker’s 
value M = 0.5 where the profile of the eigenvalue distribution, for 5 = 0, 
changes from the d — 6 law (small values of M) to a uniform distribution 
(large value of M). Thus, the phase diagram develops a distinct third phase 
between the fuzzy S 2 Nq x Sj^ o phase and the Yang-Mills matrix phase for large 
values of M which looks like a cross-over. The transition line between the 
S 2 Nq x S|r (i phase and this crossover phase is determined at a saturated line 
around 5* ~ 4.2. Another possible interpretation of this phase is that of 
a strongly coupled gauge theory on the emergent background geometry and 
this is also supported by the measurement of the radius in this region. 
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5.4.1 The 6—dimensional mass deformed Yang-Mills 
matrix model 


Generalization of the above construction is straightforward to four di¬ 
mensions, i.e. to fuzzy S 2 x S 2 , and yields immediately the 6—matrix action 


S[X,Y] = NTr 
+ NTr 
+ NTr 


\ [X a , X b ] 2 + 2 -^-e abc X a X b X c + M(X 2 ) 2 + (3X 2 a 
\\Ya, Y b ] 2 + ^e abc Y a Y b Y c + M(Y 2 ) 2 + f3Y 2 


2 [Xa,Y b ] 2 


We choose A to be a perfect square such that 


(5.26) 


N = Nq. 


(5.27) 


The parameters of the model are given by 


TTl 2 

M = — , = - a 2 /j , ii = -(4c^M - 1) , c° 2 = 


Nq ~ 1 


(5.28) 


Before we report the Monte Carlo data we discuss the one-loop effective 
potential of this theory and its phase structure. 

The background solution (absolute minimum) of this model is by con¬ 
struction given by 


x a = a<j)L a ® 1 , Y a = act) 1 ® L a . 

(5.29) 

The fuzzy four-sphere is given explicitly by 



x\ + x\ + x\ = 1 , [x a ,x b ] = 

i 

/—^r^abc'^ci 

V C 2 

(5.30) 

~U 2/2 + 2/f = 4 7 [Vai Vb] — ■ 

i 

j—^^abcVci 

V C 2 

(5.31) 

[Xa,Vb] = 0, 


(5.32) 
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where x a = L a ® 1 / yci) and y a = 1 ® L a / y 

The L a are the SU (2) generators in the irreducible representation of spin 
s = (Nq — l)/2, and the background value of <fi is 2/3. By expanding around 
this solution, as X a = 2aC^ /3, Y a = 2 aCa^/3, Ca' 1 = L$ + Aa\ we obtain 
in the limit N —>■ oo a (7(1) gauge theory on S 2 x S 2 given by the action 


S{Cm,CV] = ~Tr 


-F 

4 ab 


(1)2 


+ 7 F ab 2 + \ F ab )2 + 2m 2 <h (1)2 + 2m 2 $ (2)2 (5.33) 


By comparing with equation (20) of 147 it seems to us that m 2 = c 2 as 
in the case of the single fuzzy sphere. The definition of the components of 
curvature tensor and the two scalar fields are obviously given by 


d?=i[Cf.cfl+t.xf , (5.34) 


$ 0 ) — 


c t 


0)2 


(5.35) 


5.4.2 Quantization at one-loop 

For simplicity let us go back to the case of a single sphere and discuss the 


derivation of the effective potential there first 135 


Quantization around this background, using the background field method 
gives, after gauge fixing in the Lorentz gauge, the effective action 


r[X a ] = A[Y a ] +-Tr log H-Tr log A 2 , 


(5.36) 


where the Laplacian operator is given explicitly by the formula 

tt ab = X'Xb - ‘Zfab + 4M(Y/ - c 2 o 2 )5 ab + SMX a X b + 2(0 + 2c 2 o 2 4£)%) 


The first term in (5.36) is due to the gauge field while the second term is 


due to the ghost field. In above the notation X a and T a \, means that the 
covariant derivative X a and the curvature F a b = i[X a , X^\ + ae a b c X c act 
by commutators, i.e X a (A) = [A" a ,A], T a b(A) = [F a b,A] where AEMatN- 
Similarly, X 2 (A) = [X a , [X a ,A]). 

The UV-IR mixing behavior on the fuzzy four-sphere in the model with 


M = /3 = 0 was studied in great detail in 146 
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it can be shown that the contributions from all the terms in the Laplacian 
H except the first are subleading 135,150 . Thus, the effective potential for 
X a = acj)L a is given by 


V 

2 Co 


= a 


- 1 -m - ii — 

4 3 4 2 


+ log cj ) 2 


(5.38) 


The calculation on the fuzzy four-sphere S 2 x S 2 proceeds in exactly the same 
way 146,1501. First we fix the Lorentz gauge, then compute the effective 
action, and then substitute the background matrices X a = a<j)L a 0 1, Y a = 
acj) 1 0 L a . The calculation of the classical part of the effective potential is 
trivial. The quantum part goes along the same lines as above. In particular, 
the contributions from all the terms in the Laplacian f! are negligible except 
the first one. Thus we get the logarithmic potential 


-Tr c jTr log cj) 2 — Trcf 2 = -N 2 log <^> 2 — N 2 log q i 2 . 


(5.39) 


On S 2 we have d — 3 whereas on S 2 x S 2 we have d = 6. We get then on 
S 2 x S 2 the effective potential 


V 

2N 2 


= 2 cV 


T'y +m T-' lJ 


+ log 4>- 


ry4 


— Ol n 


l2 1 


-h m -U- 

4 3 4 2 


+ log 4> 2 


(5.40) 


where we have redefined the coupling constant by 


N 2 

iV 0 „4 


a = a n 


(5.41) 


The difference between the result on S 2 and this result lies in the replacement 
a —> «o and the replacement C 2 —> c° in the definition of /i. The analysis 
of the phase structure is therefore identical (see section 5.2). 

The equation of motion reads then 


V 

‘IN 2 


— Q(n 


— cj) 2 + m 2 0 3 — fi(j) 


2 

+ -. 


(5.42) 


In the limit m 2 —> oo, we have // —y 4?n 2 /9, and we hnd a critical line 
separating the fuzzy sphere phase solution with (j) ^ 0, from the Yang-Mills 


129 
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matrix phase solution with 0 = 0, given by the formula 

, , 3 , / 32it. 

(1 + m )0* - -(l + y 1 + —) 


1 0 2 (0* + 2/t) 



t = fi{ 1 + m 2 ). (5.43) 

More detail on the derivation of this formula and the analysis of the phase 
diagram from the effective potential can be found in section 5.2. 

The small mass limit M —> 0 gives 

a ~ N (5.44) 


which diverges with N. This is precisely the divergent asymptotic behavior 
of the critical line of the model M — 0 observed in 1149 . 

The large mass limit M —> oo (or equivalently m 2 —> oo) of these 
equations is given in terms of a by 


a ~ 3 ( M ) / . (5.45) 

Thus as opposed to the case of a single sphere the collapsed coupling constant 
is a and not a which signals the persistence of the instability of the emergent 
four-sphere geometry as we will now show with the Monte Carlo results. 


5.4.3 Monte Carlo calculation 

We perform Monte Carlo simulation of the above six dimensional matrix 
model using the Metropolis algorithm. We vary M for different values of N. 
We work with N = 4-25 and M = 0.05 - 50. 


5.4.4 The action: 

• The intersection point of the average actions (entropies) for various 
values of N leads in this case to a robust measurement of the transition 
point between the fuzzy four-sphere S 2 x S 2 phase and the matrix 
Yang-Mills phase. See (5.6). This measurement in comparison with 
the theoretical prediction value (6.34) gives an under estimation of the 
critical point. The results are included on table (5.5). 
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For small values of M we also observe a discrete jump in the entropy 
at the transition point between the fuzzy four-sphere S 2 x S 2 phase 
and the matrix Yang-Mills phase. See the second graph of(5.6). The 
physics of this discontinuity is discussed in more detail in section 5 
using the approximation of the effective potential. 


The behavior of the action which couples the two spheres, viz 

S 12 =YM 12 = -^Tr[Y a ,n] 2 , 


(5.46) 


for small and large values of M, is shown on figure (5.7). We choose 


M = 0.1, and M = 120, for N = 16. We also plot Si and S' 2 , which 
are obviously defined, for comparison. 

We observe that for small values of M the action YM 12 is comparable 
to Si, whereas for large values of M it becomes quite negligible. This 
means in particular, that for large values of M, the approximation of 
the effective potential is expected to work well. 


5.4.5 The specific heat: 

• For small values of M, the geometric phase transition from the fuzzy 
four-sphere S 2 x S 2 phase to the matrix Yang-Mills phase is marked by, 


and is measured at, the peak of the specific heat. See figure (5.8). The 
peak becomes harder to resolve for larger values of M and N. This 
is in contrast with the case of the fuzzy two-sphere S 2 where a peak 
in the specific heat is observed only for very small values of M. The 


results are included on table (5.6). 


There seems to exist a new transition in M, around M ~ 0.5 which 
is the Steinacker’s value, beyond which the peak in C v ceases from 
marking the transition from the fuzzy four-sphere phase to the Yang- 
Mills matrix phase, and the specific heat develops a minimum where 
the transition actually occurs. See figure (5.9) and the results are 
included on table (5.7). Thus, for larger values of M above M ~ 0.5 
the geometric phase transition from the fuzzy four-sphere S 2 x S 2 phase 
to the matrix Yang-Mills phase is marked by, and is measured at, the 
minimum of C v . 
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• We also observe that the value at the peak of the specific heat saturates 
at around a ~ 4.1 for larger values of M. This peak corresponds in 
this case to a transition from a fuzzy four-sphere to a crossover phase 
as we will describe further shortly. 

5.4.6 The radius: 

We also measure the radii of the two spheres given in terms of TrY 2 and 
TrY a 2 respectively by the formula (5.25) with the substitution N —> Nq. 


From the Monte Carlo data we observe no difference, beyond and above 
statistical fluctuations, between the two radii. This is obvious from the figure 


(5.10). 


We observe for small values of M a discontinuity in the radius at the 


transition point as seen neatly on figure (5.10). We note that in this case it 


becomes harder to thermalize the system in the fuzzy four-sphere phase due 
to the zero modes of the matrices X a and Y a . 

For medium and larger values of M the discontinuity is smoothed out 
as shown on figure (5.12). The transition point in this case is taken, for 


medium values of M, at the maximum reached by the radius R in the fuzzy 
four-sphere phase before decreasing to zero in the Yang-mills matrix phase. 
For larger values of M, the transition point is taken at the point where the 
radius R drops below one. In summary, 


R 


a » a* fuzzy four-sphere phase 
a « 5* Yang-Mills matrix phase. 


This consists an independent measurement of the critical point between the 
fuzzy four-sphere phase and the Yang-Mills matrix phase. The data for 
different values of N is well collapsed and thus one can obtain from the 


radius a single estimate for the critical point shown on table (5.8). 


5.4.7 The phase diagram 

The critical line, and as a consequence the phase diagram, from the mea¬ 
surements of the action (intersection point or jump), the specific heat (peak 
and minimum) and the radius (jump, maximum and dropping below 1), to¬ 


gether with the theoretical calculation given by equation (6.34), are shown 


on figure (5.13). 
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In summary, we observe roughly three phases. The fuzzy four-sphere 
phase as expected for large values of 5, and a Yang-Mills matrix phase for 
small values of 5, and this is the case for every value of the mass parameter 
M. 

Therefore, in this case the fuzzy four-sphere phase is only stable in the 
limit M —> oo since the collapsed gauge coupling constant is 5 and not b, 
in contrast to what happens in the case of a single sphere. 

As we will discuss, in the next section, the Yang-Mills phase in this case is 
characterized by different eigenvalue distributions for large and small values 
of M. There seems to exist another transition in M around the Steinacker’s 
value M = 0.5 where the profile of the eigenvalue distribution, for 5 = 0, 
changes from the d = 6 law (small values of M ) to a uniform distribution 
(large value of M ). See below for a detailed discussion. 

Also, we observe that the phase diagram develops a distinct third phase 
between the fuzzy four-sphere phase and the Yang-Mills matrix phase for 
large values of M. This looks like a crossover phase. The critical line between 
the fuzzy four-sphere phase and this new phase, as measured by the peak 
of C v , saturates around 5 = 4.2. Another possible interpretation of this 
phase is that of a strongly coupled gauge theory on the emergent background 
geometry and this is supported by the measurement of the radius in this 
region. 
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5.4. THE EMERGENT FUZZY S' 


x S 2 IS STABLE ONLY IN THE 
LIMIT M —> oo 


M 

a 

M = 0.05 

8.50 

M = 0.1 

7.20 

M = 0.5 

3.75 

M = 1.0 

2.80 

M = 10 

1.55 

M = 20 

1.17 

M = 25 

1.11 

M = 30 

1.06 

M = 35 

1.02 

M = 40 

0.98 

M = 60 

0.89 

M = 80 

0.83 

M = 100 

0.78 


Table 5.5: The critical values cc* from the intersection of the action. 


M/a 

N = 4 

N = 9 

N = 16 

N = 25 

Extrapolation 

M = 0.05 

7.50 

7.60 

8.20 

8.70 

8.5614 ±0.345 

M = 0.1 

6.10 

6.70 

6.80 

7.00 

7.0262 ± 0.1505 

M = 0.5 

4.90 

5.00 

5.20 

5.30 

5.3040 ± 0.0813 

M = 1.0 

4.60 

4.70 

4.70 

4.80 

4.7914 ±0.0360 

M = 10 

4.20 

4.30 

4.20 

4.30 

4.2851 ± 0.0534 

M = 20 

4.20 

4.20 

4.30 

4.20 

4.2482 ± 0.0492 

M = 25 

4.30 

4.30 

4.30 

4.20 

4.2419 ± 0.0449 

M = 30 

4.30 

4.20 

4.30 

4.20 

4.2149 ± 0.0534 

M = 35 

4.30 

4.20 

4.20 

4.20 

4.1666 ±0.0168 

M = 40 

4.10 

4.30 

4.30 

4.20 

4.3087 ±0.0717 

M = 60 

4.20 

4.30 

4.20 

4.20 

4.2271 ±0.0531 

M = 80 

4.10 

4.10 

4.20 

4.20 

4.2063 ± 0.0374 

M = 100 

4.20 

4.30 

4.10 

4.20 

4.1788 ±0.0849 


Table 5.6: The critical values a* from the peak in C v . 
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x S 2 IS STABLE ONLY IN THE 
LIMIT M —> oo 


M/a 

N = 4 

N = 9 

N = 16 

N = 25 

Extrapolation 

M = 0.05 

7.50 

7.60 

8.20 

8.70 

8.5614 ±0.345 

M = 0.1 

6.10 

6.70 

6.80 

7.00 

7.0262 ± 0.1505 

M = 0.5 

3.70 

3.70 

3.85 

3.90 

3.8885 ± 0.0664 

M = 1.0 

4.60 

4.70 

4.70 

4.80 

4.7914 ±0.0360 

M = 10 

2.50 

2.30 

2.30 

2.35 

2.2622 ±0.05124 

M = 20 

2.20 

2.00 

1.90 

2.00 

1.8850 ±0.0571 

M = 25 

2.10 

1.90 

1.90 

2.00 

1.8912 ±0.0717 

M = 30 

2.10 

1.80 

1.70 

1.80 

1.5404 ± 0.0233 

M = 35 

2.00 

1.80 

1.65 

1.65 

1.4917 ±0.0516 

M = 40 

1.90 

1.70 

1.60 

1.55 

1.4979 ±0.0155 

M = 60 

1.70 

1.50 

1.55 

1.50 

1.4573 ± 0.0433 

M = 80 

1.60 

1.50 

1.50 

1.45 

1.4376 ±0.0168 

M = 100 

1.60 

1.40 

1.45 

1.40 

1.3573 ± 0.0823 


Tabic 5.7: The critical values a* from the minimum in C v . 


M 

a 

M = 0.05 

8.10 

M = 0.1 

6.70 

M = 0.5 

4.60 

M = 1.0 

4.20 

M = 10 

3.50 

M = 20 

2.50 

M = 25 

2.50 

M = 30 

2.45 

M = 35 

2.40 

M = 40 

2.40 


Table 5.8: The critical values ct* from the radius. 
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5.4. THE EMERGENT FUZZY S' 


x S 2 IS STABLE ONLY IN THE 
LIMIT M —> oo 


M=40 



N=25 



Figure 5.6: The intersection point of the action in the six dimensional Yang- 
Mills matrix model. 
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5.4. 


THE EMERGENT FUZZY S' 


x S 2 IS STABLE ONLY IN THE 
LIMIT M —> oo 



a 


N=16,M=120,Sj 2 

So 


* * * * * * % > 


N 1 " 

* 0.8 

A 0.6 - 

V 0.4 - 

0.2 - 
0 




111111111 l -H+Hf 1- 


-0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 


Figure 5.7: The action YM 12 
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5.4. THE EMERGENT FUZZY S 2 x S 2 IS STABLE ONLY IN THE 

LIMIT M —> oo 


M=1 



N=16 



a 


Figure 5.8: The specific heat (peak) in the six dimensional Yang-Mills matrix 
model. 
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5.4. THE EMERGENT FUZZY S 2 x S 2 IS STABLE ONLY IN THE 

LIMIT M —> oo 


N=16 



a 


Figure 5.9: The specific heat (minimum) in the six dimensional Yang-Mills 
matrix model. 
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5.4. THE EMERGENT FUZZY S 2 x S 2 IS STABLE ONLY IN THE 

LIMIT M —> oo 


M=0.05 




0 2 4 6 8 10 

a 


Figure 5.10: The radius in the six dimensional Yang-Mills matrix model for 
small values of M. 
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5.4. THE EMERGENT FUZZY S 2 x S 2 IS STABLE ONLY IN THE 

LIMIT M —> oo 


M=0.5 



a 


M=0.5 



a 


Figure 5.11: The radius in the two sphers. 
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N=16 
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Figure 5.12: The radius in the six dimensional Yang-Mills matrix model. 
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5.4. THE EMERGENT FUZZY S' 


x S 2 IS STABLE ONLY IN THE 
LIMIT M —> oo 
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Figure 5.13: The phase diagram of the six dimensional Yang-Mills matrix 
model. 
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5.5. EIGENVALUES DISTRIBUTIONS AND CRITICAL BEHAVIOR 


5.5 Eigenvalues distributions and critical be¬ 
havior 

5.5.1 Eigenvalues distributions from Monte Carlo 

The phase transition between the fuzzy four-sphere phase and the Yang- 
Mills phase can be characterized fully by the behavior of the eigenvalue dis¬ 
tribution across the transition line. This is by far the most detailed order 
parameter at our disposal. 


5.5.2 Small M: 


The behavior of the eigenvalue distributions for small values of M, such 
as M = 0.05 and M = 0.01, is depicted on figure (5.14). The two limiting 
behaviors are as follows: 

• For large values of the gauge coupling constant we observe a point 
spectrum given by the eigenvalues of the SU ( 2 ) generators in the largest 
irreducible representation which is of size N. 


Motivated by the work 151 it was conjectured in [1521 that the joint 
eigenvalues distribution of d matrices Xi, X 2 ,...X d with dynamics given 
by a reduced Yang-Mills action should be uniform inside a solid ball of 
some radius R. See also [?,?,?]. Let p(x i,..., x d ) be the joint eigenvalues 
distribution of the d matrices X±, X- 2 , ... and X ( i. We assume that 
p(x i, ...,£ 4 ) is uniform inside a four dimensional ball of radius r. The 
eigenvalues distribution of a single matrix, say X d , which is induced by 
integrating out the other d — 1 matrices is given by 

fb-i 


p(A) = 


V d (d - 1 ) 


(r — A 


2\(d-l)/2 


(5.47) 


For small values of the gauge coupling constant a —> 0 we observe a 
very good agreement with this law with d = 6 .We also include in the 


second graph of figure (5.14), the d = 4 law and the d — 3 parabolic 
law, for comparison. The fit for N = 16, M = 0.01, 5 = 0 gives a value 
of the radius r of the distribution given by 

r ~ 2. (5.48) 

We believe that this is a universal behavior at small 5 and small M. 
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5.5. EIGENVALUES DISTRIBUTIONS AND CRITICAL BEHAVIOR 


5.5.3 Large M: 


Some data is included on figures (5.15) and (5.16) 


For large values of the gauge coupling constant we still observe a point 
spectrum given by the eigenvalues of the 577(2) generators. 

For small values of the gauge coupling constant a —> 0 we observe ap¬ 
proximately a uniform distribution. This is neatly shown for M = 30, 
N = 9 and M = 1, N = 25 on figure (5.15). The Yang-Mills matrix 


phase is then characterized, for large values of M, by this one-cut uni¬ 
form distribution as opposed to the d = 6 law which characterizes the 
Yang-Mills phase for small values of M. The transition from the fuzzy 
four-sphere to the uniform distribution goes through a new phase or a 
crossover as we will now discuss. 

A new phase or a crossover: There seems to be another phase 
appearing for large values of M between the fuzzy four-sphere phase 
and the Yang-Mills phase. This is indicated by the transition from the 
distinct point spectrum in the fuzzy four-sphere phase to a phase where 
a strong gauge held is superimposed on the four-sphere background in 
such a way that the middle peaks flatten then disappears slowly in favor 
of the uniform distribution. The last peaks to go are the maximum 
and the minimum of the SU( 2) configuration. See figure (5.16) for 


N = 16,9 and M — 10 where this transition occurs at a = 4.2 at 
the peak of C v . Recall that the peak of the spechc heat saturates for 
large values of 5 at a = 4.2. However, this new phase may only be a 


crossover transition. Indeed, on the second graph of (5.17) we plot the 
behavior of X‘l as a function of a for N = 9 and M = 30. We observe 
that the profile of the eigenvalue distribution changes drastically only 
at the transition point to the uniform distribution. 

The coupling between the two spheres: The coupling between 
the two spheres can be probed by the eigenvalue distribution of the 
commutator i[X { , Y { ], A sample is shown on figure (5.19) for A = 9 and 
M = 10. We observe deep inside the fuzzy four-sphere phase that the 
commutator i[X\,Y\] is very different from the rotationally identical 
commutators i\X\, X 2 ] and i\Y \, Y 2 ], whereas the three commutators 
behave indistinguishably from each other deep inside the Yang-Mills 
matrix phase. In the middle phase, during the crossover, the eigenvalue 
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distribution of i[X\,Y\\ approaches quickly the prohle of the other two. 
The widths, for example, become less and less different, and shrink to 
a minimum value in the limit a —>■ 0. This is shown explicitly for the 
commutator i[X\, A 2 ] for N = 9 and M = 30. 


Steinacker’s value M = 0.5 and transition point in M: See (5.18) 


The behavior of the eigenvalue distribution deep inside the Yang-Mills 
matrix phase for M = 0.5, and other medium values of M, is neither 
given by the d = 6 law, nor it is given by a uniform distribution. This 
indicates that this value is near the triple point where the transition 
line between the fuzzy four-sphere S 2 x S 2 and the Yang-Mills matrix 
phase becomes a crossover phase. 

Rotational symmetry: The eigenvalue distributions of A 3 , Y :i . and 
the commutators i[X\, A 2 ], i[Y \, Y 2 ], and the squares A 2 and Y a 2 , in the 
fuzzy four-sphere phase are shown on the Erst graph of (5.17). This 


shows explicitly the rotational symmetry between the two spheres. 


The case of a single sphere: The physics in this case is very similar to 
the case of the fuzzy four-sphere and a sample of the eigenvalue distributions 


is included on figures (5.20) and (5.21). 
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Eigenvalue X 3 , N=9, M=0.05 



Eigenvalue M=0.01, a= 0 



A, 


Figure 5.14: The eigenvalue distributions for small values of M. 
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Eigenvalue X 3 , N=9, M=30 



X 6 N=25 , M=1 



Figure 5.15: The eigenvalue distributions for large values of M. 
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Eigenvalue X 3 , N=16, M=10 



Figure 5.16: The eigenvalue distributions for large values of M. 
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Eigenvalue, N=9, M=30 



X 


Eigenvalue Xj 2 , N=9, M=30 



Figure 5.17: Rotational invariance in the eigenvalue distributions for large 
values of M. 
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Eigenvalue X 3 , N=16, M=0.5 



Eigenvalue X 3 , N=16, M=0.5 



Figure 5.18: The eigenvalue distributions for the Steinacker value M — 1/2. 
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Eigenvalue i[X 1 ,X 2 ], N=9, M=30 



Figure 5.19: The eigenvalue distributions of the commutators i[Ah,Yi], etc. 
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S2,eigenvalues for N=9,M=0.5 



A. 


S2,eigenvalues for N=4,M=0.5 



A. 


Figure 5.20: The eigenvalue distributions on a single sphere. 
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o 11 g —i-- 1 — d m i 

-3-2-10123 


X 


Eigenvalue X 3 , N=9, M=2 



X 


Figure 5.21: The eigenvalue distributions on a single sphere (more). 
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X 6 N=25 , M=15 



A, 


X 6 N=9 , M=15 



Figure 5.22: The eigenvalue distribution for X e across the transition line. 
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X 6 N=9 , M=0.001 



Figure 5.23: The eigenvalue distribution for X (i across the transition line. 
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X 6 N=16,a= 0 



X 6 N=16,a= 0 



A, 


Figure 5.24: The eigenvalue distributions for X 6 for small and large values 
of M and <5 = 0. ^57 






5.6. THE EFFECTIVE POTENTIAL REVISITED 


5.6 The effective potential revisited 


We recall the formula for the effective potential 


V 

2iW 




+ log 0 2 . 


The classical solutions are given by the condition 


V' 
2 N 2 


0 . 


(5.49) 


(5.50) 


Clearly the effective potential is not bounded at 0 = 0. However, this should 
not pose any problem since this potential is valid for large values of a where 
we know that the fuzzy sphere exists. For smaller values of a the fuzzy 
sphere configuration ceases to exist and we enter the matrix phase as we have 
seen in Monte Carlo data. The classical potential admits as solutions 0 = 
0 (local minimum) together with the fuzzy four-sphere solution 0 + (global 
minimum), and a maximum of the barrier between them denoted by 0 _ 
(local maximum). The solutions 0± exist for t > — 1/4. As we decrease 
H towards negative values the global minimum 0 + becomes degenerate with 
0 = 0 at t — ^(l + m 2 ) = —2/9 and the height at the maximum of the barrier 
becomes 5 4 /324(l + m 2 ) 3 . There exists therefore a first order transition in 
/x, at fi = —2/9(1 + m 2 ) for every fixed m 2 , which is of the same character as 
the one at m 2 = 0, from a fuzzy-four sphere phase for fi > —2/9(1 + m 2 ) to 
X a = Y a = 0 for ii < —2/9(1+m 2 ) . For our value /i = 2(2m 2 —1)/9 the global 
minimum is always 0 = 2/3 and it is separated from 0 = 0 by a barrier. The 
effective potential admits four real solutions. The largest positive solution 
is 0 + while the other positive solution gives the local maximum and it will 
determine the height of the barrier in the effective potential. At the critical 
point these two solutions merge and the barrier disappears. This is different 
from the classical solution where the barrier never disappears. The solution 
0_j_ ceases to exist at the critical value determined by the condition 


V" 
2 N 2 


0 . 


(5.51) 
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We recall also that the critical value is found by solving (5.50) and (5.51) 
and is given implicitly by 


(1 + m 2 ) 0 * — - (1 + ^1 + —) 
1 0 *( 0 * + 2/0 


a: 


o* 


8 


t = /t( 1 + m 2 ). 


(5.52) 


The critical value is sent to infinity for 0 * = — 2/i which is equivalent to 
t = —1/4. Thus quantum mechanically the fuzzy four-sphere may exist for 
t > —l/A which is to be compared with the classical prediction t > —2/9. 
In this region the classical potential is always positive and thus one should 
consider all 577 ( 2 ) representations which are degenerate with X a = Y a = 0, 
i.e. for which / N —t 0 in the limit N —» 0 with Y2i n i = N. 

However, for large a the ground state is dominated by the representation with 
the smallest Casimir. The fuzzy four-sphere is therefore not stable in this 
regime and the critical line between the fuzzy four-sphere phase and the Yang- 
Mills matrix phase asymptotes the line t = —2/9 as shown in 149 . What 


interests us the most in this section is the expansion of the solution 0 around 
the critical value 0*. This is given in 139 . The only difference between the 


result on S 2 given in 139 and our result her lies in the replacement 


cl —y ol o 


and the replacement 


c 2 —> c 2 

in the definition of ji. We get the solution 


— 0* H-k—/ „ , . ~ \/ a o ~ a o* + ••• 


~2 a/ 30* + 4/i 


(5.53) 


(5.54) 


(5.55) 


This takes the form 


4 ( 21 V ) 1 / 2 1 - 

— 0* + (J , (7 — -g-7====== V Oi — Ol* + ... 


a 


2 t /30* + 4/u 


(5.56) 


159 
















5.7. CRITICAL BEHAVIOR FROM ONE-LOOP EFFECTIVE 

POTENTIAL 


In the large N limit we have 




2 NM 
9 


(5.57) 



(5.58) 


CK, 


8N 2 M 

R 2 


162 


Thus 


(5.59) 


— + cr , <7 — 


a 


a/M 


ex — ex* + 


(5.60) 


5.7 Critical behavior from one-loop effective 
potential 

Let us start by noting the Schwinger-Dyson identity 

4 < YM 1 > +4 < YM 2 > +4 < YM 12 > +3 < CSi > +3 < CS 2 > 

+ 4 < Quar, > +4 < Quar 2 > +2 < Quad 1 > +2 < Quad 2 >= 6)Y 2 61) 

We have the obvious definitions for the various observables 


YMj = NTr 


-[X a ,X b ] 2 | , YM 2 = NTr 


\[Y a ,Y b ] 2 | , YM 12 = NTr 


i [A/,W] 2 


(5.62) 


NTr 

q £abcX a XbX c 

, CS 2 = NTr 

2ia 

^ €abcY a Yi)Y c 


O 


O 


Qua^ = NMTr(X 2 )' 2 , Quar 2 = NMTr{Y 2 ) 2 . 
Quad! = NfiTrX 2 , Quad 2 = N(3TrY 2 . 


(5.63) 

(5.64) 

(5.65) 
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5.7. CRITICAL BEHAVIOR FROM ONE-LOOP EFFECTIVE 

POTENTIAL 


We obtain immediately a formula for the average action given by 

A = V jA < CS ' > +iA < 082 > +2A *= Quad > > +2 Quad2 

We compute 


4^ <CS ' >= i^ <CSa> = 


(5.67) 


1 ~ , 1 ~ , +cd 4 0 2 

— 2 < Quad, >= Wi < Q«ad 2 >= 

We also note here the formula for the related radius (with 0 O = 2/3) 


(5.68) 


1 

— =< 


-JrXl >=< 


R 02 5 2 C O a ' - ^2 5 2 C 0 a - ^ 


As>= (!)V. 


(5.69) 


Thus 


S 3 a 4 1 3 , 2 , 

N 2 ~ 2 AN 3^ + ^ ^ 

V _ ± 

/v 2 <+' 


( 5 . 70 ) 


3 4 / Q; , 4 a 4 fi 2 

2 _ (odT _ LL/W * 


3 4 0* T 3/i 
2 3 0 * + 2 fi 


- 2 - ( 5 . 71 ) 

This is the value of the average action or entropy in the fuzzy four-sphere 
phase at the transition point exactly for M large. For M small we get instead 
the (constant) value 


(5.72) 



5.7. CRITICAL BEHAVIOR FROM ONE-LOOP EFFECTIVE 

POTENTIAL 


In the Yang-Mills matrix phase we go back to the first line of equation (5.70) 
and set (j) = 0 to get the value 3/2 at a = 0. The transition to the Yang-Mills 
phase occurs quite suddenly for small M. Thus the entropy for small M has 
a discrete jump given by 


AS 

w 


4 

3' 


The average action can also be derived using the formula 


< S > 
N 2 


d 


— -(- CL , 

2 d5“'JV 2 


( 4 ). 


where (with X a = aD ai a 4 S = S ) 

Z = exp(— F) = j dD a exp (—— ln5 4 — a 4 S). 
The free energy is given by 

2^ = 2^ + ^ lllSl + COnStant ' 

Then we can compute the specific heat by the formula 

< S 2 > — < S > 2 


N 2 


<S> 

N 2 


d ,<S> 


a 


( M2 ) ’ 


dct 4 v N 2 


We compute immediately 


C v — — + 


m**+**% 


(5.73) 


(5.74) 


(5.75) 


(5.76) 


(5.77) 


(5.78) 


From this equation we can derive immediately the divergent part of the 
specific heat to be given by 


c v = cf + 


0*(0* + 2/x) aj 2 


4(2Y)V 2 ^30* + 4 n - a* 

3 1 / 2 1 


° v + 2V8MV8 


(5.79) 
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5.8. RELATED TOPICS 


The critical exponent of the specific heat is 1/2 which is precisely the value 
obtained in two dimensions. Similarly, the critical value 5* and the coefficient 
of the singularity become vanishingly small when M —» oo. The behavior 
of the background specific heat C x f deep inside the fuzzy four-sphere phase 
is computed as follows. The solution 0 of the equation of motion for large 
values of a is found to be given by 


27 


2 Ma 4 


+ ... 


(5.80) 


By substitution in (5.78) then using 0 ~ 2/3 since we are only interested in 


the background value of the specfic heat, and also using (5.57), we get 


a B = 


l + m^ +2 ^ 

1 + ' 

5 

2 ' 


(5.81) 


5.8 Related topics 

5.8.1 Emergent gauge theory in two dimensions: 

We have established that the fuzzy sphere is completely stable in the three 
matrix model (5.18) for any value of M. This includes Steinacker’s value 


M = 1/2. We can now speak in a consistent way about constructing U{ri) 
gauge theory on fuzzy S 2 . This entails the construction of monopoles sectors 
and the direct evaluation of the partition function on the fuzzy sphere, or 
alternatively its evaluation by means of localization technique, as a sum over 

141] 


instantons contributions. This has been done for Steinacker’s value in 


and 143 respectively. The result in this case was found to be identical to the 


result on the ordinary sphere. This shows that the fuzzy sphere is acting here 
as a regularized version of the sphere, and also shows that this matrix method 
is potentially a powerful new method for gauge theory. The construction 
of fuzzy monopoles and instantons and the corresponding Ginsparg-Wilson 
fermions on this stable fuzzy sphere can also be carried out along the lines 
of (153 155 . 
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5.8.2 A stable four-sphere S 2 x S 2 and topology change: 

The main conclusion we have reached in this article is the difference be¬ 
tween fuzzy S 2 n and fuzzy x S^. The fuzzy four-sphere is only stable in 
the large M limit. A simple modification of the six matrix model may lead 
to a stable four-dimensional geometry for any M. This consists in adding 
the terms 


NP^TrXl + Tri; 2 ), 


(5.82) 


with a particular coefficient f}\. This will modify the critical line (6.34) in 
such a way that a is replaced by a. Topology change can also be easily 
obtained in the above six matrix model by setting zero the mass parameters 
M and /3 on one of the spheres. This way we will have the possibility of 
a transition between the fuzzy four-sphere S 2 N x phase and the fuzzy 
sphere phase. This is what we call a topology change. This scenario was 
previously observed on fuzzy CP' 


156 


Critical behavior: The critical behavior of the above stable fuzzy two- 
sphere Sjy and also the critical behavior of the fuzzy four-sphere x 


02 


can also be determined more carefully along the line of 157 


m 


Comparison with (147 , the 2—matrix model and instanton calculus: 


The six matrix model (5.26) with M = 1/2 is precisely the action considered 
Jl47|. Since this model is only stable for large M, the corresponding 


instanton calculus should only be expected to be valid deep inside the fuzzy 
four-sphere phase. A 2—matrix model associated with the 6—dimensional 


Yang-Mills matrix model (5.26) can also be constructed starting from an 


SO(6) formulation. This should be contrasted with the SO (3) formulation 
on the fuzzy sphere which leads to (5.4). In the case of the fuzzy four-sphere 
S 2 n x S% the matrices C and D corresponding to X a and Y a respectively 


are found to be highly constrained 147 as opposed to the single constraint 
Cq = 0 found on the fuzzy sphere. The six matrix model considered in this 


article is also closely related to the action studied in 158 . 


147 


Dirac operator: There are two seemingly different formulations of the 
Dirac operator on the fuzzy four-sphere S 2 N x S 2 N . The one presented in 
is again based on SO(6) whereas the one presented in 
SO(3) x SO( 3). 


159 is based on 
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5.8.3 Fuzzy four-sphere S 4 and fuzzy CP 2 : 

The above analysis is expected to hold without much change on fuzzy 


160 


162 


are given by 156 
1 


and fuzzy CP 2 |l63 jl66|. The analogous actions on fuzzy CP' 


S = 


g 2 N 


Tr 


— ~^[D a ,Db] + ie a bcD a DbD c 


3n M 2 o M 2 (2n + 3) 2 ^ 

-|- Tr$ H - ° Tr<L 2 H ---— —Tr$. 

4g 2 N N N 36n 

(5.83) 


5.9 Conclusion 

In this article we have studied IKKT Yang-Mills matrix models with mass 
deformations in three and six dimensions. The 3—dimensional IKKT matrix 
models considered here are very similar to the ones studied in |135 138 


However, the dynamically emergent geometry, which is given by a fuzzy two- 
sphere S^r, is found to be stable for all values of the deformation parameter 


M. This was anticipated previously for M = 1/2 in 141 . Indeed, the crit 


ical gauge coupling constant a is found to scale as in equation (5.24), i.e. 


as a ~ 1 /\f~N. The sphere-to-matrix transition line is pushed to 0 and only 
one phase survives. In this case the fuzzy sphere acts then as a regulator 
of the commutative sphere, and as a consequence, fuzzy field theory and 
fuzzy physics, based on this emergent fuzzy sphere, makes full sense for all 
values of the gauge coupling constant. We have also studied in this article 
6 —dimensional IKKT matrix models, with global SO{3) x SO( 3) symmetry, 
containing at most quartic powers of the matrices proposed in 11461. The 
value M = 1/2 of the deformation corresponds to the model of [4471. This 
theory exhibits a phase transition from a geometrical phase at low temper¬ 
ature, given by a fuzzy four-sphere S 2 N x background, to a Yang-Mills 
matrix phase with no background geometrical structure at high tempera¬ 
ture. The geometry as well as an Abelian gauge field and two scalar fields 
are determined dynamically as the temperature is decreases and the fuzzy 
four-sphere condenses. The transition is exotic in the sense that we observe, 
for small values of M, a discontinuous jump in the entropy, characteristic of 
a 1st order transition, yet with divergent critical fluctuations and a diver¬ 
gent specific heat with critical exponent a = 1/2. The critical temperature 
is pushed upwards as the scalar field mass is increased. For small M, the 
system in the Yang-Mills phase is well approximated by 6 decoupled matrices 
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with a joint eigenvalue distribution which is uniform inside a ball in R 6 . This 
gives what we call the d = 6 law given by equation (5.47). For large M, the 
transition from the four-sphere phase to the Yang-Mills matrix phase turns 
into a crossover and the eigenvalue distribution in the Yang-Mills matrix 
phase changes from the d = 6 law to a uniform distribution. 

In the Yang-Mills matrix phase the specific heat is equal to 3/2 which 
coincides with the specific heat of 6 independent matrix models with quartic 
potential in the high temperature limit and is therefore consistent with this 
interpretation. Once the geometrical phase is well established the specific 
heat takes the value 5/2 with the gauge field contributing 3/2 and the two 
scalar fields contributing 1. This should be contrasted with the case of 
in which the specific heat in the Yang-Mills matrix phase is equal to 3/4, 
coinciding with the specific heat of 3 independent matrix models with quartic 
potential in the high temperature limit, while in the geometrical phase the 
specific heat takes the value 1 divided equally, i.e. with the gauge held 
contributing 1/2 and the two scalar fields contributing 1/2. The counting on 
S 2 n is clear cut. In the sphere phase which coincides with the perturbative 
region of the theory the gauge field contributes 1/2 [167 and the scalar field 
contribute 1/2 (trivial to check in the quartic matrix model for very large 
positive values of the mass parameter 168,169]). The counting on x 
is more involved. There are here two scalar fields and one gauge field with 
4 components. Again, it is very natural to suppose that each scalar field 
will contribute 1/2 (since they are free and dimension does not enter in the 
quartic matrix model). Therefore, the gauge field will contribute 3/2. This 
is the picture we also get, at least formally, by assuming that the gauge held 
is free and Abelian (which is true deep in the sphere phase in the large N 
limit) and then fixing the gauge in the axial gauge. The 6—dimensional IKKT 
Yang-Mills matrix models studied here present thus an appealing picture of a 
4—dimensional geometrical phase emerging as the system cools and suggests 
a scenario for the emergence of geometry in the early universe. See 148 and 
references therein. 
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Chapter 6 

Emergent Geometry in the 
multitrace quartic matrix 
models 


6.1 Introduction and Motivation 


The original motivation for this work is the theory of noncommutative <f > 4 
which we now briefly describe. A scalar phi-four theory on a non-degenerate 
noncommutative Euclidean spacetime is a a three-parameter matrix model 
of the generic form 

S = Tr H (aMAM+ bM 2 + cM 4 ). (6.1) 


The Laplacian A captures precisely the underlying geometry, i.e. the met¬ 


ric, of the noncommutative Euclidean spacetime in the sense of 170,171 


This theory can be regularized non-perturbatively using N x N matrices in 
an almost obvious way, i.e. the Hilbert space H can be taken to be finite 
dimensional of size N. This theory exhibits the following three known phases: 


• The usual 2nd order Ising phase transition between disordered < M >= 
0 and uniform ordered < M >~ ljv phases. This appears for small val¬ 
ues of c. This is the only transition observed in commutative phi-four, 
and thus it can be accessed in a small noncommutativity parameter 
expansion. 


• A matrix transition between disordered < M >= 0 and non-uniform 
ordered < M >~ 7 phases with y 2 = 1^. This transition coincides, for 
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very large values of c, with the 3rd order transition of the real quartic 
matrix model, i.e. the model with a = 0, which occurs at b = — 2 a/ Nc. 
In terms of b — bN ~ 3 / 2 and c = cN~ 2 this reads 


b = — 2 a/S. 


( 6 . 2 ) 


This is therefore a transition from a one-cut (disc) phase to a two-cut 
(annulus) phase 168 169. See also 172 |l73| . 


A transition between uniform ordered < M >~ 1^ and non-uniform 
ordered < M >~ 7 phases. The non-uniform phase, in which trans¬ 
lational/rotational invariance is spontaneously broken, is absent in the 
commutative theory. The non-uniform phase is essentially the stripe 


phase observed originally on Moyal-Weyl spaces in 174,175 . 


Thus, the uniform ordered phase < $ >~ ljv is stable in the theory (6.1) 


This fact is in contrast with the case of the real quartic matrix model V = 
Tr u^bM 2 + cM 4 ) in which this solution becomes unstable for all values of 
the couplings. The source of this stability is obviously the addition of the 
kinetic term to the action. 


The non-uniform ordered phase 176 is a full blown nonperturbative man¬ 


ifestation of the perturbative UV-IR mixing effect 177 which is due to the 


underlying highly non-local matrix degrees of freedom of the noncommuta- 
tive scalar held. 

The above picture of the phase diagram holds for noncommutative phi- 
four in any dimension, and the three phases are all stable, and are expected to 
meet at a triple point. The phase structure in four dimensions was discussed 


using the Hartree-Fock approximation in 174 and studied by means of the 


Monte Carlo method, employing the fuzzy torus 178 as regulator, in 175 


In two dimensions the noncommutative phi-four theory is renormaliz- 


able 179 . The regularized theory on the fuzzy sphere 130,131 is given by 


the action ( 6 . 1| ) with a finite dimensional Hilbert space H of size N and a 
Laplacian A = [L a , [L a ,..]] where L a are the generators of £77(2) in the IRR 
of spin (N — l)/ 2 . 

The above phase structure was confirmed in two dimensions by means 
of Monte Carlo simulations on the fuzzy sphere in [180,181 . Indeed, fuzzy 
scalar phi-four theory enjoys three stable phases: i) disordered (symmetric, 
one-cut, disk) phase, ii) uniform ordered (Ising, broken, asymmetric one-cut) 
phase and iii) non-uniform ordered (matrix, stripe, two-cut, annulus) phase. 
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The phase diagram is shown on the two graphs of figure (6.1) which were 
generated using the Metropolis algorithm. 

The problem of the phase structure of fuzzy phi-four was also studied 

I82}p6 


by means of the Monte Carlo method in 


The analytic derivation 


of the phase diagram of noncommutative phi-four on the fuzzy sphere was 
attempted in 


127,187-193 


The related problem of Monte Carlo simulation of noncommutative phi 
four on the fuzzy torus, and the fuzzy disc was considered in 


196 


175, 194 


and [195] re spectively. For a recent study see 

In (1861 the phase diagram of fuzzy phi-four theory was computed by 
Monte Carlo sampling of the eigenvalues X t of the scalar field M. This was 
possible by coupling the scalar field M to a 17(1) gauge held X a on the fuzzy 
sphere which then allowed us, by employing the U(N ) gauge symmetry, to 
reduce scalar phi-four theory to only its eigenvalues. The pure gauge term is 
such that the gauge held X a is fluctuating around X a = L a . 

Another powerful method which allows us to reduce noncommutative 
scalar phi-four theory to only its eigenvalues, without the additional dynam¬ 
ical gauge held, is the multitrace approach. The multitrace approach was 
initiated in 187,188 . See also 193 for a review and an extension of this 


method to the noncommutative Moyal-Weyl plane. For an earlier approach 

The 


see 


127 and for a similar more non-perturbative approach see 189 -192 


multitrace expansion is the analogue of the Hopping parameter expansion on 
the lattice in the sense that we perform a small kinetic term expansion, i.e. 


expanding in the parameter a of ( 6 . 1 ), while treating the potential exactly. 


This should be contrasted with the small interaction expansion of the usual 
perturbation theory. The effective action obtained in this approach is a 
matrix model which can be expressed solely in terms of the eigenvalues A* 
and which, on general grounds, can only be a function of the combinations 
T 2n oc JA /- (A* — A j) 2n . To the lowest non-trivial order we get an effective 
action of the form [188 ,189,193 




* *7 U 

“^ 7 ^ 2,1 2^ (A; ~ \) 7- Y2T Va ' 1 / ^ ~ 

*7U *7U 


-Vi 


,2[^(Aj - Xj) 2 ] + .. 
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The logarithmic potential arises from the Vandermonde determinant, i.e. 
from diagonalization. The coefficients V 2 ,i, Uyi and v-i;i are given by t’ 2 .i = 
+1 , V 4.1 = 0 , V 2,2 = 1/8. Furthermore, it is not difficult to convince 
ourselves that the above action is a multitrace matrix model since it can be 
expressed in terms of various moments m n = TrM n of the matrix M. 

The original multitrace matrix model written down 187 comes with dif¬ 
ferent values of r;’s and therefore, in the commutative limit N —» 00 , it 
corresponds to a phi-four theory on the sphere modulo multi-integral terms. 

Since these multitrace matrix models depend only on N independent 
eigenvalues their Monte Carlo sampling by means of the Metropolis algorithm 
does not suffer from any ergodic problem. The phase diagrams of these 
models obtained in Monte Carlo simulations will be reported elsewhere. 

The remainder of this article is organized as follows: 


1. Section 2: We describe our proposal for how fuzzy geometry can emerge 
in generic multitrace matrix models. 

2. Section 3: We apply our proposal to an explicit example. We will 
show that if the multitrace matrix model under consideration does not 
sustain the uniform ordered phase then there is no emergent geometry. 
On the other hand, if the uniform ordered phase is sustained then there 
is an underlying or emergent geometry. In particular, we will show how 

• i) to determine the dimension from the critical exponents of the 
uniform-to-disordered (Ising) phase transition, and how 

• ii) to determine the metric (Laplacian, propagator) from the Wigner 
semicircle law behavior of the eigenvalues distribution of the ma¬ 
trix M. 

3. Section 4: We conclude by giving a straightforward generalization to 
fuzzy CP" and fuzzy T". 
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6.2 The Proposal 

We start with a general multitrace matrix model rewritten in terms of 
the moments TrM n with generic parameters B, C, D , B , C , D , A ,...as 


1/ = BTrM 2 + CTrM 4 + D 


TrM 2 


+ R' (TrM) 2 + C'TrMTrM 3 + T>' (TrM) 4 + A TrM 2 (TrM) 2 + (6.4) 

This action includes the noncommutative phi-four model on the fuzzy sphere 
(6.3) and the multitrace matrix model of 1187] as special cases. It also in¬ 


cludes as special cases the multitrace matrix models obtained by expanding 
the kinetic term on i) fuzzy CP 


188,197 


and without the harmonic oscillator term 193 


The phase diagram of the action (6.4) will generically contain the ma¬ 


on ii) Moyal-Weyl spaces with 
and on iii) fuzzy tori 178 . 


trix one-cut-to-two-cut transition line separating the two stables phases of 
disorder and non-uniform-order. However, the uniform ordered phase will 
typically be unstable as in the case of the real quartic matrix model 


V = BTrM 2 + CTrM 4 


(6.5) 


Our proposal goes as follows. We can check for a possible emergence of 


geometry in the multitrace matrix model (6.4) by following the three steps: 


1. We compute the phase diagram of the model (6.4). If the uniform 


ordered phase remains unstable as in the case of the real quartic matrix 


model (6.5) then there is no geometry and the model is just a trivial 
deformation of (6.5). In the opposite case we claim that there is an 


underlying, i.e. emergent, geometry with a well defined dimension (step 
2) and a well defined Laplacian/metric (step 3). This means that we can 
rewrite the multitrace matrix model, in the region of the phase diagram 
where the uniform ordered phase exists, in terms of a scalar function 
and a star product with a noncommutativity parameter 6 by finding the 
appropriate Weyl map. As a consequence, a small noncommutativity 
parameter expansion can be performed and the the limit 6 —> 0 can 
be taken. The disordered-to-uniform-ordered phase transition reduces 
therefore to the usual 2nd order Ising phase transition on the underlying 
geometry. 
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2. We compute the dimension of the underlying by computing the critical 
exponents of the disordered-to-uniform-ordered phase transition which, 
by universality, take specific values in each dimension. 

3. We compute the Laplaeian by computing the free behavior of the propa¬ 
gator. This is done explicitly by computing the eigenvalues distribution 
of the matrix M in the free regime, small values of C, and comparing 
with the Wigner semicircle law behavior which must hold with a specific 
radius depending crucially on the kinetic term. 


6.3 Explicit Example: The Fuzzy Sphere 


6.3.1 Phase Diagram 


We consider as an example the multitrace matrix model of 187 which 


comes with the v values u 2 ,i = —1 , fyi = 3/2 , n 2) 2 = 0. The action is given 
explicitly by 

V = BTrM 2 + CTrM 4 + D [TrM 2 ] 2 + B'(TrM ) 2 + C'TrMTrM^. 6) 

The parameters D, B' and C' are constrained as D = 31V/4, B = y/N /2 
and C = —N. The phase diagram of this model is computed by means of 


Monte Carlo elsewhere. The result is shown on figure (6.2). The details of the 


corresponding non-trivial lengthy Monte Carlo calculation will be reported 
elsewhere. As desired we have three stables phases in this particular model 
meeting at a triple point. In other words, we have established that this 
multitrace matrix model sustains the uniform ordered phase which is the 
first requirement. 


6.3.2 Dimension from Critical Exponents 

The uniform ordered phase is also called the Ising phase precisely be¬ 
cause we believe that the corresponding transition to the disordered phase 
is characterized by the universal critical exponents of the Ising model in two 
dimensions derived from the Onsager solution. These critical exponents are 
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defined as usual by the following behavior 

m/N =< \TrM\ > /N ~ (B c - Bf ~ N~ 0/u 
C v /N 2 ~ (B - B c )~ a ~ N a/u 

X =< \TrM \ 2 > - < \TtM\ > 2 ~ (H - R c )" 7 ~ Y 7/v ~ Y 2 " 7 ' 

£ ~ \B - B c \~ v ~ N. (6.7) 


There are in total six critical exponents, the above five plus the critical 
exponent 5 which controls the equation of state, but only two are truly inde¬ 
pendent because of the so-called scaling laws. The Onsager solution of the 
Ising model in two dimensions gives the following celebrated values 


198 


v — 1 , (3 = 1/8 , 7 = 7/4 , a = 0 , rj — 1/4 , 5 = 15. 


( 6 . 8 ) 


This fundamental result is very delicate to check explicitly in the Monte 
Carlo data. Since we must necessarily deal with the critical region we must 
face the two famous problems of finite size effects and critical slowing down. 
In this particular problem, the critical slowing down problem can be shown 
to start appearing in Monte Carlo simulations around N > 60 so we will 
keep below this value and employ very large statistics of the order of 2 20 
to avoid it. A more systematic solution to this problem is to employ the 


Wolf algorithm 199 which we do not attempt here. We simply employ here 
the ordinary Metropolis algorithm. The problem of finite size effects is also 
very serious for the measurement of the critical exponents since the above 
behavior (6.7) is supposed to hold only for large N. This problem can be 


avoided by not including values of N less than 20 and thus below we will 
quote for completeness N = 10 and N — 15 data but, in most cases, we will 
not take them into account in the fitting. 

Since the Ising model appears from the <h 4 theory for large values of 
the quartic coupling it is preferable to use values of C as large as possi¬ 
ble. However, we are limited from above by the appearance of the different 
physics of the transition between the disordered and non-uniform-ordered 
phases around C = 1.5. Thus, we choose C = 1.0 which is relatively large 
but well established to be within the Ising transition with an extrapolated 
critical point around B = —3.07 (see below). The critical behavior of the 
magnetization, susceptibility an d sp ecific heat around the criti cal value of 
B = —3.10 is shown on figure (6.4). We attach in table (6.1) some data 
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relevant for the computation of the critical exponents u, f3, 7 and a. The 
other critical exponents can be determined via scaling laws. 

The measurements of the critical exponents u, (3, 7 and a proceeds as 
follows: 


Critical Point and The Critical Exponent v\ By plotting the 
critical point B c obtained for each N versus N (first and second columns 


of table (6.1)) we get immediately both the N — 00 critical point and 


the critical exponent v. We obtain (see figure (6.7)) 


B c = -1.061(168).Y-°' 926(83) - 3.074(6) =>• , v = 0.926(83(6.9) 
Also we obtain 


B * = -3.074(6). (6.10) 

This prediction for v agrees reasonably well with the Onsager calcula¬ 
tion. In the following we will assume for simplicity that v — 1. The 
above fit is the only instance in which we have included N — 10 and 
N — 15 and thus we believe that the obtained value of — P* is an 
underestimation of the true critical point. 

• Magnetization and The Critical Exponent (3: The magnetization 
and the zero power are defined by 

m =< \TrM\ > , x =< \TrM\ 2 > — < \TrM\ > 2 . (6.11) 


P 0 =< TrM ) 2 > . (6.12) 

Measurements of the magnetization m/N were performed near the ex¬ 
trapolated critical point B = —3.07 for C = 1.0 but inside the uniform 
ordered phase. These are then used to compute the critical exponent 
j3 by searching for a power law behavior. 

More precisely, we measure ln(m/IV) versus In N for each value of B 
very near and around B = —3.10, fit to a straight line in the range 
20 < N < 60 and compute the slope (3, then search for the flattest line, 
i.e. the smallest slope (3. This value marks the transition from the Ising 
phase to the disordered phase. Deeep inside the Ising phase the slope 
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should approach the mean field value —1/4 which can be shown from 
the scaling behavior of the dominant configuration. After determining 
the critical value we then consider the value of B nearest to it but within 
the Ising phase and take the slope there to be the value of the critical 
exponent j3. In our example here, the flattest line occurs at B — —3.13 
with slope —0.088(10) after which the slope becomes —0.109(11) at 
B = —3.14. The slope goes f ast t o the mean held value —0.25 as we 
keep decreasing B. See figure ( |6.5[ ). Our measured value of the critical 
point P* from the magnetization and of the critical exponent f3 are 
therefore 

P* = -3.13. (6.13) 


In — = -0.109(11). In N- 1.423(43) =► 0 = -0.109(11). (6.14) 


• Susceptibility and Zero Power and The Critical Exponent 7 : 

The measurement of the critical exponent 7 is quite delicate and will 
be done indirectly as follows. We rewrite the susceptibility in terms of 
the zero power and magnetization as 


y = < \TrM\ 2 > - < \TrM\ > 2 

= N 2 P 0 — m 2 . (6.15) 

The critical exponent 7 in terms of the critical exponent 'y' of Po is 
then given by 


7 = 2 + 7 . 


(6.16) 


By using the results shown on table (6.1) at B = —3.14, plotted on 
figure (6.7), we obtain the following exponents 

lnP 0 = -0.352(10). In N — 2.289(36) =► 7 = -0.352(10). (6.17) 

Or equivalently 

In N 2 P 0 = 1.648(10). In N — 2.289(36) =^7 = 1.648(10). (6.18) 


For consistency we can check that the second term in the susceptibility 
behaves using the result (6.14) as 


In m 2 = 1.782(22). In N - 2.846(86) =► 7 = 1.782(22). (6.19) 
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Our two measurements of the critical exponent 7 agree reasonably well 
with the Onsager values. 


If we try to fit the values of the susceptibility at its maximum shown in 
third column of table ( 6 . 1 ), i.e. at the peak which keeps slowly moving 
with B , then we will obtain a very bad underestimate of the critical 
exponent 7 given by 


InXmax = 0.515(08). In N - 0.652(30) =* 7 = 0.515(08). (6.20) 


This in our mind is due in part to the dependence of B c on N and 
in another part is an indication of the critical slowing down problem 
showing up in the measurement of this second moment, i.e. the size of 
the fluctuations is observed to grow with N at the critical point but 
not at the correct rate indicated by the independent measurements of 


the zero moment and the magnetization. See figure (6.7). 


Specific Heat and The Critical Exponent a: The sepcihc heat is 
defined by 


c v =< S 2 > - < S >' 


( 6 . 21 ) 


The critical point R* as measured from the specific heat is identified 
by the intersection point of the various curves with different N shown 


on figure (6.4). We get 


R* = -3.08. 


( 6 . 22 ) 


This measurement is contrasted very favorably with the independent 
measurement obtained from the extrapolated value of B c shown in 


equation ( 6 . 10 ) but should also be contrasted with the measurement 


obtained from the magnetization shown in equation (6.13). 


By using the results shown on table (6.1) at the critical point B = 
—3.08, plotted on figure (6.7), we obtain the following exponent 


C 

In -A = 0.024(9). In N - 0.623(31) =* a = 0.024(9). (6.23) 

N 1 
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N 

P C ,P* = -3.07 

Xc 

(a)*,P* = -3.08 

B < P* = -3.13 

m<* 

10 3 (P 0 )<* 

10 

-3.20 

1.704(2) 

56.467(94) 

-3.14 

2.1776(12) 

6.256(6) 

15 

-3.16 

2.089(2) 

129.111(217) 

-3.14 

2.7750(14) 

4.315(4) 

20 

-3.14 

2.436(3) 

229.861(389) 

-3.14 

3.4423(15) 

3.571(2) 

25 

-3.13 

2.716(3) 

365.183(621) 

-3.14 

4.1759(16) 

3.220(2) 

30 

-3.12 

3.017(4) 

524.253(891) 

-3.14 

4.9772(16) 

3.042(2) 

36 

-3.11 

3.283(4) 

749.099(1267) 

-3.14 

5.8878(15) 

2.860(1) 

40 

-3.11 

3.515(4) 

941.139(1607) 

-3.14 

6.5134(14) 

2.782(1) 

50 

-3.10 

3.864(4) 

1461.597(2479) 

-3.14 

7.9250(12) 

2.576(1) 

60 

-3.10 

4.301(5) 

2144.929(3658) 

-3.14 

9.2021(11) 

2.388(1) 


Tabic 6.1: Measurements of the magnetization ( m/N )<*, the susceptibil¬ 
ity x<*, via the zero power (Po)<*, and the specific heat ( C v /N 2 )* used to 
compute the critical exponents /3, 7 and a respectively. Here C = 1.0, the 
extrapolated critical point is B = —3.07, the critical point as intersection 
point of curves of specific heat is B — —3.08, and the critical point as the 
flattest line of decrease of magnetization is B = —3.13. 


6.3.3 Free Propagator from Wigner Semicircle Law 


We can also measure the emergent geometry by measuring the free prop¬ 
agator of the theory. This will give us information on both the dimension 
and the metric since the free propagator is the inverse of the Laplacian A 
which fully encodes the underlying geometry in the sense of 170,171 . This 
goes as follows 


127 


A noncommutative phi-four on a d —dimensional noncommutative Eu¬ 
clidean spacetime Rg reads in position representation 


S = J d d x{^di<&di§ + ^m 2 $ 2 + ^<l>*). (6.24) 

The first step is to regularize this theory in terms of a finite AT —dimensional 
matrix $ and rewrite the theory in matrix representation. Then we diagonal¬ 
ize the matrix <I>. The measure becomes J7L d < hiA 2 ( < h) f dU where are 
the eigenvalues, A 2 ($) = .($j — Qj) 2 is the Vandermonde determinant 

and dU is the Haar measure. The effective probability distribution of the 
eigenvalues can be determined uniquely from the behavior of the expec¬ 
tation values < f d d 'x ( l >2n (x) >. These objects clearly depend only on the 
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eigenvalues and are computed using a sharp UV cutoff A. If we are only 
interested in the eigenvalues of the scalar matrix $ then the free theory A = 0 


can be replaced by the effective matrix model 127 


S = ^Tr<L 2 . 


CYr 


(6.25) 


This result can be traced to the fact that planar diagrams dominates over 
the non-planar ones in the limit A —> oo. This means in particular that the 
eigenvalues <l>j are distributed according to the famous Wigner semi-circle 
law with «o being the largest eigenvalue, viz 

2 


pit) =—2\/ a o - t 2 , ~a 0 <t<+a 0 . 
7rttn v 


(6.26) 


In the most important cases of d = 2 and d = 4 dimensions we have explicitly 


A 2 


«oK A ) = 4^2 ( A ' “ ln ( X + ^2 )) 1 d = 4 - 


(6.27) 


1 A 2 

Oq (m, A) = - ln(l H--) , d — 2. 


7T 


rrr 


(6.28) 


Obviously, dimension four is eliminated by the results of the critical ex¬ 
ponents. In two dimensions the regulator A originates in only one of two 
possible noncommutative spaces 


127 


1. Fuzzy Torus: As it turns the results on the fuzzy torus are different 
from those obtained using a sharp momentum cutoff due to the different 
behavior of the propagator for large momenta and as a consequence the 
resulting formula for ccq is different from the above equation (6.28). We 
obtain instead 

d 2, r 1 


«o(m, A) = 4 


(27t) 2 JA(1 — cos t i) + m 2 l 2 / 2 


d = 2. (6.29) 


l here is the lattice spacing, the noncommutativity is quantized as 6 = 
Nl 2 /n and the cutoff is 


A = 




(6.30) 


The above behavior can be easily excluded in our Monte Carlo data 
and by hindsight we know that this should be indeed so because the 
original multitrace approximation is relevant to the fuzzy sphere. 
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2. Fuzzy Sphere: The fuzzy sphere S“+ = CPjy is the simplest of fuzzy 
projective spaces CP^. In this case Af = N + 1 and the scalar held <f> 
becomes an N x N matrix </> given by (f> — \J2 tt/N a<$>. In this case the 
cutoff is given in terms of the matrix size N and the radius R of the 
sphere by 


A = 


N 

R' 


(6.31) 


Also, in this case the mass parameters B and m 2 are related by 


m 


2 


b 

MT 2 ' 


(6.32) 


By using B = B/N 3 ^ 2 and choosing a 
obtain 

A 2 _ 2t r 
m 2 y/NB 


We get then 


27 t/N, so that $ 


ckq (m, A) 


- ln(l + 
7 r 


2tt 

Vnb 


)• 


4 >, we 


(6.33) 


(6.34) 


In the limit B —> oo we get the one-cut 5 2 = 2 N/B of the Gaussian 
matrix model BTrM 2 , viz B = 2J\f/a^. This can also be obtained by 
taking the limit B —> oo of the one-cut (deformed Wigner semicircle 
law) solution 


P(A) 


1 

Ntt 


(2CA 2 + B + C6 2 )VS 2 - X 2 


5 2 


1 

3 C 


{-B + VB 2 + 12 NC) 
(6.35) 


of the quadratic matrix model BTrM 2 + CTrM 4 . 
This result was also generalized in 


191 . The eigenvalues distribution 


of a free scalar held theory on the fuzzy sphere with an arbitrary kinetic 
term, viz S = Tr(MK.M + BM 2 ) /2, where /C(0) = 0 and fC is diagonal 
in the basis of polarization tensors T^, is always given by a Wigner 
semicircle law with a radius 


R 2 = 6 2 = 


«5 — 


4 m 

N 


N -1 


. m = X 


1=0 


21 + 1 
K{l) + B' 


(6.36) 


179 












6.4. GENERALIZATION AND CONCLUSION 


Some Monte Carlo results are shown on figures (6.9) and (6.10). These are 
obtained in Monte Carlo runs with 2 20 thermalization steps and 2 18 ther- 
malized configurations where each two configurations are separated by 2 4 
Monte Carlo steps in order to reduce auto-correlation effects. We consider 
IV = 20 — 40, C = 0.05 - 0.35 and 5 = 0-5. 

It is not difficult to convince ourselves that the mass parameter B is 
precisely the mass squared in this regime. For each value of ( N,C,B) we 
compute the eigenvalues distribution p( A) and fit it to the Wigner semicircle 
law ( 6.26[ ) (see figure ( |6.9P ). We obtain thus a measurement of the radius 
of the Wigner semicircle law 5 2 = «q = K 2 . We have checked carefully 
that in this regime the Wigner semicircle law is the appropriate behavior 
rather than the one-cut solution (6.35) as evidenced by the first graph in 
figure (p.9). The measurement of the radii 5 2 for various values of B is then 


plotted and compared with the expected theoretical behaviors (6.34) as well 
as wit h the B —y oo behavior 5 2 = 2N/B (see figure (6.10)). The agreement 
with (6.34) is very reasonable with some deviation for small values of B as 


we approach the non-perturbative region where the uniform ordered phase 
appears at some B < 0. This discrepancy for small values of B is already 
seen on figure (6.9) when we fit the distributions to the Wigner semicircle 
law. However, this effect is reduced as we decrease the value of C. 

In summary we conclude that we are indeed dealing with the geometry 
of the fuzzy sphere and, given hindsight, we know that this should be true. 


6.4 Generalization and Conclusion 

The emergence of geometry in the very early universe is a problem of 
fundamental importance to our understanding of quantum gravity and cos¬ 
mology. In this letter, we have proposed a novel scenario for the emergence of 
geometry in random multitrace matrix models which depend on a single her- 
mitian matrix M with full unitary U(N) invariance and without any kinetic 
term. Thus, the model under consideration has no geometry a priori pre¬ 
cisely because of the absence of a kinetic term. On the other hand, previous 
proposals of emergent geometry required the input of several matrices with 
some rotational symmetry group besides the U(N) gauge symmetry 

Our proposal consists in checking whether or not the uniform ordered 
phase is sustained by the multitrace matrix model under consideration. If 
yes, then the dimension of the underlying geometry, in the region of the phase 
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6.4. GENERALIZATION AND CONCLUSION 


diagram where the uniform ordered phase is stable, can be inferred from 
the values of the critical exponents of the Ising phase transition. Whereas, 
the metric/Laplacian of this geometry can be inferred from the behavior of 
the free propagator encoded in the Wigner semicircle law behavior of the 
eigenvalues distribution of the matrix M in the weakly coupled regime. An 
explicit example is given in which the geometry of the fuzzy sphere emerges, 
with all the correct properties, in the phase diagram of a particular multitrace 
matrix model containing multitrace terms depending on the moments m\ = 
TrM, m 2 = TrM 2 and 7773 = TrM 3 in a particular way 

This idea can be generalized in a straightforward way to all higher fuzzy 
projective spaces CP n and fuzzy tori T' 1 by tuning appropriately the coef¬ 
ficients of the multitrace matrix model and/or including higher moments in 
the multitrace matrix model . 
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Figure 6.1: The phase diagram of noncommutative phi-four theory on the 
fuzzy sphere. In the first figure the fits are reproduced from actual Monte 
Carlo data 186 . Second figure reproduced from 180 with the gracious 


permission of D. O’Connor. 
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Model l,N=10-50 



Figure 6.2: The phase diagram of the multitrace matrix model of 187 . The 


Ising and matrix transition data points are not shown but we only indicate 
their extrapolated fits whereas the N = 25, N = 36 and N — 50 stripe data 
points are included explicitly. 
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Figure 6.3: The Ising critical behavior. 
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Figure 6.4: The Ising critical behavior. 
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Figure 6.5: The critical exponent /3. 
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Figure 6 . 6 : The critical exponents u, 7 and a. 
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Figure 6.7: The critical exponents is, 7 and a. 
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Figure 6.8: The semicircle law as a function of B. 
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Figure 6.9: The semicircle law as a function of B. 
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Figure 6.10: The behavior of the radius of the Wigner semicircle law 5 2 = a jj 
as a function of B. 











Chapter 7 

Gauge Theory on The 
Noncommutative Torus 


The study of massless field theories on a torus is of great interest in the 
noncommutative case because the compactness of the spacetime gives a nat¬ 
ural infrared regularization of the theory. One may therefore analyse more 
carefully the ultraviolet behaviour and also the new light degrees of freedom 
which are responsible for the UV/IR mixing. From a more mathematical 
point of view, the noncommutative torus constitutes one of the original ex¬ 
amples in noncommutative geometry 200-203 which captures the essential 
topological changes which occur when one deforms a compact space. It is per¬ 
haps the most basic example which still contains a rich geometrical structure. 
In this section we shall describe some basic aspects of the noncommutative 
torus with particular emphasis on the properties of vector bundles defined 
over them. From the study of the global properties of gauge theories defined 
on this space. 


7.1 The Noncommutative Torus 

Most of what we have said about noncommutative quantum field theory 
is true when 1Z D is replaced by a D dimensional torus T D , with only subtle 
changes that we shall now explain. Let be the D x D period matrix of 
"1° which is a vielbein for its metric, i.e.J^ S ab . Here and in the 

following the indices i,j ,... will label spacetime directions while a,b ,... will 
denote indices in the frame bundle of T l> . The matrices parametrize 
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the moduli of D dimensional tori and they may be regarded as maps from 
the frame bundle to the tangent bundle of T L) . They define the periods of 
the directions of T D , 


i 

x 1 ~ x l + , a — l,...,D 

a 


(7.1) 


for each i = 1 When is not proportional to Sf , the identifica¬ 

tions (7.1) for a 7 ^ i describe how the torus is tilted in its parallelogram 
representation. 

Smooth functions on the torus must be single-valued, which implies that 
the corre- sponding Fourier momenta k are quantized as 


h 



m a G Z 


(7.2) 


Therefore, to describe the deformation of the function algebra, one cannot 
use the unbounded operators x l obeying ( 1 . 1 ). 


[x % , x j ] = id lj 


(7.3) 


Instead, one must restrict to the proper subalgebra of the algebra of 
noncommutative 1Z D that is generated by the Weylbasis of unitary operators 

Z a = (7.4) 

which generate the algebra 

Z a Z b = e~ 2nieab Z b Z a , (7.5) 

where 

e- = 2»(£ )V(f )‘ U6) 


are the corresponding dimensionless noncommutativity parameters. The 
commutation relations (7.5) define the algebra of functions on the noncom¬ 
mutative torus. Formally, if L = Z D is the lattice of rankD (with bilinear 
form Gij ) which generates the torus as the quotient space T 1 " 1 = 1Z D /C, 
then the projective regular representations £© in (7.5) of the lattice group C 
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are labelled by an element Q ab of the second Hochschild cohomology group 
H 2 (C , U( 1)) This latter characterization can be generalized to describe other 
sorts of noncommutative compactihcations of 1Z D 

Any function on T D can be expanded as a Fourier series 

f{x) = Y, frne 2 ^~^ maxi (7.7) 

m£Z D 


204 


The corresponding Weyl algebra is generated by the operators () and Weyl 
quantization takes the form of the map 

W[f] = f d D xf(x) A(x), (7.8) 

where the integration is taken over T 1 " 1 and 


A(x) = 


| det AA 


D 

E 

rh^Z D cl—1 


nW n 




b m be ~2ni(Y,~ 


L )?m a x l 


(7.9) 


a<b 


is a periodic field operator, 


i 

A (x + i) — A(o;), a — l,...,D, (7-10) 

a 


with i a unit vector in the i-th direction of spacetime. Like on 1Z V , we 
may introduce anti-Hermitian, commuting linear derivations <9j which on the 
noncommutative torus are defined by their actions on the Weyl basis, 


[di,Z a ] = Z a - 


The basis (7.9) then has the requisite property 

[4, A(i)l = -d,A(x) 


(7.11) 


(7.12) 
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7.2 Topological Quantum Numbers 

A U(N) noncommutative Yang-Mills theory on the torus T D can be con¬ 
structed in much the same way as we did in the previous section. If we restrict 
to gauge field configurations which are single-valued functions on T D , then 
everything we have said goes through with- out a hitch, with single-valued 
star-unitary functions g(x) parametrizing the star-gauge transformations 

A(x) -> g(x) * Ai(x) *g{x) ] - ig(x) * <9^(Y) t (7.13) 


The only difference which arises is that, like in the commutative case, there 
are extra observables associated with the non-trivial homotopy of the torus. 


The most general star-gauge invariant observable is still given by 67 


0(C V ) 


d D x tr N (U(X-,C v ))*e ik( ‘ v)xi , 


(7.14) 


but now there is a larger set of line momenta. Because the momenta are now 


quantized as in (7.2), the identification of the translation vector v in 


Jki(v)x l 


-k g(x ) -k e 


—iki(v)x l _ 


= g{x + v), 


(7.15) 


is ambiguous up to an integer translation of the periods of T D , and the 
relationship 


ki{v) = (9 Q ijV 3 , 


is now modified to 


l 

v l 9 tj kj(v, n ) + ^ n a 

a 


(7.16) 


(7.17) 


for arbitrary integer-valued vectors n a . When 9 = 0, the relationship (7.17) 
reproduces the well-known result that the only open line observables in ordi¬ 
nary Yang-Mills theory are those which are associated with loops that wind 
n a times around the a-th non-contractible cycle of the torus. Therefore, we 
obtain the analog of Polyakov lines in noncommutative Yang-Mills theory 


associated with the different homotopy classes of the torus 67 


More interesting things happen, however, when we consider gauge held 
configurations of non-vanishing topological charge on the noncommutative 
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torus. An elegant way to keep track of the quantum numbers associated 
with topologically non-trivial gauge fields is through their Chern numbers. 
In the commutative case, these would be represented by the integers 

mn)= (ftTNF n /(27f) n (7.18) 

defined in terms of the curvature two-form F of some gauge connection of a 
U (. N ) gauge bundle E over T D , and suitably integrated over cycles of the 
torus. For n — 0 they produce the rank N of the vector bundle E, for n — 1 
they yield the fluxes Q ab of the gauge fields through the surface formed 
by the a-th and 6-th cycles of T 1 ’ , and for n — 2 they give the instanton 
number k of the bundle E when D — 4. We can collect these integers into 
the inhomogeneous Grassmann form 

d, ^ 

Ch„(£) = JV + Y. (")(S)a 1 ...a 2n p“‘...p M " (7-19) 

n\ 

n= 1 

where here and in the following we will assume that the spacetime torus 
has even dimension D = 2d. We have introduced a set p a , a = 1 of 

anticommuting Grassmann variables, 

P a p b = ~p b p a , (7.20) 

which can be thought of as local generators of the cotangent bundle of . 
The quantity cho(E) then defines an integer cohomology class of the ordinary 
torus T d . Given these integers which characterize the given bundle E, there 
is an elegant formula for the noncommutative Chern character 

yv\F] 

cho,{E ) = Tr (g) tr^exp—— (7-21) 

27T 

which characterizes the corresponding gauge bundle over the noncommuta¬ 
tive torus. Here F is the noncommutative curvature two-form of the bundle 
with local components 


F a b — ^ f^ y j 

a b 


(7.22) 
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where F tJ is defined by 

Fij = diAj — djAi — i(Ai-k Aj — Aj * Ai) 

= djAj - djAj - i[A h A j] 

+ V(<9 fc A<M,) - dkAjdtAi) + 0(6 2 ) (7.23) 


It can be regarded as an element of the ordinary cohomology ring H even (A D , 77) 


of even degree differential forms on the torus. The quantity (7.21) can be 


written in terms of (7.19) through the Elliott formula 205 


c M £) = ex P ( Qe-AA 


ch 0 (E ) 


(7.24) 


with O regarded as a two-cycle of the homology group H 2 (T D ,1Z) (206,207 


The coefficients of p ai ...p a2n in the expansion of (7.24) define the n — th 


noncommutative Chern numbers of the given noncommutative gauge the¬ 
ory. They represent the topological invariants of the corresponding deforma¬ 
tion E —> Eq from a commutative to a noncommutative gauge bundle. In 
the commutative limit 0 = 0,cho(E) generates the ordinary integer-valued 
Chern numbers. But for 0^0 they are non-integral in general. 

For example, in two dimensions we find 


ch e (E) = (N ~ Q0) + Qp 1 p 2 , 


where Q is the magnetic flux through T“ 


(7.25) 
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Conclusion 


In this article we have studied IKKT Yang-Mills matrix models with mass 
deformations in three and six dimensions. 

The 3—dimensional IKKT matrix models considered here are very sim¬ 
ilar to the ones studied in 135,138 . However, the dynamically emergent 
geometry, which is given by a fuzzy two-sphere S 2 N , is found to be stable for 
all values of the deformation parameter M. This was anticipated previously 
for M = 1/2 in 141 . In deed, the critical gauge coupling constant a is found 
to scale as in equation (5.24), i.e. as a ~ 1/y/N. The sphere-to-matrix 
transition line is pushed to 0 and only one phase survives. 

In this case the fuzzy sphere acts then as a regulator of the commutative 
sphere, and as a consequence, fuzzy held theory and fuzzy physics, based 
on this emergent fuzzy sphere, makes full sense for all values of the gauge 
coupling constant. 

We have also studied in this article 6—dimensional IKKT matrix models, 
with global SO(3) x SO( 3) symmetry, containing at most quartic powers 
of the matrices proposed in ([146 . The value M = 1/2 of the deformation 
corresponds to the model of |147|. This theory exhibits a phase transition 
from a geometrical phase at low temperature, given by a fuzzy four-sphere 
S 2 n x S 2 N background, to a Yang-Mills matrix phase with no background 
geometrical structure at high temperature. 

The geometry as well as an Abelian gauge held and two scalar helds 
are determined dynamically as the temperature is decreases and the fuzzy 
four-sphere condenses. 

The transition is exotic in the sense that we observe, for small values of M, 
a discontinuous jump in the entropy, characteristic of a 1st order transition, 
yet with divergent critical huctuations and a divergent specihc heat with 
critical exponent a = 1/2. The critical temperature is pushed upwards as 
the scalar held mass is increased. For small M, the system in the Yang-Mills 
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phase is well approximated by 6 decoupled matrices with a joint eigenvalue 
distribution which is uniform inside a ball in R 6 . This gives what we call 
the d = 6 law given by equation (5.47). For large M, the transition from the 
four-sphere phase to the Yang-Mills matrix phase turns into a crossover and 
the eigenvalue distribution in the Yang-Mills matrix phase changes from the 
d = 6 law to a uniform distribution. 

In the Yang-Mills matrix phase the specific heat is equal to 3/2 which 
coincides with the specific heat of 6 independent matrix models with quartic 
potential in the high temperature limit and is therefore consistent with this 
interpretation. Once the geometrical phase is well established the specific 
heat takes the value 5/2 with the gauge field contributing 3/2 and the two 
scalar fields contributing lQ 

This should be contrasted with the case of S 1 2 N in which the specific heat in 
the Yang-Mills matrix phase is equal to 3/4, coinciding with the specific heat 
of 3 independent matrix models with quartic potential in the high tempera¬ 
ture limit, while in the geometrical phase the specific heat takes the value 1 
divided equally, i.e. with the gauge field contributing 1/2 and the two scalar 
fields contributing 1/2. 

The counting on S 2 N is clear cut. In the sphere phase which coincides with 
the perturbative region of the theory the gauge field contributes 1/2 [167 
and the scalar field contribute 1/2 (trivial to check in the quartic matrix 
model for very large positive values of the mass parameter 168 169]). 


The counting on x 


C2 

°N 


is more involved. There are here two scalar 


fields and one gauge field with 4 components. Again, it is very natural to 
suppose that each scalar field will contribute 1/2 (since they are free and 
dimension does not enter in the quartic matrix model). Therefore, the gauge 
held will contribute 3/2. This is the picture we also get, at least formally, 
by assuming that the gauge held is free and Abelian (which is true deep in 
the sphere phase in the large N limit) and then fixing the gauge in the axial 
gauge. 

The 6—dimensional IKKT Yang-Mills matrix models studied here present 
thus an appealing picture of a 4—dimensional geometrical phase emerging as 
the system cools and suggests a scenario for the emergence of geometry in 
the early universe. See 


148 and references therein. 


1 Recall that in the 3d Yang-Mills matrix model the specific heat takes the value 1 in 

the geometrical phase which is attributed in this case to the normal scalar field since there 
is no propagating gauge degrees of freedom in 2 dimensions. 
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The model presents thus an appealing picture of a geometrical phase 
emerging as the system cools and suggests a scenario for the emergence of 
geometry in the early universe. 
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Appendix A 


A.l Susceptibility and Specific Heat 


A. 1.1 Susceptibility 

We consider <t> 4 on the fuzzy sphere coupled to a constant magnetic held 
H given by the action 


5 = Tr(a$[L a , [L a , $]]) + 6d> 2 + c$ 4 + H$ 
The magnetization and the susceptibility are defined by 


(A.l) 


magnetization = — < TV$ > 


1 d 

NdH 


In Z 


(A.2) 


susceptibility = < (Tr<f>) 2 >< Tr ( l> > 2 
d 2 


dH 2 


In Z 
d 


On the fuzzy sphere we have 


= —A-—magnetization 
oH 


N 


2 R T r , iO fl 2 R 2 

x a = ^L a , [x a ,x b ] = —e abc x c , 9—, Tr = ^ R2 

The regularized noncommutative plane is then defined by 


d 2 x. 


(A.3) 


(A.4) 


1 1 f 

x 3 = R, [xi,x 2 \ = id, di = ——e.ijLj = --e^Xj, / d 2 x = 2n6Tr. (A.5) 
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We have e V2 = 1. The above action becomes, including a rescaling of the held 
—> (j) = \jNa/2n$ given by the equation 


S = 2ir6Tr(^(j)didi(j) + ^m 2 0 2 + * Ao 1 + h<j)). 


(A.6) 


m 2 = 


o A = 


47TC 


h = 


N H 


(A.7) 


aR 2 ’ Na 2 R 2 ' 1 V 2na2r 2 

The commutative limit is 6 —y 0. By using a lattice in this limit we have 


S = l 2 lattice + \m 2 (fa + jXc/)* + h(j) n ). 


We compute in this limit on the lattice 


magnetization = — < Trp > 

a r 2 i 2 i 


47T.R 2 W 2 


> 


(A.8) 


(A.9) 


A. 1.2 Specific Heat 

The specific heat is defined by 


c v = 


d 2 


In z 


dp 2 

= <S 2 > - <S> 2 . 

The inverse temperature is introduced in the usual way as 

Z = f dMexp(—pS[M]). 


(A.10) 


(A.ll) 


The calculation of the effective potential proceeds as before with the replace¬ 
ment a —» a/3. The partition function in the quartic multitrace approximation 


is 


Z = 


dAA 2 (A)exp(-/3U 0 ) + p / dAA 2 (A)exp(- / dU 0 )(-U 2 )(A.12) 


+ p 2 / dAA 2 (A)exp(-/3U 0 )(-U 4 +-U 2 2 ) 


(A.13) 
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A straightforward calculation yields 

C v =< (V + V 4 ) 2 > - < (V + V 4 ) > 2 -2 < (V 4 + 2 Vl + 2U 2 U 4 ) 2 > .(A. 14) 
The last term could make this approximation of the specific heat negative. 


This actually happens in the approximation of 208 
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A.2 Grosse-Wulkenhaar Model 

The multitrace approach can also be applied to a regularized noncom- 
mutative $2 on the Moyal-Weyl plane in the matrix basis 193 with action 
given by 


S = Tr 


N 


1 


U 


-m M + —M 4 + a EM + Jui 


_2 ' N 

Two cases are of importance to us here: 
1 


6 






(A.15) 


The noncommutative theory without a harmonic oscillator term. In 
this case the effective action takes the form 

S'effe = bTr N M' 2 + cTr N M 4 + d(Tr N M 2 ) 2 + b x ( Tr N M ) 2 + Cl ( Tr N M) A 
+ diTr]yM 2 (TrNM ) 2 + eTrjvAfTrjvM 3 . 

(A.16) 


The parameters are given by 
' 2 cLV 


m atv u 

b= ~ + — ’ C= N 


a 2 N 


a a 

h = ~2 ’ Cl = 24N 2 


24 
, d 1 = - 


, d = - — 

’ 12 

a 2 a 2 

12A' ’ e = IT' 


(A. 17) 


If we assume the symmetry M —>■ — M then all odd moments vanish 
identically and we end up with the action 


^ = bTr N M 2 + cTr N M A + d{Tr N M) . 


(A. 18) 


2. At the self-dual point we have D 2 = 1, and thus y/cu = 0, and as a 
consequence the effective action reduces to the multitrace model 


S'effe = bTr N M 2 + cTr N M A + d(Tr N M 2 ) 2 . 
The parameters b, c and d are given by 


, m 2 aN u a 2 N 

b = -1-, c =- 

2 2 5 N 24 


d=*. 

24 


(A.19) 


(A.20) 


Both the actions (A. 18) and (A.19) do not contain odd moments and thus 
the corresponding phase diagrams are expected to not contain the uniform 
ordered phase with all matrix-like behavior as consequence. 
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Appendix B 


B.l Metropolis Algorithm for Yang-Mills Ma¬ 
trix Models 


B.1.1 Metropolis Accept/Reject Step: 

The basic Yang-Mills action of interest is 

Sym = ~Tr[X^X u f 

d d 

= E (w. -Y'- 2 >- 

fl= 1 i/=fl +1 


We perform the variation 

y a x' x = X A + AX a 


where 


(X\)nm d8 n i8mj+ d 8 n j8 ml 


The corresponding variation of the action is 


(B.l) 


(B.2) 


(B.3) 


AS = S(X') - S(X) (B.4) 

The Metropolis accept/reject step is based on the probability distribution 

P[X\ = min(l,exp(— AS)) (B.5) 
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B.l. METROPOLIS ALGORITHM FOR YANG-MILLS MATRIX 

MODELS 


B.1.2 Auto-Correlation Time: 

In any given ergodic process we obtain a sequence (Markov chain) of 
field/matrix configurations 4>i,4>2, ■ We will assume that <pi are ther- 

malized configurations. Let f some (primary) observable with values /,; = 
f((j>i) in the configurations <pi respectively. The average value < / > of / 
and the statistical error 5f are given by the usual formulas 

1 T 

</> = < B - 6 ) 

i=1 

Sf - ^ (B.7) 

The standard deviation (the variance) is given by 

a 2 =< f 2 > - < / > 2 (B.8) 

The above theoretical estimate of the error is valid provided the thernralized 
configurations 4>i , are statistically uncorrelated, i.e. independent. 

In real simulations, this is certainly not the case. In general, two consec¬ 
utive configurations will be dependent, and the average number of config¬ 
urations which separate two really uncorrelated configurations is called the 
auto-correlation time. The correct estimation of the error must depend on 
the auto-correlation time. 

We define the auto-correlation function R, and the normalized auto¬ 
correlation function pj for the observable / by 

1 T ~ j 

r » = £(/.- < / >)()/.+j+ < / >)■ (B.9) 

i= 1 


Pj = 



(B.10) 


These function vanish if there is no auto-correlation. Obviously To is the 
variance a 2 , viz To = cr 2 . In the generic case, where the auto-correlation 
function is not zero, the statistical error in the average < / > will be given 
by 



Vt 


V^int 


(B.ll) 
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B.l. METROPOLIS ALGORITHM FOR YANG-MILLS MATRIX 

MODELS 


The so-called integrated auto-correlation time Ti nt is given in terms of the 
normalized auto-correlation function pj by 

1 X 

'Tint = 2 + Pj (B.12) 

3=1 

The auto-correlation function Tj , for large j, can not be precisely determined, 
and hence, one must truncate the sum over j in r int at some cut-off M, in 
order to not increase the error 5 t int in r int by simply summing up noise. 
The integrated auto-correlation time r int should then be defined by 

1 M 

T int = - + J^ Pj ( B - 13 ) 

3=1 

The value M is chosen as the first integer between 1 and T such that 


M A 4r; nt + 1 

The error 5T int in T mt is given by 


. AM + 2 

eTlnt \l j-, Tint 


(B.14) 


(B.15) 


This formalism can be generalized to secondary observables F which are 
functions of n primary observables f a , viz F = F(f 1 , f 2 ,..., f n ). 


B.1.3 Errors 

We use the Jacknife method to estimate the errors. Given a set of T = 2 P 
( with P some integer ) data points f(i) we proceed by removing z elements 
from the set in such a way that we end up with n — T/z sets ( or bins). The 
minimum number of data points we can remove is z — 1 and the maximum 
number is z — T — 1. The average of the elements of the ith bin is 

f(J) ~ /((* _ X ) 2 + 4')) , i = 1, n. (B.16) 

3=1 ' 

For a fixed partition given by z the corresponding error is computed as follows 

e(z) = £(< y(j) >, - < f >) 2 . (B.17) 

\ i= 1 j=l 


< y(j) >i= 


T-z 


As 

— i 
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B.2. THE HYBRID MONTE-CARLO ALGORITHM 


We start with z — 1 and we compute the error e(l) then we go to z = 2 and 
compute the error e(2). The true error is the largest value. Then we go to 
z — 3, compute e(3), compare it with the previous error and again retain the 
largest value and so on until we reach z = T — 1. 


B.2 The Hybrid Monte-Carlo Algorithm 


The hybrid monte carlo algorithm is a combination of the molecular dy¬ 
namics method and the metropolis algorithm. First we introduce a fictitious 
time r and 3 bosonic matrices P a = P a (r). These P a play the role of con¬ 
jugate momenta for the 3 matrices X„ = X a (r). We will have a classical 
dynamical system described by (XC a (r), P a (r)) with Hamiltonian 

H = \TrPZ + S, ff [X] = tr r P„ 2 + SB[X]-Tfllog5Ai. (B.18) 
The partition function is given by 

Z= J [dX a ][dP a ]e- H ^ p y (B.19) 


By integrating out P a we obatin the original partition function, viz 


Z= j[dX a ]e~ s ^ x] . 


(B.20) 


The Hamiltonian classical equations of motion are 


d{P a )jj _ dH 
dr d(X a )ij 


dS eS 

d(X a )ij 



[Af,, [X a , X b T 2 iot(i a b c 


x,x, 


TRM- l 


dM 

d(X a ) i3 ' 


ji 

(B.21) 


djXjjj = dH 
dr d(P a )ij 


(B.22) 
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